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BREAKING  ICE  WITH  EXPLOSIVES 


Malcolm  Mellor 


INTRODUCTION 

Although  explosives  have  been  used  to  break 
floating  ice  sheets  for  at  least  200  years,  systematic 
design  procedures  for  ice  blasting  are  still  evolv¬ 
ing.  In  many  cases,  simple  field  tests  can  soon  es¬ 
tablish  optimum  procedures  for  the  prevailing 
conditions  but,  taken  in  isolation,  the  results  of 
such  tests  are  of  limited  value  for  predicting  blast 
effects  under  different  conditions. 

About  10  years  ago,  an  attempt  was  made  to  de¬ 
velop  design  curves  which  could  be  used  for  study¬ 
ing  the  potential  of  ice  blasting  as  an  aid  to  ship 
navigation  in  ice-covered  waters  (Mellor  1972).  All 
available  data  from  field  tests  were  compiled,  and 
the  dependence  of  “crater  radius”  on  charge  size, 
charge  depth,  ice  thickness,  ice  type,  and  explosive 
type  was  considered.  There  were  insufficient  data 
for  complete  consideration  of  all  variables,  and 
some  simplifications  had  to  be  made.  Variations 
of  ice  type  and  explosive  type  were  ignored,  and 
cube  root  scaling  was  used  to  account  for  the  vari¬ 
ation  of  charge  size  when  using  point  charges. 
Scaled  crater  radius  was  related  to  scaled  charge 
depth  and  to  scaled  ice  thickness  by  means  of  mul¬ 
tiple  regression  analysis,  and  design  curves  were 
drawn.  These  curves  (“MM  72  curves”)  proved 
instructive,  and  they  gave  good  predictions  for 
near-optimum  blasting  conditions.  However,  the 
curves  could  not  be  relied  upon  to  predict  behav¬ 
ior  at  large  values  of  scaled  charge  depth,  or  large 
values  of  scaled  ice  thickness,  since  most  of  the 
test  data  used  for  the  analysis  referred  to  blasting 
conditions  close  to  optimum. 


In  1981,  the  MM72  curves  were  tested  against 
the  results  of  major  field  work  by  the  Canadian 
Armed  Forces  (Fonstad  et  al.  1981).  According  to 
comparisons  with  the  original  data,  the  curves 
tended  to  underpredict  scaled  crater  radius  for 
large  values,  and  to  overpredict  for  small  values; 
they  did,  however,  give  good  prediction  for  the 
Canadian  tests.  In  1982  some  further  U.S.  Army 
tests  were  made  on  very  thin  ice,  and  for  the  first 
time  a  few  tests  were  deliberately  made  under  con¬ 
ditions  far  from  optimum. 

In  this  report  the  general  problem  is  reexam¬ 
ined,  and  the  design  curves  are  revised  on  the  basis 
of  currently  available  data. 


GENERAL  BEHAVIOR  OF 
UNDERWATER  EXPLOSIONS 

When  explosives  are  used  to  break  a  semi-infi¬ 
nite  solid  medium,  detonation  of  charges  on  the 
surface  is  very  inefficient.  Experience  suggests 
that  the  same  is  true  of  charges  fired  in  air  on  top 
of  a  floating  ice  sheet.  We  are  therefore  concerned 
largely  with  the  detonation  of  charges  in  water  be¬ 
neath  the  ice  cover. 

When  a  concentrated  charge  of  high  explosive  is 
detonated  well  below  the  surface  of  deep  water,  it 
propagates  a  spherical  shock  wave  and  creates  a 
gas  bubble. 

The  shock  propagates  at  high  velocity  and  its 
pressure  decays  with  distance,  largely  because  of 
spherical  spreading.  From  similitude  considera¬ 
tions,  radial  distances  from  charges  of  different 


Figure  1.  Variation  of  peak  pressure  with  sealed  distance 
for  a  point  charge  of  TNT  in  deep  water.  (From  a  relation 
given  by  Swisdak  11978].) 


size  can  be  scaled  with  respect  to  charge  radius,  or 
with  respect  to  the  cube  root  of  charge  weight  (as¬ 
suming  approximate  constancy  of  charge  density 
and  specific  energy  and,  if  necessary,  taking  into 
account  differences  of  detonation  pressure  and 
bubble  characteristics  for  explosives  of  different 
type).  Spherical  spreading  causes  the  shock  wave 
energy  to  decay  with  the  square  of  the  radius  of 
the  source  and,  since  wave  energy  is  proportional 
to  the  square  of  wave  amplitude,  this  leads  to  ex¬ 
pectation  that  amplitude  might  be  inversely  pro¬ 
portional  to  radius.  There  is,  in  addition,  some 
dissipation  and  dispersion,  and  this  increases  the 
rate  of  pressure  decay  with  distance,  probably  by 
an  exponential  factor.  In  explosions  technology, 
the  variation  of  peak  shock  pressure  with  scaled 
distance  is  usually  expressed  by  an  empirical  equa¬ 
tion  derived  from  log-log  plots  of  test  data,  for 
example  (Swisdak  1978): 


where  Pmn  =  peak  shock  pressure 
P,  =  reference  “pressure” 

R  =  radial  distance  from  the  source 


IV  =  charge  weight 
a  =  an  exponent  of  order  unity. 

Since  W'A  is  proportional  to  charge  radius  R, 
( WVi/R )  can  be  treated  as  a  dimensionless  quanti¬ 
ty  for  a  given  explosive  type.  For  TNT,  P,  =  52.4 
and  a  =  1.13  when  Pmax  is  in  MPa,  W  is  in  kg, 
and  R  is  in  metres  (see  Fig.  1).  The  range  of  applic¬ 
ability  of  eq  1  is  3.4  <  <138  MPa.  The  num¬ 

bers  are  not  much  different  for  other  high-density 
solid  explosives.  Figure  2  gives  a  conversion  of 
(R/W>a)  to  a  true  dimensionless  radius,  taking  ac¬ 
count  of  variations  in  charge  density. 

The  detonation  produces  gas  at  high  tempera¬ 
ture  and  pressure,  and  this  creates  a  bubble  in  the 
water.  The  gas  bubble  expands  against  hydrostatic 
pressure,  but  because  of  inertial  effects  the  expan¬ 
sion  does  not  cease  until  the  bubble  pressure  has 
dropped  well  below  the  external  water  pressure. 
Eventually  it  collapses,  again  with  inertial  over¬ 
run,  until  the  bubble  pressure  is  well  in  excess  of 
water  pressure.  This  process  gives  rise  to  succes¬ 
sive  bubble  pulsations.  While  these  pulsations  are 
occurring,  a  bubble  in  deep  water  rises  by  virtue  of 
its  buoyancy.  However,  proximity  to  surfaces 
complicates  the  translational  motion;  the  bubble 
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Figure  2.  Conversion  of  scaled  distance.  The  vertical  axis  gives 
distances  normalized  with  respect  to  charge  radius;  the  horizontal 
scale  gives  distances  scaled  with  respect  to  the  cube  root  of  charge 
weight.  The  conversion  is  given  for  four  values  of  charge  density. 
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Figure  3.  Spalling  of  the  water  surface  by  an  underwater  explo¬ 
sion.  (After  Young  1973.) 


has  a  tendency  to  move  towards  solid  boundaries 
and  away  from  free  boundaries.  For  present  pur¬ 
poses,  translational  motion  can  be  ignored.  The 
maximum  bubble  radius  /?bm  for  the  first  pulse  in 
deep  water*  can  be  expressed  as 

His.  =  - H -  (2) 

where  A-  is  a  constant  for  a  given  explosive  type,  H 
is  charge  depth  (head  of  water)  and  is  the  at- 

*  An  underwater  explosion  is  usually  considered  to  be  "deep" 
when  the  charge  depth  is  greater  than  A^.  For  present  pur¬ 
poses,  water  can  be  regarded  as  "deep”  when  the  depth  is 
greater  than  4Abn). 


mospheric  pressure  head  ( *  10  m  water).  With  R 
and  H  in  metres  and  W  in  kg,  the  value  of  K  for 
TNT  is  3.50  and  thus 

^bm  _  3.5 

W'A  ~  (H+  10)'7’  '  (3) 

When  an  underwater  charge  is  detonated  at 
moderate  depth,  the  water  surface  first  receives  an 
impulse  from  the  shock  wave.  The  incident  (com¬ 
pressive)  shock  is  reflected  from  the  water/air  in¬ 
terface  as  a  rarefaction  (tensile)  wave,  and  the  sur¬ 
face  of  the  water  “spalls  off,”  creating  a  zone  of 
cavitation  behind  the  reflected  wave  (Fig.  3). 
However,  the  main  surface  disturbance  is  caused 
by  the  ejection  of  water  displaced  by  the  expand- 
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ing  gas  bubble.  If  the  scaled  charge  depth  is  not 
too  great,  the  gas  bubble  causes  an  eruption  of 
water,  spray  and  gases  (Fig.  4). 

The  general  appearance  of  the  overall  water¬ 
spout  changes  with  time,  and  its  character  varies 
with  the  scaled  depth  of  the  charge.  Some  limits 
for  different  types  of  surface  behavior  have  been 
deduced  by  Young  (1971),  as  shown  in  Figure  5.  A 
charge  which  is  sufficiently  far  below  the  surface 
does  not  produce  any  surface  eruption  or  any 
emergence  of  explosion  products.  The  limit  for 
this  kind  of  total  containment  has  been  deduced 
(Fig.  5)  to  be  a  charge  depth  dc  of  at  least  40  /?bm, 
where  /?bm  is  the  maximum  bubble  radius  of  the 
first  pulse,  as  given  by  eq  3.  For  slightly  shallower 
charges,  with  depths  perhaps  in  the  range  40  > 
dc/Rbm  >  25,  there  is  no  surface  eruption,  but  lo¬ 
calized  upwelling  of  explosion  products  occurs.  As 
charge  depth  decreases  into  the  range  25  >  dc/Rbm 
>  7,  the  explosion  forms  a  mound  on  the  water 
surface.  As  dc  decreases  through  this  range,  the 
mound  tends  to  change  from  a  hump  of  turbulent 
water  to  a  squat  dome  of  foamy  water  and  spray. 
Further  decrease  of  dc  changes  the  surface  mound 
to  a  low,  rounded  cloud  of  spray  and  eventually, 
say  about  dc/Rbm  =  4,  this  cloud  starts  to  develop 
plumes  shooting  up  around  its  perimeter.  Within 
the  range  4  >  dQ/Rbm  >  1,  plumes  are  well  devel¬ 
oped,  looking  like  violently  ejected  projectiles  of 
water  which  break  up  to  form  individual  trails  of 
spray.  When  the  charge  is  at  very  shallow  depth, 
with  dc  <  Rbm,  the  bubble  bursts  directly  through 
the  surface  and  forms  a  coherent  vertical  column 
of  spray.  At  the  heart  of  the  column  there  may  be 
a  dense  central  plume,  known  as  the  jet.  Smoke 
may  emerge  from  the  column,  and  if  the  charge  is 


at  very  small  depth  ( dc  <  0.2  /?bm)  there  may  be  a 


visible  “smoke  crown”  (black  smoke  in  the  case 
of  an  oxygen-deficient  reaction). 

The  terminology  for  water  thrown  into  the  air 
by  an  underwater  explosion  is  not  always  clear  or 
consistent.  This  writer’s  understanding  of  U.S. 
terminology  is  as  follows. 

Upwelling.  Water  and  explosion  products  rising 
to  the  surface  at  a  point  and  flowing  radially  out¬ 
ward.  There  is  no  perceptible  disturbance  of  the 
surface  level. 

Mound.  A  compact  hump  on  the  water  surface 
consisting  of  turbulent  water,  foamy  water,  or 
foamy  water  plus  spray.  There  is  no  ejection  of 
discrete  globs  or  jets  of  water. 

Plume.  The  term  “plume”  should  probably  be 
reserved  for  the  visible  trail  created  by  the  high¬ 
speed  ejection  of  a  jet  or  glob  of  water.  This 
plume  is  made  up  of  a  dense  central  core  of  water 


and  a  diffused  boundary  layer  of  spray.  Plumes 
tend  to  shoot  out  in  a  ring  near  the  perimeter  of 
the  waterspout,  either  vertically  upward  or  inclined 
radially  outward  like  the  jets  of  a  fountain.  There 
may  also  be  one  or  more  plumes  in  the  center  of 
the  cluster.  Some  people  use  the  term  “plume”  to 
describe  the  entire  outburst  created  by  an  under¬ 
water  explosion,  irrespective  of  its  form.  This  lat¬ 
ter  usage  undercuts  the  usefulness  of  the  term,  and 
seems  to  ignore  the  basic  meaning  of  the  word 
“plume.” 

Jet.  The  jet,  or  water  jet,  is  the  core  of  the  cen¬ 
tral  plume.  It  can  sometimes  be  seen  emerging 
from  the  top  of  the  general  spra?  cloud,  reaching 
scaled  heights  of  30  m/kg‘/s  (75  ft/lb1/j)  or  more 
when  the  charge  depth  is  fairly  ill. 

Column.  The  waterspout  th  n  up  by  a  very 
shallow  underwater  explosioi  ’ds  to  have  a 
more  or  less  vertical  columna  ^im.  The  term 
“column”  is  applied  to  eithe  entire  water¬ 
spout,  or  to  the  outer  sheath  )  vaterspout, 
recognizing  that  the  column  tyi  ,  nas  a  central 
jet.  The  column  may  have  a  slight  flare  at  the  base 
during  the  earlier  stages  of  its  development. 

Base  surge.  When  the  waterspout  collapses  and 
falls  back  into  the  water,  it  creates  an  annular,  or 
toroidal,  cloud  of  spray  and  drives  this  spray  radi¬ 
antly  outward  as  a  “base  surge.” 

For  charges  lying  on  the  bed  of  shallow  water, 
the  maximum  diameter  of  the  base  of  the  water 
column  Z?max  is  given  by  the  following  empirical 
relation  (Swisdak  1978): 


where  d  is  charge  depth  (m)  and  Dmm,  and  W  are 
in  m  and  kg  respectively.  The  limits  of  applicabili¬ 
ty  for  this  empirical  relation  are  0.08  <  (d/W'A)  < 
0.88.  Empirical  relations  are  also  available  for  the 
overall  height  and  diameter  of  the  waterspout 
from  underwater  TNT  explosions  (Swisdak  1978), 
but  these  are  probably  derived  from  large  explo¬ 
sions  going  to  large  water  depths  (charge  depth  is 
scaled  as  m/kgl/4). 

Studies  of  waterspouts  from  underwater  explo¬ 
sions  (McAnaily  and  Rand  1972,  Outlaw  and 
Strange  1974)  have  given  empirical  relations  for 
the  diameter  of  an  “effective  column”  and  for  the 
internal  jet  of  a  waterspout.  The  effective  column 
is  a  denser  water  column  inside  the  masking  spray 
cloud,  while  the  jet  is  the  high  velocity  vertical 
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Figure  6.  Effect  of  charge  depth  on:  I)  the  maximum 
diameter  at  the  base  of  the  water  column,  Dmax,  2)  the 
‘ 'effective''  column  diameter,  Deff,  3)  the  water  jet 
diameter,  Dj.  (Data  for  Dmax  from  Swisdak  [1978]; 
data  for  Deff  and  Dj  from  Me  Anally  and  Rand  [1972] 
and  Outlaw  and  Strange  [1974].) 


core  of  the  waterspout.  For  deep  water  (and  as  a 
lower  limit  for  shallow  water)  the  effective  column 
diameter  Deff  and  the  jet  diameter  D}  are 


Oeff  /  \  o  1 

—  =  29.6  (— \ 
rc  \rc  ) 

Di  /  rc  \  ° 

~  '  ,4'4  hr) 


There  obviously  ought  to  be  limits  set  on  these  re¬ 
lations.  The  original  data  plots  cover  the  range  2  < 
d/rc  <11. 

Figure  6  displays  the  various  empirical  relations 
for  column  and  jet  diameters. 

In  considering  charge  depth,  it  is  of  some  in¬ 
terest  to  define  a  depth  d„  where  charge  depth  d  is 
equal  to  the  theoretical  maximum  bubble  radius 
for  the  first  pulse,  /?bm.  From  eq  3: 

3  5  WVl 

d.  =  Rbm  =  (7) 

(H+  10) /] 

where  d„  Rbm  and  H  are  in  metres  and  IF  is  in  kg. 
In  Figure  7  this  relation  (with  slightly  different 
constants)  is  plotted  on  logarithmic  scales,  and  it 
is  compared  with  a  simple  cube  root  relation  be¬ 
tween  depth  and  charge  weight.  For  the  range  of 
depths  and  charge  sizes  normally  employed  in  ice 


blasting,  the  simple  cube  root  relation  is  a  good 
approximation.  For  very  big  charges  (over  1  ton) 
at  appropriate  depths,  scaling  with  respect  to  the 
one-fourth  root  of  charge  weight  gives  a  better  ap¬ 
proximation.  .  he  approximation  involved  in  cube 
root  scaling  of  water  depth  effects  is  emphasized 
further  in  Figure  8,  where  Rb  fW'A  is  plotted 
against  absolute  depth.  Also  marked  on  the  graph 
are  charge  weights  for  which  d  =  Rbm. 

The  vertical  velocity  of  the  water  jet  can  be  ob¬ 
tained  from  sequential  photographs.  The  velocity 
of  the  jet  tip  decreases  with  time,  down  to  zero  at 
the  maximum  jet  height.  The  initial  jet  velocity 
(Kj)0  can  be  estimated  by  tracking  back  to  time 
zero;  values  obtained  in  this  way  by  McAnally  and 
Rand  (1972)  and  by  Outlaw  and  Strange  (1974)* 
are  shown  in  Figure  9. 

If  a  horizontal  plate  lies  above  the  ejected  water¬ 
spout,  energy  and  force  are  transmitted  to  the 
plate.  McAnally  and  Rand  (1972)  and  Outlaw  and 
Strange  (1974)  made  observations  on  a  target 
plate,  and  by  suitable  analysis  derived  values  for 
the  force  transmitted  by  the  water  jets  of  1-lb  TNT 
charges  (Fig.  10).  The  numbers  are  quite  consis¬ 
tent,  indicating  that  force  is  a  maximum  at  a 
charge  depth  of  d/rc  *  5,  i.e.  d/W'A  *  0.7 
fb/lb 1/3  *  0.27  m/kg 1/5 .  Scaling  factors  for  the 
water  jet  force  were  also  considered. 

•  Equation  4.3  in  Outlaw  and  Strange  (1974)  seems  incorrect, 
and  Figure  4.3  gives  a  ‘‘shallow  water"  curve  that  is  not  fully 
supported  by  experimental  data. 
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Figure  7.  Depth  d,  which  is  equal  to  the  theoretical  maximum  bubble 
radius  for  the  first  pulse,  Rbm .  The  relation  which  takes  account  of  hy¬ 
drostatic  pressure  is  compared  with  a  simple  cube  root  approximation. 
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Figure  8.  Theoretical  maximum  bubble  radius  Rbm  scaled 
with  respect  to  the  cube  root  of  charge  weight  and  plotted 
against  absolute  charge  depth.  The  arrows  indicate  the  values  of 
Rbm/W  /'  at  which  d  =  Rbm  for  various  charge  weights. 
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Figure  9.  Initial  velocity  of  the  emerging  water  jet  plotted  against  charge  depth  for  un¬ 
derwater  explosions.  < Data  from  McAnallv  and  Hand  11972]  and  Outlaw  and  Strange  1 1974].) 


(UN)  dbM 


35,000 

— i — i — i — i — i — i — i — i — i — i — i — 

150 

— 

1  lb  TNT  Cltorgn 

30.000 

-  — 

O  0 

•  °. 

100 

_ 

0  1  !  l  1 

F 

0  •  .  ' 

P 

20.000 

-  •  0  0  _ 

•  o 

15,000 

_ 

0 

50 

— 

10,000 

5,000 

_ 

°  Deep  Water  (McAnolly  ond  Rood,  1972) 

•  Shallow  Water  (Outlaw  and  Strange,  1974) 

■  i  i  i  i  i  i  i  i  i  i  1 

0 

2  A  6  6 iO  IS 

X 

_ .  i _ i _ i _ i _ i _ 

I _ i _ I _ i _ I _ i _ 1_ 

0  05  10  15 

d/w'/s  (ll/lb'/5) 


Figure  10.  Force  exerted  on  a  horizontal  plate  by  the 
waterspout  from  a  I -lb  charge  of  TNT  fired  at  various 
depths  in  the  water. 
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REGRESSION  ANALYSIS  FOR 
ICE-BLASTING  DATA 

A  basic  requirement  in  ice  blasting  is  to  predict 
the  size  of  the  “crater”  produced  by  a  single 
charge.  The  dependent  variable,  which  is  what  we 
want  to  predict,  is  crater  radius.  The  independent 
variables,  representing  the  input  data  for  the  pre¬ 
diction,  are  charge  size,  explosive  type,  charge 
depth,  and  ice  thickness. 

Because  the  number  of  data  sets  available  for 
regression  analysis  is  rather  small,  it  is  necessary 
to  somehow  reduce  the  number  of  variables. 

Explosive  type  is  the  first  candidate  for  dele¬ 
tion,  since  experience  suggests  that  crater  size  does 
not  vary  greatly  for  different  types  of  common  ex¬ 
plosives.  Even  extreme  variations  of  explosive 
type  may  have  relatively  little  effect,  since  tests 
with  low-pressure  gas-blasting  devices  (Mellor  and 
Kovacs  1972)  indicated  that  these  devices  broke 
just  about  the  same  area  of  ice  as  did  high  explo¬ 
sive  charges  of  comparable  energy  content.  There 
are  equivalence  factors  for  adjusting  the  charge 
weight  of  a  given  explosive  so  that  it  performs  in 
some  respect  equal  to  a  reference  explosive.  How¬ 
ever,  these  factors  differ  depending  on  whether 
the  adjustment  is  for  peak  shock  pressure,  bubble 
energy,  total  energy,  impulse,  time  constant, 
energy  flux  density,  or  whatever.  For  production 
of  the  MM72  curves  and  for  the  present  exercise, 
variation  of  explosive  type  was  ignored.  If  a  sim¬ 
ple  overall  model  for  ice  blasting  should  eventual¬ 


ly  emerge,  the  data  can  be  reanalyzed  with  charge 
weights  adjusted  by  an  appropriate  factor.  If  cube 
root  scaling  is  used,  the  adjustment  may  not  have 
a  large  effect,  since  the  final  factor  is  the  cube  root 
of  a  number  which  is  usually  not  far  from  unity. 

Another  way  to  reduce  the  number  of  variables 
is  to  invoke  some  physics,  so  as  to  express  one 
variable  in  terms  of  another.  In  the  present  con¬ 
text,  an  obvious  expedient  is  cube  root  scaling, 
which  derives  from  similitude  considerations.  One 
possible  objection  to  this  is  that  gravity  effects  are 
not  covered  by  cube  root  scaling,  as  can  be  seen  by 
referring  back  to  the  discussion  of  bubble  size  as  a 
function  of  charge  depth.  However,  for  the  ranges 
of  charge  size  (mostly  <  70  kg)  and  charge  depth 
(mostly  <  10  m)  which  are  represented  in  the  field 
data,  cube  root  scaling  of  charge  depth  seems  to 
be  a  perfectly  acceptable  approximation.  The 
point  is  emphasized  by  Figure  1 1 ,  which  is  a  re¬ 
plotted  version  of  Young’s  (1971)  scheme  for 
characterizing  waterspouts.  Figure  11  gives  the 
proposed  dividing  lines  between  different  forms  of 
surface  eruption  in  terms  of  scaled  charge  depth 
for  a  range  of  absolute  charge  sizes.  Over  the  typi¬ 
cal  range  of  charge  sizes  used  in  icebreaking,  say 
1-100  lb  (0.5-50  kg),  the  value  of  dc/Wv>  for  tran¬ 
sition  from  one  type  of  waterspout  to  another 
does  not  vary  much. 

Perhaps  to  some  people  a  more  troubling  factor 
is  that  the  horizontal  dimensions  for  flexural 
breakage  of  ice  sheets  might  not  necessarily  de¬ 
pend  on  the  magnitude  of  applied  force  or  normal 
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Figure  II.  Re-plot  of  Figure  5  to  show  variation  of  waterspout  charac¬ 
teristics  with  seated  charge  depth. 
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displacement.  In  static  cases,  flexural  dimensions 
are  controlled  largely  by  ice  thickness,  by  the  elas¬ 
tic  properties  of  the  ice,  and  by  the  reaction  of  the 
supporting  water.  Thus,  if  the  underwater  explo¬ 
sion  had  an  effect  similar  to  that  of  a  concentrated 
static  uplift  force,  the  fracture  radius  might  be  in¬ 
dependent  of  the  magnitude  of  that  force.  For  the 
purpose  of  regression  analysis,  cube  root  scaling 
of  linear  dimensions  was  adopted,  but  it  was 
recognized  that  in  some  cases  the  crater  radius 
might  conceivably  include  an  increment  attribut¬ 
able  to  flexure,  and  controlled  by  ice  thickness  and 
ice  properties.  In  particular,  flexural  breaks  by 
charges  of  minimum  effective  weight  or  at  max¬ 
imum  effective  depth  were  of  concern. 

By  ignoring  variations  of  explosive  type  and  by 
scaling  the  crater  radius,  the  charge  depth,  and  the 
ice  thickness  with  respect  to  charge  weight,  we  are 
left  with  only  two  independent  variables: 

Scaled  crater  radius,  Y  =  Rc/ Wl/j 
Scaled  charge  depth,  Xt  =  dc/  w'A 
Scaled  ice  thickness,  X2  -  t/WVl . 

The  chosen  regression  equation  is  a  polynomial 
with  cross-products  and  terms  up  to  the  third 
power: 

Y  =  h  +  b.X'+b^+b,# 

+  b<XlX1  +  biX]  +  btXl 

+  biXlX2  +  b,X,Xl  +  b,X].  (8) 

The  basic  data  used  for  the  analysis  are  given  in 
Appendix  A,  together  with  source  references.  The 
scaled  input  data,  i.e.  Y,  X„  X2,  are  tabulated  in 
Appendix  B,  and  the  weighting  for  replicate  sets  is 
explained.  Appendix  C  gives  the  computer  output 
for  the  initial  regression  analysis. 

The  analysis  assumes  that  X,  and  X2  are  exact 
values  and  that  all  error  is  in  determinations  of  Y. 
The  effective  number  of  data  sets  is  291,  and  the 
initial  number  of  regression  variables  (the  ”b” co¬ 
efficients)  is  10.  With  the  variables  in  English  units 
(see  Appendix  B),  the  best  values  of  the  coeffi¬ 
cients  are  as  follows: 

bo  4.8722  b,  -0.019293 

b ,  0.14566  bo  0.0079326 

b,  0.32645  b-,  0.0021915 

b>  -0.15269  b,  -0.00052360 

b,  -0.0015176  bo  0.00026389 


With  these  coefficients,  the  standard  error  of  Y 
about  the  regression  plane  is  1.268  ft/lb v>.  The 
multiple  correlation  coefficient  r  =  0.7066  (r2  = 
0.4993),  and  the  F-test  value  is  31.136  for  281  de¬ 
grees  of  freedom  with  nine  parameters. 

To  examine  the  relative  significance  of  each 
term  in  the  polynomial,  stepwise  regression  was 
carried  out,  with  coefficients  being  deleted  suc¬ 
cessively  on  the  basis  of  T-test  values  (Appendix 

C) .  Coefficients  dropped  out  in  the  following 
order  (least  significant  first): 

bo,  b\,  bo,  by,  bo,  bo,  b2,  by. 

Dropping  b,  and  b,  resulted  in  a  very  slight  im¬ 
provement  in  the  standard  error  of  estimate  (to 
1.265  ft/Ib'7*)  without  significant  change  in  the 
correlation  coefficient.  For  final  calculation  of  the 
design  curves,  b ,  and  b ,  were  deleted  (Appendix 

D) .  In  the  1972  work,  the  same  two  coefficients 
were  dropped. 

When  the  original  work  was  done  in  1972,  it  was 
obvious  that  the  available  data  were  ill-condi¬ 
tioned  for  regression  analysis,  since  most  tests  had 
been  made  with  conditions  close  to  optimum.  In 
subsequent  field  testing,  the  tendency  to  optimize 
was  maintained.  Experimenters  avoided  the  use  of 
charges  that  are  too  small  for  good  breakage, 
charges  that  are  unnecessarily  large,  and  charge 
depths  that  are  too  great  for  good  breakage.  This 
means  that  the  regression  is  likely  to  be  unreliable 
for  large-scaled  charge  depth,  and  for  extremes  of 
scaled  ice  thickness.  Another  source  of  uncertain¬ 
ty  is  inconsistency  in  defining  and  measuring 
crater  radius.  Some  craters  have  a  dear  demarca¬ 
tion  between  the  highly  fragmented  ice  of  the  cen¬ 
tral  crater  and  the  surrounding  intact  ice.  Other 
craters  have  a  central  zone  with  heavy  fragmenta¬ 
tion,  surrounded  by  ice  which  has  been  flexed  and 
cracked,  but  which  is  still  in  place.  The  transition 
from  obvious  flexural  damage  to  the  zone  of  insig¬ 
nificant  radial  and  circumferential  cracks  may  be 
hard  to  define,  especially  when  the  ice  has  a  thick 
snow  cover.  We  might  add  to  this  list  a  suspicion 
that  observers  occasionally  suffer  from  the  “fish¬ 
erman  syndrome,”  wanting  to  get  the  best  possi¬ 
ble  values  without  actually  cheating. 

Figure  12  shows  how  well,  or  badly,  the  regres¬ 
sion  equation  represents  the  data.  For  each  report¬ 
ed  value  of  scaled  crater  radius  Ya,  the  regression 
equation  calculates  a  value  Yp,  using  the  values  of 
X,  and  X2  which  correspond  to  Ya.  In  Figure  12 
the  predicted  value  Yp  is  plotted  against  Ya.  The 
earlier  tendency  for  tne  1972  regression  equation 
to  underpredict  for  the  largest  craters  (Fonstad  et 


Figure  12.  Values  of  scaled  crater  radius  predicted  by  the  regression  relation  plotted  against  the 
corresponding  observed  values. 


al.  1981)  is  no  longer  apparent,  but  for  very  small 
craters  the  new  regression  equation  still  appears  to 
overpredict.  The  scatter  about  the  1:1  line  in  Fig¬ 
ure  12  is  large,  and  it  is  worth  looking  at  the  data 
points  more  closely  to  see  what  might  cause  the 
greatest  departures  from  the  regression  plane. 

In  Appendixes  C  and  D,  values  of  Ya  and  Yp  are 
tabulated  for  each  data  set.  Alongside  each  pair  of 
Y  values,  the  residual  representing  Ta  -  Yp  is 
given.  Another  value,  in  the  last  column,  gives  the 
residual  divided  by  the  standard  error  of  Y  about 
the  regression  plane.  In  Table  1,  these  residuals 
are  used  to  compare  the  data  from  various  sources 
and  for  various  types  of  explosive.  For  each  dis¬ 
tinct  group  of  data  sets,  the  table  gives  the  relative 
proportions  of  positive  and  negative  residuals,  i.e. 
the  percentage  of  data  sets  for  Yt  >  Yp  and  for  Ya 
<  Yp  respectively.  The  table  also  gives  the  percent¬ 
ages  where  Yt  is  above  or  below  the  standard  error 
of  Yp.  Finally,  the  table  shows  what  proportion  of 
the  data  has  values  of  Ym  above  or  below  a  value 
that  is  outside  a  range  more  than  twice  the  stan¬ 
dard  error  of  Yp. 

Because  some  of  the  residuals  are  very  small, 
the  overall  ratio  of  positive  to  negative  residuals  is 
not  very  illuminating.  A  better  test  is  provided  by 
the  relative  proportions  of  Yt  values  which  fall 


outside  the  positive  and  negative  limits  set  by  the 
standard  error.  Looking  down  the  appropriate 
columns  of  Table  1 ,  it  can  be  seen  that  for  some  of 
the  data  sources  the  Ta  values  agree  well  with  Tp, 
having  most  results  inside  the  limits  of  the  stan¬ 
dard  error,  and  a  balance  between  positive  and 
negative  values  of  the  most  extreme  residuals.  In 
this  category  are  the  results  of  Van  der  Kley  for 
guncotton  (nitrocellulose),  American  TNT,  and 
dynamite;  the  Frankenstein  and  Smith  results  for 
ANFO;  the  results  of  Fonstad  et  al.  for  DM- 12; 
and  Nikolayev’s  results  for  trotyl  (Russian  TNT). 
Balanced  results  with  somewhat  greater  scatter  in¬ 
clude  those  of  Van  der  Kley  for  Dutch  TNT,  and 
of  Mellor  for  TNT.  A  strong  imbalance  with  Ka  < 
Yp  appears  in  the  entries  of  Van  der  Kley  for  gun¬ 
powder,  of  Fonstad  et  al.  for  the  blasting  agents 
Amex  II  (ANFO)  and  Hydromex  (AN/TNT  slur¬ 
ry),  and  of  Wade  for  ANFO.  Since  all  of  these  en¬ 
tries  represent  low  velocity  blasting  agents  rather 
than  true  high  explosives,  they  suggest  that  such 
substances  might  systematically  produce  smaller 
craters  than  high  explosives.  The  only  contradic¬ 
tion  is  the  single  test  result  (or  estimate)  for  ANFO 
by  Frohle,  which  can  probably  be  discounted  in 
this  context.  Looking  for  strong  imbalances  with 
Yt  >  Yp,  and  ignoring  the  single  shot  described  by 
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Frohle,  the  candidates  are  the  results  by  Purple 
(C-4),  Barash  (TNT  and  HBX-3),  and  the  2nd  En¬ 
gineer  Battalion  (TNT).  They  all  represent  high 
velocity  high  explosives,  but  there  are  other  data 
sets  for  TNT,  C-4,  and  DM- 12  which  show  good 
balance  and  good  agreement  with  the  regression 
relations. 

Another  characteristic  which  can  be  examined  is 
the  distribution  of  “wild  points”  which  have  de¬ 
partures  of  twice  the  standard  error  or  more. 
“Wild  points”  for  gunpowder  (Van  der  Kley)  and 
for  AN  blasting  agents  (Fonstad  et  al.)  are  nega¬ 
tive,  and  consistent  with  the  general  distributions 
for  the  data  group.  In  the  opposite  sense  (positive 
residuals),  the  trend  for  “wild  points”  remains 
consistent  with  the  overall  distribution  for  the 
HBX-3  data  of  Barash,  and  for  the  C-4  data  of 
Purple.  In  all  the  other  groups  there  is  no  obvious 
relation  of  “wild  points”  to  explosive  type  or  to 
overall  distributions. 

To  push  the  exploration  a  bit  further,  all  data 
sets  with  a  residual  greater  than  1 .8  times  the  stan¬ 
dard  error  were  identified,  and  the  basic  test  rec¬ 
ords  were  examined.  To  check  for  a  possible 
breakdown  of  cube  root  scaling,  these  “wild 
points”  were  arranged  in  three  groups  according 
to  absolute  charge  weight.  For  charges  s  1  kg, 
there  were  only  two  data  sets  with  positive  resid¬ 
uals,  but  six  sets  with  negative  residuals.  However, 
two  of  the  latter  were  for  low  velocity  blasting 
agents.  For  charges  between  1  and  10  kg,  there 
were  three  positive  and  four  negative  residuals. 
Two  of  the  positive  residuals  were  for  HBX-3  and 
the  other  was  for  C-4.  Two  of  the  negative 
residuals  were  for  gunpowder,  and  two  were  for 
dynamite  when  the  ice  carried  a  very  thick  snow 
cover.  For  charges  heavier  than  10  kg,  there  were 
five  positive  and  five  negative  residuals.  Four  of 
the  negative  residuals  were  for  gunpowder,  and 
one  was  for  C-4.  Two  of  the  positive  residuals 
were  for  ANFO,  two  for  TNT  and  one  for  C-4.  In 
short,  there  is  no  clear  relation  between  the  charge 
weight  and  the  distribution  of  “wild  points.” 

The  same  data  sets  were  examined  for  possible 
dependence  on  absolute  ice  thickness,  scaled  ice 
thickness,  and  scaled  charge  depth,  but  there  were 
no  obvious  relationships. 

One  thing  which  is  fairly  obvious,  both  from 
Figures  13  and  14,  is  the  inability  of  the  regression 
equations  to  predict  reliably  the  limiting  charge 
depth  at  which  crater  radius  goes  to  zero.  There  is 
no  mystery  about  why  this  should  be  so,  since  only 
six  data  sets  give  zero  values  for  Ya  (two  from  Bar¬ 
ash,  four  from  Mellor,  all  for  TNT).  The  scarcity 
of  data  sets  giving  small  values  of  Ya  allows  the 


regression  equation  to  be  controlled  by  more  heav¬ 
ily  populated  parts  of  the  data  domain. 

Not  much  can  be  said  about  the  possibility  of 
errors  in  the  input  data.  The  most  likely  sources  of 
error  are  inconsistencies  in  defining  crater  radius, 
and  imprecise  measurement  of  ice  thickness. 
There  are  indications  that  some  data  sets  might 
have  been  affected  in  this  way,  but  for  present 
purposes  they  have  to  be  accepted  at  face  value. 
Inspection  of  the  data  in  Appendix  A  gives  some 
idea  of  the  possibilities  for  input  error.  Crater  di¬ 
mensions  given  by  Wade  were  measured  to  the 
nearest  10  ft,  giving  an  uncertainty  of  ±8%  to 
±  17%.  The  1977  data  from  the  2nd  Engineer  Bat¬ 
talion  suggest  an  uncertainty  of  ±  8%  in  crater  di¬ 
ameter. 

No  consideration  has  been  given  yet  to  ice  type. 
The  effective  mechanical  properties  of  the  ice  vary 
with  absolute  thickness  of  the  ice  cover,  and  with 
the  temperature  of  the  upper  layers.  With  thick  ice 
and  low  air  temperatures,  the  ice  tends  to  be  elas¬ 
tic  and  brittle.  With  thin  ice  and/or  high  air  tem¬ 
peratures,  the  ice  may  be  relatively  soft,  or  even 
mushy.  Another  oversight  is  failure  to  account  for 
snow  cover  on  the  ice.  Even  though  the  strength  of 
snow  may  be  negligible,  its  mass  could  be  signifi¬ 
cant  where  small  charges  are  involved.  For  exam¬ 
ple,  16  in.  of  snow  at  a  mean  density  of  0.25 
Mg/m’  (specific  gravity  0.25)  is  equivalent  to  4  in. 
of  water,  and  for  a  1-lb  charge  this  might  have  a 
similar  effect  to  increase  of  the  charge  depth  by  an 
increment  of  0.33  ft  (i.e.  0.33  ft/lb l/l). 


GENERAL  FEATURES  OF 
THE  REGRESSION  CURVES 

The  curves  shown  in  Figures  13  and  14,  referred 
to  as  the  MM82  curves,  are  broadly  similar  to  the 
MM72  curves,  but  there  are  some  significant  dif¬ 
ferences. 

Looking  first  at  the  effect  of  charge  depth  (Fig. 
13),  it  can  be  seen  that  for  all  ice  thicknesses  the 
crater  radius  is  greatest  where  charge  depth  is  close 
to  zero;  the  crater  size  decreases  as  charge  depth 
increases.  The  curves  are  considered  unreliable  for 
(dQ/  WVi)  >  7  ft/lb ,/j,  since  there  are  very  few  data 
sets  for  large  charge  depth.  The  curve  shown  for 
(t/wv>)  =  0  is,  of  course,  unsupported  by  actual 
data.  The  curves  for  (t/Wl/!)  =  35  and  40  in./lb'7’ 
have  only  weak  support  from  actual  data.  The  dis¬ 
tinctive  optimum  depth  pattern  of  the  MM72 
curves  has  virtually  disappeared,  and  for  most 
practical  purposes,  zero  charge  depth  is  optimum. 

Turning  to  the  effect  of  ice  thickness  (Fig.  14), 
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there  is  a  very  clear  indication  of  optimum  ice 
thickness  at  all  charge  depths.  There  is  a  slight 
shift  of  optimum  ice  thickness  as  charge  depth 
varies,  but  for  most  practical  purposes  it  seems 
sufficient  to  accept  a  constant  value  of  (t/w'A)  = 
10.5  in./lb'4  as  the  optimum  thickness  when 
charge  depth  is  in  the  range  0  to  5  ft/lb  l/\  The 
MM72  curves  gave  optimum  thickness  of  (t/w'A) 
=  10  in. /lb 1/5  for  small  charge  depths.  For  scaled 
thicknesses  greater  than  11  in. /lb'4,  crater  radius 
decreases  for  all  charge  depths.  The  curves  are 
considered  unreliable  for  thicknesses  greater  than 
30  in. /lb 1/5  because  of  the  shortage  of  data  at  large 
values  of  scaled  thickness. 

For  a  charge  of  optimum  weight  placed  at  opti¬ 
mum  depth,  the  MM82  curves  predict  a  maximum 
crater  radius  of  almost  6.6  ft/lb 1/1 .  The  MM72 
curves  were  a  little  bit  more  optimistic,  predicting 
about  6.9  ft/lb 1/1 . 


USE  OF  THE  REGRESSION  CURVES  AS 
DESIGN  CURVES  FOR  ICE  BLASTING 

The  regression  curves,  Figures  13  and  14,  can  be 
used  to  predict  the  effects  of  under-ice  explosions, 
or  to  select  charge  sizes  and  charge  spacings  for  ice 
demolition.  For  a  start,  we  consider  only  single 
charges. 

In  typical  circumstances,  the  user  can  measure 
or  estimate  the  ice  thickness  /.  If  the  aim  is  to  max¬ 
imize  the  damage  from  a  blast,  the  user  can  plan 
on  firing  the  charge  at  essentially  zero  depth,  i.e. 
directly  beneath  the  ice  cover.  In  order  to  calculate 
the  optimum  charge  weight  IV,  it  can  be  assumed 
that  greatest  efficiency  will  be  achieved  when 
{t/w'A)  =  10.5  in./lb,/j,  so  that 

Wopt  =  (f/10.5)5  (9) 

where  fVopt  is  in  lb  and  t  is  in  inches.  If  the  user 
prefers  to  work  with  SI  units,  the  optimum  condi¬ 
tion  is  given  by 

t/W'A  =  0.347  m/kg 1/5 


or, 


WQVl  =  (t/0.347)1  kg 

when  t  is  in  metres  and  W  is  in  kg.  The  procedure 
can  be  illustrated  by  a  numerical  example. 

Example.  The  measured  ice  thickness  on  Lake 
Jokich  is  1716  in.  What  size  of  explosive 
charge  will  be  most  efficient  for  breaking  the 


ice?  What  will  be  the  diameter  of  the  crater 
made  by  a  charge  of  optimum  weight? 

Ice  thickness  t  =  17. 5  in. 

Optimum  condition  is: 

-4-  =  10.5  in./lb'4 

w7* 


lSA 


i 

10.5 


=  ,b'/, 


10.5 


W  = 


=  4.63  lb. 


If  the  explosive  is  packaged  in  pound  quanti¬ 
ties,  the  probable  practical  choice  would  be  a 
5-lb  charge.  If  the  explosive  is  packaged  in 
kilogram  quantities,  a  2-kg  charge  would  be 
close  enough. 

For  prediction  of  the  crater  radius,  the  ac¬ 
tual  charge  weight  is  used.  The  best  result  is 
likely  to  be  obtained  with  the  charge  in  the 
water  immediately  below  the  ice  cover. 
Reading  from  the  curves  in  Figures  13  and  14 
for  scaled  charge  depths  close  to  zero  and 
t/W>A  close  to  optimum,  the  highest  value  of 
the  predicted  crater  radius  R  is  about  6.57 
ft/lb'4,  i.e. 


~  =  6.57  ft/lb '4. 

W*A 

Taking  W  =  5  lb, 

R  =  6.57  x54  x  =  11.23  ft 

and  the  predicted  crater  diameter  is  22.5  ft  for 
typical  high  explosives. 

In  some  circumstances  the  use  of  optimum  size 
charges  may  be  a  false  economy.  For  example,  if 
the  ice  is  very  thin  the  optimum  charge  weight  will 
be  small  and  the  size  of  a  single  crater  will  be 
small.  This  means  that  in  order  to  blast  a  broad 
area  of  ice,  many  shotholes  have  to  be  drilled  and 
many  individual  charges  have  to  be  placed  and 
connected  to  firing  lines.  Provided  that  plenty  of 
explosive  is  available,  it  may  be  quicker  and 
cheaper  to  use  a  few  big  charges  instead  of  many 
small  charges.  The  general  idea  can  be  illustrated 
by  a  numerical  sample. 

Example.  The  ice  on  the  Slim  Jim  River  is  6  in. 
thick.  Calculate  the  optimum  charge  weight 
for  best  explosive  energy  efficiency,  and 
estimate  the  crater  size  for  a  single  charge  of 
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optimum  weight.  Consider  some  alternative 
charge  sizes,  and  estimate  the  corresponding 
effects. 

Taking  (t/w'A)  =  10.5  in./lb'7’  as  the  op¬ 
timum  condition,  and  noting  that  t  =  6  in.. 


=  0.19  1b. 


This  weight  is  about  one-fifth  of  the  weight  of 
a  standard  1-lb  block  of  military  TNT.  If  the 
cardboard  sheath  of  a  1-lb  TNT  charge  is  re¬ 
moved,  five  small  slabs  of  cast  TNT  are  re¬ 
vealed.  The  optimum  charge  weight  just  calcu¬ 
lated  is  equal  to  only  one  of  these  small  slabs. 

If  one  of  these  optimum  charges  is  placed 
immediately  beneath  the  6-in.  ice  cover,  the 
predicted  value  of  scaled  crater  radius  is  6.57 
ft/lb'7’,  and  so  the  actual  crater  radius  is: 

R  =  6.57  W*A 

=  6.57x0.57 

=  3.75  ft. 

This  means  that  the  crater  diameter  is  7.5  ft. 

A  crater  7.5  ft  in  diameter  is  a  good  return 
for  a  few  ounces  of  explosive,  but  it  is  an  unim¬ 
pressive  hole  in  the  ice  cover  of  a  wide  river.  As 
alternatives  to  the  optimum  charge,  consider 
the  effects  of  a  5-lb  charge  and  a  20-lb  charge, 
each  fired  immediately  below  the  ice. 

Taking  R  =  6.57  tfA  ft/lb'/j,  the  predicted 
crater  radii  for  a  5-lb  charge  and  a  20-lb  charge 
are: 

51b  R  =  6.57x1.71  =  11.23  ft 
201b  R  =  6.57x2.7/  =  17.83  ft. 

The  respective  diameters  are  22  ft  and  36  ft. 

The  approximate  areas  of  ice  broken  by 
charges  of  0.19,  5,  and  20  lb  are  44,  400  and 
1000  ft1  respectively.  The  returns  in  terms  of 
demolished  area  per  unit  charge  weight  are: 
233  ft  Vlb  for  an  optimum  (0. 19-lb)  charge,  79 
ft’/lb  for  a  5-lb  charge,  and  50  ft’/lb  for  a 
20-lb  charge.  However,  the  work  involved  in 
drilling  and  loading  one  shothole  is  about  the 
same  for  any  of  these  charges,  and  it  may  be 
more  efficient  operationally  to  use  a  relatively 
small  number  of  big  charges  rather  than  a  very 
large  number  of  little  charges. 


There  may  be  situations  in  which  the  charge  will 
not  be  at  optimum  depth  for  some  operational 
reason.  In  such  cases,  the  design  curves  might  be 
entered  for  finite  values  of  scaled  charge  depth. 
The  procedure  can  be  illustrated  by  a  numerical 
example. 

Example.  An  ice  bridge  is  to  be  built  across 
Maird  Inlet,  where  the  mean  water  depth  is  4 
m.  The  ice  bridge  will  not  have  continuous  re¬ 
inforcement,  and  its  design  thickness  is  0.8  m. 
Demolition  charges  will  be  laid  on  the  sea  bed 
beneath  the  bridge.  Estimate  the  size  of  a 
single  charge  which  will  be  capable  of  cratering 
the  ice  bridge  efficiently. 

Looking  at  the  curves  in  Figure  14,  the  opti¬ 
mum  value  of  scaled  ice  thickness  (peak  of  the 
curve)  will  probably  be  slightly  greater  than  the 
value  10.5  in. /lb 1/3  used  in  previous  examples, 
but  not  much.  We  can  take  the  value  11 
in. /lb 1/3 ,  but  the  small  difference  from  10.5 
in. /lb 1/3  is  hardly  worth  fussing  over.  Using  the 
design  ice  thickness  of  0.8  m  (31.5  in.)  to  esti¬ 
mate  the  best  charge  weight  from  the  relation 


(t/W>A)  = 

11  in./lb 1/5 

w  = 

(n)  ' 

/31.5N  3 

\Tl) 

— 

23.5  lb. 

For  practical  reasons  we  round  this  value  up  to 
25  lb,  and  then  check  the  probable  effects  of  a 
25-lb  charge.* 

The  charge  will  be  set  approximately  4  m  be¬ 
low  the  water  surface,  and  the  ice  will  extend 
almost  0.8  m  below  the  water  surface.  Thus  the 
depth  of  the  charge  below  the  base  of  the  ice 
(dc)  will  be  about  3.2  m  (10.5  ft).  The  scaled 
charge  depth  (dc/w'A)  will  thus  be 

dc  =  JOS 

Wi/'  25 1/3 

=  3.56  ft /lb 1/1 

and  the  scaled  ice  thickness  for  a  25-lb  charge 
will  be 

I  =  3L5 

IT*  25' 1/1 

=  10.8  in./lb'7’. 


•We  could  equally  well  round  this  off  to  a  10-kg  charge  and 
continue  the  calculation  with  that  value. 
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In  Figure  13,  a  vertical  line  drawn  through  the 
value  3.S6  on  the  horizontal  axis  cuts  the  curve 
for  10  in./lb'/j  at  a  F-value  of  about  5.4 
ft/lb l/j.  It  cuts  the  curve  for  15  in. /lb Vi  at  a  Y- 
value  of  about  5.1  ft/lb >A.  The  F-value  for 
t/W>A  =  10.8  in./lb'/j  is,  by  interpolation, 
about  5.35  ft/lb  /j.  Finally  the  predicted  radius 
of  the  crater  Rc  is  given  by 

~  =  5.35  ft/lb 1/5 

WA 

i.e. 

Rc  =  5.35  W*A 
=  5.35x25 1/1 
=  15.64  ft, 

which  means  a  predicted  crater  diameter  of  31 
ft  (9.5  m)  for  a  single  25-lb  (1 1 .3-kg)  charge  ly¬ 
ing  on  bottom. 

A  similar,  but  less  closely  defined,  problem  is 
given  in  the  following  example. 

Example.  Mines  are  to  be  laid  on  a  river  bed  to 
permit  demolition  of  the  winter  ice  cover  while 
hostile  forces  are  crossing.  Each  mine  contains 
150  lb  of  explosives,  and  the  maximum  record¬ 
ed  winter  ice  thickness  is  21  inches.  Calculate 
the  maximum  water  depth  for  effective  opera¬ 
tion  of  the  mines. 

For  this  problem,  assume  that  there  will  be 
no  enhancement  of  the  explosion  effect  by 
having  the  charge  lying  on  the  river  bed,  and 
make  the  arbitrary  assumption  that  the  mine 
has  to  produce  a  crater  that  is  within  25  Vo  of 
the  maximum  possible  size.  The  scaled  ice 
thickness  is  very  small: 

/  =  _21_ 

fV>A  150 1/4 

=  3.95  in./lb 1/1 . 

A  scaled  crater  radius  25V*  smaller  than  the 
maximum  value  of  6.57  ft/lb 1/5  has  the  value 
4.93  ft/lb  /! .  Referring  to  Figure  14  and  draw¬ 
ing  lines  through  X,  =  3.95  and  Y  =  4.93,  the 
intersection  point  lies  just  below  the  curve  for 
X,  *  3  ft/lb .  Making  an  arbitrary  (and 
slightly  conservative)  choice  of  X,  =3  ft/lb'^ 
for  the  maximum  operational  water  depth,  the 
absolute  depth  value  d  is 

d  -  3  x  (150) 1/1 
-  15.94 
•  16  ft. 


The  mines  should  therefore  be  laid  in  places 
where  the  winter  river  depth  is  not  more  than 
16  ft. 

Another  use  for  the  design  curves  is  prediction 
of  the  cratering  effects  from  an  underwater  charge 
when  all  input  data  are  given,  as  in  the  example 
below. 

Example.  What  size  of  crater  can  be  expected 
if  a  40-lb  charge  is  fired  7  ft  below  the  base  of  a 
19-in. -thick  ice  cover? 

Scaling  the  charge  depth  and  the  ice  thick¬ 
ness: 

dc  =  _7_ 

fy*  40 1/3 

=  2.05  ft/lb‘/j 

t  _  _19_ 

W*  40 1/3 

=  5.56  in./ft 1/3 . 

It  is  convenient  to  use  the  curves  in  Figure 
14,  drawing  a  vertical  line  through  the  X ,  value 
of  5.56  in./lb‘/j.  This  cuts  the  curve  for  X,  =  2 
ft/lb 1/3  at  a  T-value  of  about  5.8  ft/lb 1/3 .  The 
small  difference  between  (d^/  H/J/l)  =  2.05  and 
(dc/WiA)  =  2.0  is  not  worth  bothering  about. 
Thus, 

=  5.8  ft/lb'/3 

tfA 

and 

Rc  =  5.8x40 1/3 
=  19.84  ft. 

The  expected  crater  diameter  is  approxi¬ 
mately  40  ft. 

The  design  curves  are  based  on  test  data  for  ice 
that  is  mostly  no  more  than  1  m  thick,  but  in  the 
absence  of  data  for  thick  ice  some  cautious  extrap¬ 
olations  might  be  in  order,  as  in  the  following  ex¬ 
ample. 

Example.  A  submarine  has  to  surface  through 
ice  which  is  2.3  m  thick,  and  it  is  considered 
necessary  to  pre-break  the  ice  by  releasing  a 
buoyant  charge  and  firing  it  in  contact  with  the 
underside  of  the  ice.  What  size  of  charge 
would  be  appropriate?  What  area  of  ice  would 
be  broken  by  a  single  charge? 

Ice  thickness  is  90.6  in.,  and  the  optimum 
condition  for  a  charge  fued  at  zero  depth  be¬ 
low  the  ice  is  (t/W1'')  =  10.5  in. /lb 1/3 .  Thus 
the  optimum  charge  weight  is 
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=  641  lb 


=  291  kg. 

The  predicted  crater  radius  R  for  an  opti¬ 
mum  charge  at  zero  depth  is  given  by  (R/w‘ /’) 

=  6.57  ft/lb 1/5 : 

R  =  6.57  x  641 1/1 
=  56.6  ft. 

The  charge  would  therefore  break  ice  over  a 
roughly  circular  area  113  ft  (34.5  m)  in  diam¬ 
eter,  and  some  cracking  of  the  ice  could  be  ex¬ 
pected  over  a  wider  area. 

A  safe  standoff  distance  for  the  submarine 
could  be  determined  from  a  graph  of  the  type 
shown  in  Figure  1.  Taking  an  accepted  value 
for  safe  overpressure  on  the  hull,  a  corre¬ 
sponding  value  of  scaled  range  R/Wi/l  can  be 
read  off,  and  this  scaled  value  can  be  con¬ 
verted  to  absolute  distance  by  multiplying  by 
the  cube  root  of  the  charge  weight  (641  lb,  or 
291  kg). 

In  the  regression  analysis,  variation  of  explosive 
type  was  ignored,  although  it  is  fairly  clear  that 
different  explosives  produce  different  results.  The 
design  curves  apply  to  some  undefined  “average” 
explosive.  They  are  likely  to  overpredict  for  pro¬ 
pellants  and  blasting  agents  of  low  velocity  and 
low  specific  energy.  They  are  likely  to  underpre¬ 
dict  for  explosives  with  very  high  specific  energy. 
For  common  explosives  such  as  TNT,  C-4,  good 
ANFO,  sensitized  or  aluminized  slurries,  nitro- 
methane,  PETN,  RDX,  and  nitrocellulose,  the 
specific  energy  of  explosion  (“heat  of  explosion”) 
is  likely  to  be  in  the  range  0.9  to  1 .4  kcal/g  (3.8  to 
5.9  kJ/g),  and  we  can  perhaps  accept  this  as  the 
range  of  specific  energy*  to  which  the  curves  app¬ 
ly.  Within  this  range,  there  is  a  possible  variation 
of  energy  about  the  mean  by  ±22%,  but  because 
linear  dimensions  such  as  crater  radius,  charge 
depth  and  ice  thickness  scale  with  the  cube  root  of 
charge  weight,  they  also  scale  with  respect  to  the 
cube  root  of  specific  energy.  The  variation  of  a 
linear  dimension  with  variation  of  specific  energy 
might,  to  a  first  approximation,  be  given  by  the 
cube  root  of  the  energy  ratio.  In  other  words,  if 
the  reference  value  of  specific  energy  is  1.15 
kcal/g  (4.8  kJ/g)  and  one  extreme  value  is  1.4 


kcal/g  (5.9  kJ/g),  linear  dimensions  would  differ 
by  a  factor  of  (1.4/1.15),/\  or  1.07.  Since  the 
curves  cannot  predict  to  an  accuracy  of  7%,  there 
is  little  point  in  making  adjustments  within  the 
range  outlined  above. 

If  a  propellant  or  explosive  with  very  low  specif¬ 
ic  energy  is  being  used,  it  might  be  worth  making 
an  adjustment.  For  example,  the  heat  of  explosion 
of  black  powder  is  0.6-0.7  kcal/g  (2. 5-2. 9  kJ/g), 
and  ANFO  with  a  very  small  oil  content  might  be 
about  0.6  kcal/g  (2.5  kJ/g)  or  less.  Taking  the 
energy  of  explosion  as  0.6  kcal/g  (2.5  kJ/g) 
against  a  reference  value  of  1.15  kcal/g  (4.8  kJ/g), 
the  adjustment  factor  is  (0.6/1.15)'/j,  or  0.81.  In 
other  words,  the  crater  radius  would  perhaps  be 
about  20%  less  than  the  radius  predicted  by  the 
design  curves. 

Going  to  the  other  extreme,  some  special  highly 
aluminized  formulations  might  have  values  for  the 
energy  of  explosion  up  to  2  kcal/g.  At  this  ex¬ 
treme,  the  adjustment  factor  for  crater  radius 
might  be  about  (2/1.15)'/j  =  1.2,  which  means  a 
crater  radius  20%  bigger  than  that  predicted  by 
the  curves. 


ROW  CHARGES  AND 
PATTERN  CHARGES 

So  far,  the  discussion  has  been  limited  to  the  ef¬ 
fects  of  single  charges,  but  in  dealing  with  real 
problems  it  may  be  necessary  to  use  a  row  of 
charges,  multiple  rows  of  charges,  or  some  other 
kind  of  pattern.  Since  there  have  been  few  system¬ 
atic  tests  with  row  charges  and  patterns,  it  is 
necessary  to  work  largely  from  the  data  for  single 
charges  in  preliminary  design. 

If  a  number  of  charges,  each  of  weight  W,  are 
set  in  a  row  with  spacing  s,  the  results  are  predict¬ 
able  from  Figures  13  and  14  if  s  is  big  enough  for 
the  charges  to  act  independently  (Fig.  15).  If,  by 
contrast,  s  is  relatively  small  and  adjacent  charges 
interact,  it  is  possible  for  the  width  of  the  channel 
formed  by  the  connected  craters  to  be  greater  than 
the  diameter  of  a  crater  formed  by  an  independent 
single  charge.  This  is  sometimes  called  “row 
charge  enhancement,”  although  in  some  cases  the 
row  charge  may  be  less  efficient  than  a  single 
charge.  In  the  limit,  where  s  becomes  very  small, 
the  row  of  charges  acts  like  a  linear  charge,  so  that 
the  blast  effects  spread  cylindrically  rather  than 
spherically,  and  linear  dimensions  scale  with  the 
square  root  of  unit  charge  weight  instead  of  with 
the  cube  root  of  total  charge  weight. 


'See  section  on  specific  energy  (p.  26)  for  further  discussion. 


imsse&Beaxsassi&usi 


Interacting  Charges  Vr^*  l.5,w«3.4Rcl 

*  $•£  j  s  ’* 


Interacting  Charges  S/rci=|.0,w ®3.4RC1 


x  x  >r  x  x  >\  ">x  xr-  x  x  x  x 

1 1  f  i  h  i  i  r 

Wi/vvy_vvv'vVy 


Figure  15.  Geometry  of  row  charges. 
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Figure  16.  Geometry  of  pattern  charges. 


The  value  of  5  for  charge  interaction  has  not 
been  firmly  established,  but  a  reasonable  condi¬ 
tion  for  the  upper  limit  of  interaction  is  (s//?cl)  = 
2,  where  Rc]  is  the  radius  of  the  crater  formed  by  a 
single  independent  charge  of  weight  W.  With 
greater  spacing  the  row  charge  is  likely  to  form  a 
chain  of  craters  rather  than  a  continuous  channel. 
Limited  testing  by  Kurtz  et  al.  (1966)  suggested 
that  the  maximum  channel  width  w  was  achieved 
with  (s//?cl)  =  2,  but  the  broken  ice  was  not  well 
fragmented  at  this  spacing.  With  (s/f?cl)  =  2  the 
channel  width  w  was  about  1.9  £>cl,  where  Dc,  = 
2/?c).  With  (s/Rc j)  =  1.5  and  (s/flcl)  =  1,  the 
channel  width  w  was  about  1.7£>cl.  For  planning 
purposes  it  is  probably  best  to  allow  for  some 
overlap  of  the  craters  that  would  be  formed  by 
single  charges,  i.e.  take  (s/Rcl)  as  less  than  2  but 
not  less  than  1.5.  The  resulting  channel  width 
might  then  be  about  1.7Dcl  =  3.4/?cl. 

When  a  pattern  of  charges  is  laid  out  to  break  a 
wide  area  of  ice,  the  logical  arrangement  is  to  have 
the  charges  at  the  node  points  of  a  square  net,  in 
which  the  spacing  s  between  charges  in  a  row  is 
equal  to  the  spacing  between  adjacent  rows  (Fig. 
16).  If  the  aim  is  to  break  a  wide  strip  of  ice,  the 
charges  may  be  laid  out  in  a  “5-spot”  pattern,  but 
this  is  just  a  square  grid  rotated  through  45  °  rela¬ 
tive  to  the  axis  of  the  strip  (Fig.  16). 


If  the  mesh  size  s  is  very  small,  then  the  pattern 
acts  like  a  sheet  of  explosive,  producing  essentially 
plane  wave  propagation  in  one  dimension  only. 
The  trick  is  to  find  the  value  of  s  which  gives  the 
most  efficient  fragmentation  of  the  complete  area. 

From  what  has  already  been  said  about  the  row 
charges,  a  first  guess  might  be  that  (s/7?cl)  =  2 
represents  a  reasonable  condition.  However,  for  a 
multi-row  pattern  this  is  probably  too  conser¬ 
vative,  because  a  single  row  with  (s/Rcl)  =  2  will 
break  over  a  half-width  of  about  1.9/?cl,  albeit 
with  poor  fragmentation.  The  writer  has  blasted  a 
pattern  with  (s/Rci)  =  2.3  and  judged  the  result  to 
be  an  overkill.  This  suggests  that  the  mesh  size 
should  be  bigger  than  2.3/?cl ,  but  probably  no 
more  than  3.8/?cl .  In  the  absence  of  firm  test  data, 
it  is  suggested  that  a  mesh  size  of  s  =  3/?c,  can  be 
adopted  for  planning  purposes.  This  should  do  the 
job  without  being  unduly  wasteful. 

Example.  Design  a  pattern  charge  capable  of 
breaking  river  ice  up  to  0.5  m  thick,  with  good 
fragmentation.  The  approximate  width  of  the 
river  is  120  m,  and  the  required  width  of  the 
broken  channel  is  45  m. 

The  first  job  is  to  determine  the  required  size 
for  a  single  charge,  and  to  estimate  the  crater¬ 
ing  effect  of  this  charge  when  it  acts  alone.  If 
we  were  concerned  solely  with  ice  0.5  m  (19.7 
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in.)  (hick,  the  optimum  charge,  placed  at  zero 
depth  below  the  ice,  would  be  as  given  by: 

i/WA  =  10.5  in./lb* 

W  =  (r/10.5)’ 

=  (19.7/10.5)’ 

=  6.6  lb 


Whether  or  not  the  lines  of  the  grid  are  parallel 
and  perpendicular  to  the  center  line  of  the  river 
is  immaterial  in  principle,  but  in  practice  the 
grid  lines  should  either  be  aligned  with  the  river 
direction,  or  else  set  at  45°.  For  the  charge 
spacing  (grid  size)  s,  we  take 

s  =  3/?e, 

=  3x3.2 


=  3  kg. 


=  9.6  m. 


Assuming  metric  packaging  and  taking  W  =  3 
kg  (6.62  lb),  the  predicted  crater  radius  in  ice 
0.5  m  thick  would  be,  from  Figure  14: 


WA 


6.57  ft/lb  '/j 


or, 


/?cl  =  6.57  x  6.62  'A 
=  12.3  ft 
=  3.76  m. 

In  this  problem,  (he  charge  design  has  to 
cope  with  any  ice  thickness  up  to  0.5  m,  and  it 
can  be  seen  from  Figure  14  that  the  3-kg  charge 
will  create  a  smaller  radius  in  thin  ice  than  it 
will  in  ice  0.5  m  thick,  which  is  optimum  for 
this  charge  size.  The  thinnest  ice  which  is  of 
concern  is  4  in.,  since  thinner  ice  will  not  safely 
support  men  to  place  the  charges,  and  also  4 
in.  of  ice  or  less  can  be  broken  up  by  bridge 
boats.  With  a  charge  weight  of  6.62  lb,  the 
scaled  thickness  of  4  in.  of  ice  is  2.13  in./lb'/j. 
Checking  the  design  curves,  Figures  13  and  14, 
for  a  charge  depth  of  zero  and  a  scaled  ice 
thickness  of  2.1  in./lb'7’,  the  predicted  scaled 
crater  radius  is  5.6  ft/lb  .  Thus  the  actual 
crater  radius  in  thin  ice  should  be: 

*cl  =  5.6x6.62 1/3 
=  10.5  ft 
=  3.2  m. 

The  charge  weight  has  to  be  sufficient  to 
break  the  thickest  ice  efficiently,  and  3  kg  is  a 
suitable  choice.  However,  the  value  of  *d 
which  will  decide  the  spacing  of  the  charges  has 
to  be  the  lowest  value  expected  for  the  possible 
range  of  ice  conditions,  and  for  a  3-kg  charge 
this  is  3.2  m,  rather  than  the  3.76  m  expected  in 
the  thickest  ice. 

To  obtain  uniform  spacing,  the  charges  are 
laid  out  in  a  square  grid  with  a  mesh  size  s. 


If  the  grid  lines  run  directly  across  the  river, 
the  space  between  cross-stream  rows  is  9.6  m. 
To  break  a  channel  of  width  45  m,  five  rows 
will  be  needed.  In  each  row  there  will  be  about 
12  or  13  charges,  giving  a  total  requirement  for 
60-65  charges,  or  180-195  kg  of  explosive. 

If  the  grid  lines  run  at  45°  to  the  cross¬ 
stream  direction,  then  in  each  of  the  cross¬ 
stream  rows  the  space  between  charges  increas¬ 
es  to  9.6  x  >J2  =  13.58  m.  The  transverse  spac¬ 
ing  between  adjacent  rows  decreases  to  9.6/ V2 
=  6.79  m.  The  direct  distance  between  charges 
remains  exactly  the  same.  To  cover  the  re¬ 
quired  width  of  45  m,  seven  cross-stream  rows 
are  required,  each  of  them  having  about  nine 
charges.  The  total  requirement  is  about  63 
charges,  or  about  189  kg  of  explosive,  which  is 
essentially  the  same  as  with  the  net  laid  out  the 
other  way.  In  cases  where  the  required  channel 
width  is  not  very  much  bigger  than  the  spacing 
between  cross-stream  rows,  one  grid  orienta¬ 
tion  might  be  slightly  more  economical  than 
the  other  because  the  row  spacing  divides  more 
neatly  into  the  channel  width. 


RESPONSE  OF  FLOATING  ICE  SHEETS 
TO  UNDERWATER  EXPLOSIONS 

The  mechanical  effects  of  explosives  on  floating 
ice  have  not  been  studied  in  detail,  and  these  notes 
are  in  the  nature  of  exploratory  speculations. 
Three  things  are  considered:  1)  direct  damage  to 
the  ice  by  the  shock  wave,  2)  gross  displacement  of 
the  ice  by  water  erupting  in  response  to  the  gas 
bubble,  3)  planar  flexure  of  the  ice  by  a  concen¬ 
trated  uplift  force. 

The  initial  shock  wave  from  a  concentrated 
charge  can  be  assumed  to  propagate  through  the 
water  spherically,  attenuating  as  discussed  earlier. 
When  this  wave  reaches  the  ice  it  will  be  reflected 
from  the  lower  surface  of  the  ice,  and  it  will  be  re¬ 
fracted  and  reflected  within  the  ice  layer.  Shock 
measurements  were  made  in  water  just  below  an 
ice  cover  in  a  U.S.  Navy  study  (Barash  1966b). 
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Figure  17.  Shock  wave  paths  for  an  underwater  explosion  with  overly¬ 
ing  ice.  (After  Barash  1966b.) 


The  pressure-time  records  could  be  interpreted  in 
accordance  with  the  scheme  of  ray  paths  shown  in 
Figure  17,  subject  to  the  usual  limitations  of  criti¬ 
cal  angles. 

Two  potential  modes  of  failure  for  the  ice  are: 
1)  internal  cracking  under  impact  from  a  compres¬ 
sion  wave,  and  2)  spalling  of  the  upper  surface  in 
tension  when  the  wave  is  reflected  from  the  ice/air 
interface.  The  amplitude  of  a  plane  wave  which  is 
needed  to  crack  a  wide  sheet  of  ice  is  not  known. 
The  quasi-static  strength  of  ice  under  uniaxial 
compression  reaches  about  10  MPa  at  the  highest 
loading  rates  used  in  testing  laboratories,  but  with 
effective  lateral  confinement  this  value  could  easi¬ 
ly  be  doubled.  From  Figure  1  it  can  be  seen  that, 
for  TNT  in  water,  shock  pressures  in  excess  of  10 
MPa  can  be  expected  at  ranges  up  to  4.3  m/kg  7' 
( *  1 1  ft/lb ,/\  or  about  80  charge  radii).  Pressures 
in  excess  of  20  MPa  can  be  expected  at  ranges  up 
to  2.3  m/kg'/j  (*  6  ft/lb ,/\  or  about  44  charge 
radii).  The  quasi-static  strength  of  ice  under  uni¬ 
axial  tension  is  about  2  MPa  at  the  highest  loading 
rates  which  have  been  studied,  so  that  incident 
shock  waves  which  are  too  weak  to  cause  crushing 
could  conceivably  cause  surface  spalling.  Another 
possibility  is  that  the  ice  might  be  shattered  by 
some  kind  of  interaction  between  compression 
and  tension,  i.e.  by  failure  in  a  stress  field  which 
has  both  tensile  and  compressive  principal  stress¬ 
es.  Studies  of  failure  criteria  for  ice  under  multi- 
axial  stress  suggest  that  failure  in  a  mixed  tension- 
compression  state  can  occur  with  the  tensile  and 
compressive  principal  stresses  both  below  the  uni¬ 
axial  failure  values. 


If  it  is  postulated  that  craters  are  formed  solely 
by  direct  shock-wave  shattering  of  the  ice,  a  sim¬ 
plistic  model  for  charge  depth  effects  can  be  ob¬ 
tained.  Assuming  that  shock  waves  propagate  in 
the  surface  water  layers  as  they  do  through  deep 
water,  the  crater  radius  Rc  can  be  expressed  in 
terms  of  the  maximum  range  for  damaging  shock 
waves 

Rc  =  (P,-tPc)  (11) 

where  dc  is  charge  depth,  and  f.  varies  with  the  ice 
thickness.  dc  and  (,  can  all  be  scaled  with 
respect  to  WA,  but  this  is  unnecessary  if  eq  11  is 
written  as: 

Rc/t,  =  (1  -  (cfc/f«)2l^-  (12) 

This  relation  is  shown  in  Figure  18.  It  predicts  that 
crater  radius  will  be  a  maximum  at  zero  charge 
depth,  which  is  not  unrealistic,  and  that  Rc  will 
decrease  with  increasing  charge  depth.  The 
limiting  value  of  dc,  where  Rc  drops  to  zero,  is 
itself  equal  to  the  maximum  value  of  Rc  (for  dc  = 
0).  The  maximum  scaled  crater  radius  RQ/  W'A  is 


If  R^/W'A  =  6.5  ft/lb l/\  then  iJW'A  =  6.5 
ft/lb A,  and  this  implies  that  the  limiting  shock 
wave  pressure  is  about  18  MPa  (see  Fig.  1).  If 
Rc/W'a  =  UW'A  =  4  ft/lb ,/\  the  limiting  shock 
pressure  is  about  31  MPa.  These  stress  levels  are 
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Figure  18.  Variaiion  of  crater  radius  with 
charge  depth  according  to  the  simple 
model  described  by  eq  12. 


rather  high  in  comparison  with  accepted  values  for 
the  uniaxial  compressive  strength  of  ice. 

While  the  preceding  considerations  do  not  lead 
to  any  very  definite  conclusions,  they  do  suggest 
that  internal  fracture  of  the  ice  by  the  shock  wave 
could  be  a  factor  for  virtually  the  entire  crater 
zone  when  a  high  velocity  explosive  like  TNT  is 
used.  On  the  other  hand,  low  velocity  explosives 
or  propellants  with  low  detonation  or  deflagration 
pressures  are  not  likely  to  produce  much  internal 
shattering  of  the  ice.  Because  the  latter  materials 
can  certainly  break  ice  by  heaving  and  flexing  the 
sheet,  it  is  obviously  necessary  to  consider  pro¬ 
cesses  other  than  shock  wave  shattering. 

Going  to  the  other  extreme  of  rates  and  time 
durations,  we  can  consider  how  the  ice  might 
behave  if  the  explosive  were  to  push  it  slowly  and 
gently  upward,  forming  a  crater  by  flexure.  For 
static  loads  pushing  up  or  down  on  a  floating  ice 
sheet,  the  scale  of  the  flexural  deformation  can  be 
given  in  terms  of  a  “characteristic  length”  X, 
defined  as 

(4 \EI  1/4  _  (  4Et'  \ 1/4 
V*  \  1 2/r(  1  -  v1)) 

=  (fJW'“ 

where  £  =  Young’s  modulus 

I  =  moment  of  inertia  of  the  plate’s 
cross  section 

k  =  foundation  modulus  (unit  weight  of 
water) 

v  =  Poisson’s  ratio 
t  =  ice  thickness. 


For  a  load  applied  to  a  circular  area  of  radius  a, 
the  radius  to  the  extreme  circumferential  crack  Rtc 
is  (Assur  1961) 

Rec  =  2.06X  +  a/3.  (14) 

Gold  (1971)  summarized  field  data  and  deduced 
a  representative  relation  for  X  in  the  form  of  eq  13 
above,  with  substitution  of  E/(\  +  vf  =  7.65 
GPa.  Using  this  value,  together  with  the  substitu¬ 
tion  k  =  10  kN/ms, 

X  =  22.41  tVt  m  (15) 

when  r  is  in  metres.  Assuming  a  compact  loading 
area,  and  thus  ignoring  a/3  in  eq  14, 

£ec*46.3t3/4  m.  (16) 

If  Rtc  is  identified  with  the  crater  radius  Rc  for 
explosive  events,  a  quasi-static  flexural  mechan¬ 
ism  would  give  the  same  size  of  crater  for  any  size 
of  effective  charge,  with  Rc  being  determined  sole¬ 
ly  by  the  ice  thickness.  For  ice  1  m  thick,  Rc  would 
be  46.3  m.  For  thicknesses  of  0.5  m  and  0.1  m,  the 
values  of  Rc  would  be  27.5  m  and  8.2  m  respective¬ 
ly.  These  values  are  much  larger  than  typical  val¬ 
ues  for  real  craters  made  by  explosive  charges  of 
reasonable  size.  In  short,  there  is  neither  qualita¬ 
tive  nor  quantitative  similarity  between  quasi¬ 
static  flexural  breakage  and  explosive  cratering. 

Actually,  common  sense  should  tell  us  that 
static  theory  is  unlikely  to  be  applicable  to  very 
rapid  loadings.  A  cantilever  given  a  swift  karate 
chop  will  not  necessarily  break  at  the  beam  root. 
A  heavy  rock  tossed  from  a  bridge  will  punch  a 
small  hole  through  ice  instead  of  flexing  and 
cracking  a  wide  area.  Thus  we  should  probably  go 
back  to  a  consideration  of  how  the  ice  might  be  af¬ 
fected  by  the  eruption  of  a  waterspout. 

An  easy  case  to  think  about  is  the  situation 
where  the  scaled  thickness  of  the  ice  is  very  small, 
so  that  the  surface  skim  of  ice  has  little  effect  on 
the  venting  of  the  waterspout.  In  this  case,  we 
would  expect  the  crater  diameter  to  be  at  least 
equal  to  the  diameter  of  the  water  column.  In  Fig¬ 
ure  19,  waterspout  dimensions  and  bubble  diam¬ 
eter  are  compared  with  crater  diameters  for  zero 
scaled  ice  thickness  (according  to  the  regression 
equation).  For  small  scaled  charge  depths,  crater 
diameter  is  20%  greater  than  the  bubble  diameter. 
For  very  small  charge  depths,  crater  diameter  is 
50%  bigger  than  the  diameter  of  the  water  col¬ 
umn,  but  the  crater  and  the  water  column  have  the 
same  diameter  at  d/r  *  14.5  (i.e.  d/WVi  *  0.77 
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Figure  20.  Waterspout  characteristics  from  ordinary  underwater  explosions  (see  Figures  5 
and  II)  compared  with  ice  crater  predictions. 


m/kg'7’  *  1.94  ft/lb '/j)-  This  is  actually  quite  close 
agreement,  since  the  relation  for  water  column 
diameter  refers  to  a  finite  distance  above  water 
level,  where  it  may  have  narrowed  somewhat. 

As  a  matter  of  interest,  Figure  20  shows  the 
boundaries  delineated  by  Young  (1971)  for  differ¬ 
ent  types  of  waterspouts  from  underwater  explo¬ 
sions,  and  it  compares  these  boundaries  with  the 
curves  relating  crater  radius  and  charge  depth. 
The  range  of  charge  depths  over  which  ice  can  be 
broken  is  the  range  corresponding  to  columns  and 
plumes  for  ordinary  underwater  explosions. 
Charges  which  would  create  mounds  on  an  open 
water  surface  are  too  deep  to  blow  craters  in  ice. 

If  the  scaled  ice  thickness  is  small  but  finite, 
ejection  of  the  waterspout  will  not  be  significantly 
impeded  by  the  ice,  but  there  will  be  a  tendency 
for  ice  to  “peel  back”  at  the  rim  of  the  crater. 
Thus  it  is  reasonable  to  expect  the  diameter  of  the 
crater  to  be  somewhat  bigger  than  that  of  the 
water  column  for  this  condition. 

When  there  is  an  outer  annulus  of  ice  broken  by 
flexure,  we  must  consider  whether  or  not  the 
width  of  that  annulus  is  likely  to  scale  with  respect 
to  the  charge  weight.  By  analogy  with  static  flex¬ 
ure,  the  width  of  the  flexed  rim  might  depend 
largely  on  the  ice  thickness,  and  not  on  the  charge 
yield.  It  has  already  been  shown  that  flexure  by 
“point  loads”  gives  a  large  radius  to  the  extreme 
circumferential  crack.  For  the  flexure  of  a  long 
floating  beam  or  for  a  semi-infinite  ice  sheet,  the 
distance  from  the  loaded  free  edge  to  the  critical 
crack  is  quite  similar  to  RK  for  the  radially  sym¬ 
metrical  case.  For  a  semi-infinite  sheet  or  a  wide 
beam,  the  critical  length  X  is  as  given  by  eq  13,  and 
the  distance  to  the  critical  crack  is  (7r/2)X  =  1.57X 
instead  of  2.06X  (Mellor,  in  press).  There  is  no 
need  to  repeat  the  arithmetic  that  was  given  for  the 
radially  symmetric  flexure;  the  width  of  annulus 
predicted  by  static  theory  is  far  bigger  than  the  ob¬ 
served  cracked  rims  of  explosion  craters.  It  seems 
more  reasonable  to  expect  that  peripheral  flexure 
will  be  caused  by  “heave”  or  “fallback”  at  the 
base  and  rim  of  the  waterspout,  and  the  dimen¬ 
sions  of  this  water  disturbance  are  likely  to  scale 
with  the  charge  weight. 

Looking  at  the  opposite  extreme  of  ice  thick¬ 
ness,  there  must  obviously  be  some  ice  thickness 
limit  for  crater  formation.  With  a  charge  at  zero 
depth,  i.e.  in  contact  with  the  base  of  the  ice,  a 
very  thick  ice  cover  might  behave  in  a  similar  way 
to  a  semi-infinite  ice  mass  when  a  charge  is  embed¬ 
ded  inside  it.  Test  data  exist  for  explosive  cratering 
of  a  semi-infinite  ice  mass  (Fig.  21),  and  they  in¬ 
dicate  that  crater  radius  drops  to  zero  when  the 


Figure  21.  Variation  of  crater  radius  with 
charge  depth  for  explosions  in  a  semi-infinite 
ice  mass.  ( Data  from  Livingston  1 1960/.) 


charge  depth  reaches  about  5  ft/lb  'A .  By  analogy, 
it  might  therefore  be  expected  that  crater  radius  in 
floating  ice  would  drop  to  zero  when  the  scaled  ice 
thickness  is  around  5  ft/lb 'A  (60  in./lb,/*),  and 
when  charge  depth  is  zero.  The  regression  curve 
for  (Rc/W'a)  versus  (t/WYl)  with  (dQ/Wv')  =  0 
does  not  give  this  result,  but  extrapolation  of  the 
curve  from  its  range  of  validity  would  certainly 
not  conflict  with  the  idea  that  the  limiting  value  of 
(T/w'a)  is  around  60  in. /lb'7’. 

If  the  direct  explosive  attack  (shock  wave  plus 
gas  expansion)  reaches  a  limit  at  scaled  ice  thick¬ 
nesses  somewhere  near  60  in./lb'7\  the  remaining 
question  is  whether  some  kind  of  flexural  break¬ 
age  occurs  at  such  ice  thicknesses.  One  way  of 
looking  at  this  question  is  to  pose  it  specifically  for 
a  1-lb  charge  with  ice  60  in.  thick,  since  there  are 
data  for  the  water  jet  force  imparted  to  overlying 
plates  by  1-lb  underwater  charges  of  TNT  (Fig. 
10).  With  ideal  conditions,  the  force  imparted  to  a 
plate  by  a  1-lb  charge  is  about  25,000  Ibf.  For 
short-term  static  failure  of  an  ice  plate  by  flexure, 
the  failure  force  P  is  in  the  range  (Gold  1971) 

P/t1  =  50  to  I0001bf/in.J  (17) 

with  the  largest  values  probably  most  appropriate 
for  very  rapid  failure.  With  ice  of  60  in.  thickness, 
P  might  therefore  be  in  the  range  180,000  to 
3,600,000  lbf  for  flexural  failure.  This  is  vastly  in 
excess  of  the  water  jet  force  for  1  lb  of  TNT. 

The  only  remaining  question  for  thick  ice  is 
whether  the  explosive  might  crack  ice  by  creating  a 
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wave.  However,  small  amplitude  flexure  is  only 
likely  to  crack  the  ice — it  is  not  likely  to  form  a 
recognizable  crater. 


I 


where  V  is  the  volume  broken  by  charge  weight  H'. 
In  effect,  the  powder  factor  is  the  inverse  of  speci¬ 
fic  energy.  When  applied  to  breakage  of  floating 
ice, 


SPECIFIC  ENERGY  AND 
“POWDER  FACTOR” 


7 rR^t 

Powder  factor  =  -jrr  .  (20) 

rr 


The  concept  of  specific  energy  is  very  useful  for 
assessing  the  efficiency  of  cutting  and  breaking 
processes.  The  specific  energy  £s  is  the  energy  util¬ 
ized  to  cut  or  break  unit  volume  of  material.  For 
an  explosive  charge,  the  energy  can  be  expressed 
as  kW,  where  W  is  the  charge  weight  and  k  is  a 
characteristic  specific  energy  content  for  the  ex¬ 
plosive.  For  breakage  of  floating  ice,  the  volume 
of  material  fragmented  by  an  explosive  charge  can 
be  taken  as  the  crater  area  multiplied  by  the  ice 
thickness.  Thus, 


(18) 


where  Rc  is  crater  radius  and  t  is  ice  thickness. 

In  mining  and  rock  blasting,  a  traditional  term 
is  “powder  factor,”  which  means  the  volume  of 
rock  broken  by  unit  weight  of  explosive,  i.e. 


Powder  factor 


(19) 


If  appropriate  and  consistent  units  are  used  for 
/?c,  t  and  W ,  specific  energy  and  powder  factor 
can  be  written  in  terms  of  the  scaled  variables  as 


k/ 7 r 

(R<./wv,y(t/wv') 


(21) 


Powder  factor  =  ir(Rc/W'Ay(t/W'A).  (22) 


When  applied  to  explosives,  the  term  specific 
energy  is  open  to  some  degree  of  interpretation. 
U.S.  and  Canadian  commercial  handbooks  do  not 
mention  the  term,  and  neither  do  most  of  the  (few) 
textbooks  in  the  field.  Specific  energy  can  be  taken 
as  the  theoretical  work  done  by  expansion  of  the 
gases  produced  by  the  explosion  (see  Meyer  1981). 
Alternatively,  it  can  be  taken  as  the  heat  of  explo¬ 
sion,  which  is  the  total  energy  liberated  in  the  reac¬ 
tion.  The  latter  interpretation  is  used  in  this 
report,  partly  because  values  for  heat  of  explosion 
are  relatively  easy  to  obtain  from  reference  books, 
and  partly  because  it  more  truly  represents  the 
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Figure  22.  Predicted  specific  energy  for  ice  blasting  when 
charges  are  at  optimum  depth.  For  this  plot,  specific 
energy  is  based  on  the  heat  of  explosions. 
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energy  actually  consumed.  The  work  of  gas  expan¬ 
sion  is  roughly  20%  of  the  heat  of  explosion. 

Taking  *  as  the  heat  of  explosion,  values  for  ty¬ 
pical  explosives  range  from  about  0.6  kcal/g  (2.5 
kJ/g)  for  ANFO  with  low  oil  content  and  for  pro¬ 
pellant  powder,  to  almost  2  kcal/g  (8.3  kJ/g)  for 
heavily  aluminized  high  explosive  formulations. 
For  cast  TNT  and  nitromethane,  the  heat  of  ex¬ 
plosion  is  about  1.1  kcal/g  (4.6  kJ/g).  For  RDX 
(Cyclonite)  and  HMX  (Octogen),  k  *  1.3  kcal/g 
(5.4  kJ/g);  for  PETN,  C-3/C-4  (plasticized  RDX) 
and  MS  80-20  slurry,  k  *  1.4  kcal/g  (5.9  kJ/g). 

By  taking  values  of  (RC/W'A)  and  (t/W'A)  from 
the  regression  curves,  and  assuming  a  value  of  k, 
then  £s  can  be  plotted  against  (t/W'A)  for  opti¬ 
mum  charge  depth.  This  plot  (Fig.  22)  gives  the 
best  possible  value  of  £s  according  to  the  regres¬ 
sion.  It  brings  out  the  point  that,  whereas  maxi¬ 
mum  crater  radius  is  achieved  with  t/W'A  *  11 
in./lb'7’,  specific  energy  is  minimized  when  t/WVi 
*  20in./lb'/3. 

For  crater  blasting  in  a  semi-infinite  ice  mass, 
the  best  values  of  £s  are  in  the  range  90-270 
lbf/in.2  (0.6- 1.8  MJ/m!).  A  representative  best 
value  for  optimum  charge  depth  might  be  150 
lbf/in.2  (1  MJ/m1). 

To  put  the  absolute  values  of  specific  energy  in 
context,  we  might  note  first  of  all  that,  with  opti¬ 
mum  conditions,  blasting  is  about  as  efficient  as 
the  best  possible  mechanical  cutting  by  drag-bit 
machines  (e.g.  large  and  well-designed  saws,  mill¬ 
ing  drums,  or  rotary  drills).  £s  for  blasting  is  also 
comparable  to  the  specific  energy  for  icebreaking 
ships  (Mellor  1980).  What  this  means  is  that,  in 
terms  of  energy,  blasting  is  a  very  efficient  way  to 
break  ice.  If  specific  energy  for  explosives  is  based 
on  the  work  of  gas  expansion  instead  of  the  heat 
of  explosion,  specific  energy  for  ice  breaking 
looks  even  more  favorable. 

Just  to  complete  the  picture,  it  is  worth  making 
an  estimate  for  the  specific  energy  of  ice  breakage 
by  slow  flexure.  For  a  static  load  P  acting  on  a 
small  circular  area,  Assur  (1961)  gives  the  vertical 
deflection  at  failure,  wf,  as 


W(  - 


(23) 


where  k  is  the  foundation  reaction  and  X  is  the 
characteristic  length  (see  eq  13).  Thus  the  work 
done  in  flexing  the  ice  to  failure  is 


Work  =  Pw,  =  jgj  . 


(24) 


Assur  (1961)  also  gives  the  radius  to  the  extreme 
circumferential  crack  (see  eq  14).  Ignoring  the 
small  effect  of  the  loading  area,  the  radius  to  the 
extreme  crack,  RK,  is: 

Rec  =  2.06X.  (25) 

If  the  volume  of  ice  broken  Vb  is  taken  as  the 
volume  bounded  by  the  extreme  circumferential 
crack,  then 

Vb  =  rr(2.06X)2/.  (26) 

Finally,  the  specific  energy  for  the  slow  flexural 
breakage  is 

E  =  _ £ _ 

s  8*X27r(2.06X)2/ 


P2 
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106.7*/ 

X4 

P2 

12*(1  _  v2) 

106.7*/ 

£/' 

0.1125 

•  (f)‘ ' 

£/(  1  -  v2) 

Gold  (1971)  summarized  field  data  which  indicat¬ 
ed  that  P/t 2  was  in  the  range  50  to  1000  Ibf/in. 2 
(0.35  to  6.9  MN/m2),  with  the  effective  value  of 
£/(  1  -  v2)  equal  to  1.11x10*  lbf/in. 2  (7.65 
GN/m!).  Substitution  of  these  values  puts  £s  for 
slow  flexure  in  the  range  2.5  x  10‘4  to  0.1  lbf/in.2 
(1.7  J/m'  to  0.7  kJ/m’).  Such  values  for  £s  are 
completely  outside  the  range  for  any  practical 
methods  which  have  been  applied  to  the  cutting 
and  breaking  of  ice.  It  is  interesting  to  note  that 
these  values  are  orders  of  magnitude  lower  than 
the  specific  energy  for  icebreaking  ships,  prompt¬ 
ing  one  to  question  whether  conventional  plate 
flexure  analyses  are  of  direct  relevance  to  ship  ac¬ 
tion. 


SUMMARY  AND  CONCLUSIONS 

For  the  range  of  conditions  and  variables  that 
apply  to  typical  ice  blasting  operations,  test  data 
can  be  scaled  in  a  simple  but  rational  manner  in 
order  to  produce  manageable  design  curves  for  ice 
blasting. 

All  available  test  data  have  been  compiled  and 
analyzed  to  provide  design  curves  which  predict 
crater  size  as  a  function  of  charge  weight,  ice 
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thickness,  and  charge  depth.  The  curves  are  be¬ 
lieved  to  provide  a  reliable  guide  for  optimum  per¬ 
formance  blasting  with  typical  explosives.  Their 
application  to  engineering  problems  is  illustrated 
by  means  of  worked  sample  problems. 

The  design  curves  could  be  improved  by  addi¬ 
tional  data  and  more  refined  analysis.  There  is  a 
need  for  test  data  derived  from  blasts  at  large 
scaled  depths,  and  from  blasts  under  very  thick 
ice.  Systematic  investigation  of  explosives  with 
widely  differing  characteristics  would  be  useful. 

Interaction  of  charges  to  create  channels  or  to 
break  broad  areas  of  ice  has  not  received  much 
study,  but  some  interim  working  guidelines  have 
been  deduced.  These  need  to  be  checked  by  field 
testing. 

The  optimum  specific  energy  for  ice  breaking  by 
explosives  is  surprisingly  favorable,  being  com¬ 
pared  to  the  best  attainable  values  for  ice-cutting 
machines,  and  not  much  different  from  good  val¬ 
ues  for  ice-breaking  ships.  Since  explosives  repre¬ 
sent  energy  which  requires  little  capital  equipment 
for  utilization,  they  may  be  of  interest  for  the  pro¬ 
tection  of  offshore  structures  against  ice  forces. 

The  basic  mechanics  of  ice  blasting  is  still  not 
well  understood.  With  typical  high  explosives  the 
ice  is  likely  to  be  damaged  over  most  of  the  crater 
area  by  the  incident  shock  wave  and  its  reflections. 
However,  the  breakout  and  displacement  of  ice 
appears  to  be  controlled  largely  by  the  eruption  of 
water  impelled  by  the  first  phase  of  gas  bubble  ex¬ 
pansion.  The  characteristics  of  ice  breakage  and 
displacement  can  probably  be  related  to  the  char¬ 
acteristics  of  waterspouts  created  by  ordinary  un¬ 
derwater  explosions.  The  mechanics  of  quasi¬ 
static  flexure  in  floating  ice  sheets  does  not  seem 
to  have  much  relevance  to  the  formation  of  explo¬ 
sive  craters.  In  future  studies  it  might  be  worth  ex¬ 
ploring  the  use  of  theoretical  maximum  bubble 
radius  as  a  scaling  parameter,  taking  water  level  as 
the  datum  for  charge  depth. 

The  basic  data  used  in  the  study  have  been  ap¬ 
pended  to  the  report,  since  some  of  the  original 
sources  are  not  easily  accessible.  Scaled  data  and 
details  of  the  computer  analysis  are  also  given. 
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APPENDIX  A:  BASIC  DATA  ON  ICE  BLASTING 


A1 .  Van  der  Kley 
A2.  Van  der  Kley 
A3.  Van  der  Kley 
A4.  Van  der  Kley 
A5.  Van  der  Kley 
A6 •  Wade 
A7.  Kurtz  et  al. 

A8.  Frankenstein  and  Smith 

A9 .  Purple 

A10.  Barash 

All.  Mellor  and  Kovacs 

A12.  2nd  Engineer  Battalion,  USA  CE 

A13.  Mellor 

A14.  Fonstad  et  al. 

A15.  Nikolayev 


Gunpowder 

Dynamite 

Guncotton 

American  TNT 

Dutch  TNT 

ANFO 

C-4 

ANFO 

C-4 

TNT  and  HBX-3 

Dynamite 

TNT 

TNT 

DM-12,  dynamite,  blasting  agents 
TNT 


Tabl« 

Al  .  Data  of  Van 

der  Kley, 

Gunpowder. 

Table 

A2 .  Data 

of  Van  der  Kley, 

Dynamite. 

Charge 

Ice 

Charge 

Crater 

Charge 

Ice 

Charge 

Crater 

No.  of 

weight 

thickness  depth 

rad  1  us 

No.  of 

weight 

th 1 ckness 

depth 

rad  1  us 

tests 

(kq) 

(m) 

(m) 

(m) 

tests 

(kg) 

<m> 

(m> 

(m? 

5 

7.0 

0.25 

1.25 

2.50 

4 

2.50 

0.15 

0.85 

3.75 

5 

12.0 

0.25 

1 .75 

4.00 

4 

2.50 

0.25 

0.25 

3.50 

5 

12.0 

0.35 

1.65 

3.50 

8 

5.00 

0.25 

0.25 

3.75 

5 

15.0 

0.35 

2.15 

3.75 

12 

5.00 

0.25 

1.75 

4.25 

5 

25.0 

0.35 

2.15 

6.00 

30 

2.50 

0.35 

0.65 

4.00 

5 

12.0 

0.45 

1.55 

3.00 

3 

5.00 

0.35 

1.65 

4.50 

5 

25.0 

0.45 

2.05 

7.00 

16 

2.50 

0.45 

0.55 

4.00 

5 

12.0 

0.55 

1.45 

3.50 

12 

2.50 

0.55 

0.45 

3.25 

5 

25.0 

0.55 

1.95 

7.00 

7 

2.50 

0.65 

0.35 

3.50 

5 

12.0 

0.65 

1.35 

3.50 

5 

2.50 

0.85 

0.15 

2.625 

3 

25.0 

0.65 

1.85 

6.00 

3 

25.0 

1.90 

0.60 

4.25 
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Table  A3.  Data  of  Van  der  Kley,  Guncotton. 


No.  of 

tests 

Charge 

weight 

(kg) 

Ice 

thickness 

(m) 

Charge 

depth 

(m) 

Crater 

radius 

(m) 

6 

0.28 

0.63 

0.32 

1.20 

1 

0.28 

0.52 

0.43 

1.25 

4 

0.28 

0.47 

0.48 

1.425 

6 

0.28 

0.95 

0.00 

1.105 

3 

0.56 

0.63 

0.32 

2.00 

2 

0.56 

0.52 

0.43 

1.90 

6 

0.56 

0.47 

0.48 

1.79 

3 

0.84 

0.63 

0.47 

2.00 

1 

0.84 

0.52 

0.58 

2.50 

2 

0.84 

0.47 

0.63 

2.50 

Table 

A5 .  Data 

of  Van  der 

Kley,  Dutch  TNT. 

Charge 

Ice 

Charge 

Crater 

No.  of 

weight 

thickness 

depth 

radius 

tests 

(kq) 

(m) 

(m) 

(m) 

1 

0.25 

0.40 

1.00 

0.575 

1 

0.25 

0.35 

0.00 

1.25 

1 

0.50 

0.36 

0.00 

1.725 

1 

1.00 

0.39 

0.00 

2.30 

1 

2.00 

0.35 

0.00 

3.075 

1 

0.25 

0.30 

0.00 

1.20 

1 

0.50 

0.30 

0.00 

1.95 

1 

1.00 

0.30 

0.00 

2.65 

1 

8.50* 

0.40 

2.70 

4.75 

1 

0.50 

0.40 

2.00 

0.40 

1 

1.00 

0.40 

2.00 

0.50 

1 

2.00 

0.40 

2.00 

3.00 

1 

3.00 

0.40 

2.00 

3.50 

1 

8.50* 

0.30 

1.50 

5.25 

1 

0.50 

0.40 

1.00 

2.10 

1 

1.00 

0.40 

1.00 

3.15 

1 

2.00 

0.40 

1.00 

3.50 

1 

3.00 

0.40 

1.00 

4.625 

1 

0.50 

0.30 

0.50 

1.95 

1 

1.00 

0.30 

0.50 

2.45 

1 

2.00 

0.30 

0.50 

3.00 

•  M-26  mine. 


Table 

A4 .  Data 

of  Van  der 

Kley, 

American 

TNT. 

Charge 

Ice 

Charge 

Crater 

No.  of 

weight 

thickness 

depth 

radius 

tests 

(kg) 

(m) 

(m) 

(m) 

1 

0.454 

0.09 

0.00 

1.525 

1 

0.454 

0.11 

0.50 

1.70 

1 

0.908 

0.10 

0.00 

2.11 

1 

0.906 

0.14 

0.50 

2.10 

1 

0.908 

0.13 

2.00 

1.25 

Table  A6.  Test  results  reported  by  Wade  (1966), 
ANFO,  a  1 1  charges  80  I b. 

Charge  depth  was  measured  from  snow  surface. 

Depths  have  been  adjusted  so  that  they  are  referred 
to  the  bottom  of  the  Ice  cover. 


Charge 

depth 

(ft) 

Water  depth 

below  Ice 

(ft) 

Approx.  Ice 
th 1 ckness 
(In.) 

Crater  dlam/ 

cracked  dlam 

(ft) 

31 

31 

17 

50/100 

14 

20 

20 

30/100 

6 

18 

17 

60/80 

9 

18 

16 

30/100 

9 

18 

16 

30/100 

9 

18 

16 

30/100 

A  160- lb  charge  of  ANFO  was  fired  on  top  of  the 
Ice,  producing  a  crater  30  ft  In  dlam. 
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Table  A8.  Results  of  tests  on  Mississippi  River  by 
Frankenstein  and  Smith  (1966). 


Explosive:  ANFO 

Note:  No  snow  cover,  melting  Ice,  low  strength 


Charge 

weight 

(lb) 

Water  depth 
below  top  of  Ice 
(ft) 

Charge 

depth 

(ft) 

Ice 

thickness 

(In.) 

50 

14 

11 

19.5 

50 

13 

10.5 

19.5 

33.3 

11 

1.0 

18.0 

Table  A9.  Oata  reported  by  R.A. 
Purple  (1965). 

Explosive:  C-4,  all  charges  5  lb. 
Ice  thickness:  4-6  I n 


Charge 

Crater 

Size  of  Ice 

depth 

dlam. 

fragments 

(ft) 

(ft) 

(ft) 

10 

10 

2 

5 

16 

1.5 

2 

22 

1 

30 

1 

-0.5* 

24 

1.5 

Crater 

dtam 

(ft) 

42.4 

40.1 

32.0 


•On  top  of  Ice,  capped  by  sandbags 


Table  A10.  Ice  blasting  tests  reported  by  R.M.  Barash  (1966),  lake  Ice, 


Charge  depth 
below  top 
of  Ice 
(ft) 


I  ce 

th I ckness 
(In.) 


Water 
depth 
from  top 
of  Ice 
(ft) 


Average 
dlam  of 
open  water 
area 
(ft) 


Average  More  or  less 
dlam  of  than  half  of 
broken  Ice  original  Ice 
area  returned 

(ft) 


Average 
diameter  of 
buck  I  Ing  or 
severe  crack 
(ft) 


*d  -  Ice  thickness  given  In  col.  3,  and  converted  to  feet. 
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Table  All.  Tests  by  M.  Mel  lor  and  A.  Kovacs  (1972) 


Results  of  dynamite  shots  under  lake  Ice  In  New  Hampshire. 

Date  24  February  1971,  mean  snow  cover  approx.  7  In.  All  charges  1  lb  of  40< 
gelatin  dynamite. 


Shot 

no. 

Charge 

depth* 

(ft) 

1  ce 

thickness 

(In.) 

Mean  dlam 

of  hole  or 

cracked  zone 

(ft) 

Remarks 

i 

0 

16 

10.5/12.9 

10.5-ft-dlam  open  hole  with 
depressed  rim  of  12.0  ft  dlam. 

2 

1.5 

19 

14.5 

No  open  hole.  Ice  thoroughly 
broken,  but  fragments  fell  back, 

3 

3 

17 

17.0 

No  open  hole.  Fragments  fell 
back. 

4 

4.5 

19 

4.1/33 

4.0-ft-dlam  open  hole  with  a 
33-ft-dlam  circumferential 
crack.  Flyrock  travel  50  ft  or 
more. 

*  Depth  below  bottom  of  Ice. 


Results  of  dynamite  shots  on  lake  Ice  In  Alaska. 

Date  24  March  1971,  snow  cover  16.5  In.,  Ice  thickness  2.71  ft,  water  depth 
10-12  ft,  charges:  military  dynamite. 

Charge  Charge 

Shot  weight  depth 

no. _ (ft) _ (ft) _ Effect  of  shot  _ 


1  2  2  Circumferential  cracks  to  21  ft  dlam,  slight 

depression  Inside  this  area.  10-ft-dlam  central 
area  domed  and  fragmented. 

2  4  2.5  Circumferential  cracks  to  diameter  of  34.5  ft. 

Slight  depression  28  ft  dlam.  Central  hump  to 
9  ft  dlam.  Open  hole  5  ft  dlam. 


Circumferential  cracks  and  slight  depression 
25  ft  dlam.  Hole  completely  choked  with  Ice 
fragments  10  ft  dlam. _ _ 


Table  A13.  Data  recorded  in  unpublished  CRREl  field  report  by  M,  Mel  lor  dated  24  February  1982 


4  1.81  0.33  0.1  5.0  0.13  13.1  4  14.75  4.50 

2  0.907  0  0  5.0  0.13  13.1  4  17  3.96 

10  18.14  6  1.83  5.75  0.15  13.1  4  51  15.55 

0.6  0.272  0.33  0.10  5.25  0.13  13.1  _ 4  _ 8.5  2.59 


Table  A14.  Data  from  tests  at  Drummond  Lake,  B.C,  (Fonstad  et  a  I,, 
data  given  on  p,  1618-1620  of  vol.  III,  Proceedings  of  POAC  81), 


1981 


Shot 

no. 

Placement 

hole  dlam 

(mm) 

Charge 
we  1 gh  t 
(kq) 

Ice 

thickness 

(m) 

Placement 

depth 

(m> 

Water 

depth 

(m) 

Crater 

radius 

(m) 

i 

152.4 

1.00 

0.439 

0.19 

- 

3.12 

2 

152.4 

2.00 

0.439 

0.15 

- 

4.06 

3 

152.4 

1.25 

0.340 

0.12 

2.28 

2.79 

4 

152.4 

1.25 

0.355 

0.13 

3.45 

2.95 

5 

(1) 

0.50 

0.368 

-0.29 

- 

0.55 

6 

(1) 

1.00 

0.387 

-0.31 

- 

0.74 

7 

(1) 

2.25 

0.343 

-0.27 

- 

1.26 

8 

(1) 

4.50 

0.381 

-0.30 

- 

2.14 

9 

152.4 

0.50 

0.318 

0.00 

4.33 

2.20 

10 

152.4 

1.00 

0.305 

0.00 

4.34 

2.31 

11 

152.4 

2.25 

0.343 

0.00 

4.41 

3.26 

12 

152.4 

4.50 

0.356 

0.00 

4.50 

4.19 

13 

304.8 

0.50 

0.330 

0.00 

4.98 

2.36 

14 

304.8 

1.00 

0.381 

0.00 

4.55 

2.61 

15 

304.8 

2.25 

0.305 

0.00 

4.50 

2.99 

16 

304.8 

4.50 

0.305 

0.00 

4.50 

4.22 

17 

457.2 

0.50 

0.324 

0.00 

4.83 

1.58 

18 

457.2 

1.00 

0.381 

0.00 

4.67 

2.20 

19 

457.2 

2.25 

0.362 

0.00 

4.72 

3.19 

20 

457.2 

4.50 

0.375 

0.00 

4.83 

4.52 

21 

Refrozen 

0.50 

0.381 

0.00 

5.05 

2.10 

22 

(152.4) 

1.00 

0.381 

0.00 

4.88 

2.60 

23 

(152.4) 

2.25 

0.356 

0.00 

4.83 

3.77 

24 

(152.4) 

4.50 

0.381 

0.00 

4.93 

4.53 

25 

152.4 

0.50 

0.381 

0.14 

4.85 

1.78 

26 

152.4 

1.00 

0.318 

0.11 

4.74 

2.31 

27 

152.4 

2.25 

0.343 

0.12 

4.23 

3.44 

28 

152.4 

4.50 

0.330 

0.12 

4.17 

3.84 

29 

304.8 

0.50 

0.318 

0.11 

4.10 

1.92 

30 

304.8 

1.00 

0.305 

0.11 

4.24 

2.58 

31 

304.8 

2.25 

0.406 

c  • 

4.52 

3.37 

32 

304.8 

4.50 

0.381 

0.14 

5.11 

4.27 

33 

457.2 

0.50 

0.368 

0.13 

4.99 

1.91 

34 

457.2 

1.00 

0.381 

0.14 

4.60 

2.40 

35 

457.2 

2.25 

0.381 

0.14 

4.27 

3.26 

36 

457.2 

4.50 

0.387 

0.14 

4.29 

3.96 

37 

Refrozen 

0.50 

0.381 

0.14 

4.14 

1.80 

38 

(152.4) 

1.00 

0.381 

0.14 

4.19 

2.98 

Table  A14  (cont'd) 


Shot 

no, _ 

Placement 

hole  dlam 

(mm) 

Charge 

weight 

(kg) 

Ice 

th 1 ckness 

(m) 

Placement 

depth 

(m) 

Water 

depth 

(m) 

Crater 

rad  I  us 

(m) 

78 

152.4 

4.54 

0.394 

0.82 

4.38 

3.75 

79 

152.4 

2.27  (4) 

0,387 

0.53 

4.34 

3.35 

80 

152.4 

0.57  (5) 

0.356 

0.74 

4.27 

1.09 

81 

152.4 

1.13 

0.387 

0.73 

4.19 

1.71 

82 

152.4 

2.27 

0.356 

0.74 

4.75 

1.12  (6) 

83 

152.4 

4.45 

0.387 

0.81 

4.74 

1.11  (6) 

84 

152.4 

4.45  (5) 

0.343 

0.77 

4.25 

1.15  (6) 

85 

At  Ice 

9.00 

0.395 

-0.39 

- 

1.86 

86 

Surface 

18.00 

0.395 

-0.39 

“ 

2.99 

Un 1  ess 

otherwise  specified,  the 

exp  1  os  1 ve 

used  was  DM-12, 

a  German 

PETN  - 

based  explosive  with  oil  and  grease  plasticizer. 

All  placement  depths  and  depths  of  water  are  measured  from  the  Ice/water 
Interface. 

NOTES: 

1.  In  shallow  (3-In)  hole  at  Ice  surface,  tamped  with  snow. 

2.  Ice  bordering  crater  bent  up,  crater  radius  large,  suspect  measurement 
error, 

3.  Explosive  for  this  series  was  CIL  40 %  Forclte. 

4.  Explosive  for  this  series  wbs  CIL  Amex  II. 

3.  Explosive  for  this  series  was  CIL  Hydromex. 

6.  Crater  radius  small,  suspected  partial  detonation  of  explosive  due  to 

water  seepage  Into  the  charges. 


Data  from  tests  5-8  and  82-86  not  used  for  *N82  regression  analysts 


APPENDIX  B:  SCALED  INPUT  DATA 

A1,:  scaled  charge  depth  (ft/lb v')  X2:  scaled  ice  thickness  (in./lb'7’)  Y:  scaled  crater  radius  (ft/lb />) 


The  data  of  Van  der  Kley  in  Appendix  A  give  mean  results  for  multiple  tests  in  some  cases.  In  this 
table  each  mean  result  from  a  multiple  test  is  entered  \Jn  times,  where  n  is  the  number  of  replications 
of  the  multiple  test. 


Xl 

x2 

Y 

4.ccc: 

lb. 

P  •  3  c  v*  u 

4  .GOOC 

19.230 

2 . 3  0  0  C 

C  •  C  o  0  0  0 

1  b  .  8  6  C 

5 .  c :  o  c 

0.  GGOOC 

13. 750 

5.4800 

0  •  C  0  0  0  ci 

1 1 . 8  C  C 

5 . 8  C  0  0 

0  .  C  0  C  0  c 

8 . 4  1  C  0 

6.1600 

C  •  C  ti  C  0  0 

13.it: 

4 . 8  1  P  C 

0 . 0  C  0  0  c 

10.470 

6.2000 

C  •  G  0  0  0  C 

1  C  .  4  7  C 

6 . 2  0  0  C 

C  .  G  G  G  0  0 

8.500C 

6 . 69  0  0 

3  •  3  3  0  o 

5.53CC 

5.8800 

6 . 3  b  a  c 

15.27G 

1.2700 

5.350G 

12.120 

1.2600 

4.0000 

9.6C00 

6.0000 

3  .  5  0  0  u 

8.4000 

6.1500 

1  .  b  5  0  G 

4.0700 

6.4800 

3 .1800 

15.270 

6.68CG 

2.52CC 

12.120 

7.5200 

2.0C0G 

9.6CCC 

7.G1 CO 

1 .7500 

8.4000 

8.1100 

1  .  59  G  0 

10.470 

6.200C 

1.2600 

8 . 3  0  0  C 

6.1700 

1  .  C  0  0  G 

6.5800 

6.0G0C 

1.2300 

29.100 

4.6200 

1  .2300 

29. 10C 

4.6200 

1.8550 

24. 05C 

4.8100 

1 . 8  4  0  C 

21. 7C  C 

5.4800 

1.8400 

21.700 

5. 48  0C 

0.00000 

43.9CC 

4.25C0 

O.COOCO 

43.900 

4. 2500 

Q.9S1  CG 

23.2CC 

6 . 1 5  0  C 

0\98100 

23.2G3 

6.1300 

1  .’20G 

19.150 

5.82CC 

1 .4700 

17.300 

5.49CC 

1 .4700 

1 7.3GC 

5.49CC 

1  .2500 

20. 15 C 

5.3330 

1  .5450 

16.66C 

6.66CC 

1  .6750 

15.C3C 

6.6t  CO 

0 . 0  0  C  0  0 

3.5400 

5 . 0  G  C  C 

1  .6400 

4.33CC 

5.5800 

0.00000 

3.1250 

5.49CC 

1 . 3  c  0  C 

4.3700 

5.4700 

5.2100 

4  .  C  7  C  l 

3.26C0 

1.6500 

3.9550 

3.3000 

1  .6500 

3.9550 

3.3000 

1.9500 

3.9'  .50 

4.4600 

Xl 

X2 

Y 

1.7300 

6 . 0  2  0  m 

3.3500 

1 .7300 

6 . 0  2  0  C 

3.3500 

1.7700 

4.6600 

6 . 0  6  0  C 

1.7700 

4 . 6  6  0  C 

6  .  C  t C  C 

1.6200 

7 .380  0 

3.9100 

1 . 6  2  0  G 

7 . 3  8  0  C 

3.9100 

1  .680C 

5 . 7  1  0  G 

6 . 0  6  C  C 

1 .6800 

5 . 7 1  0  G 

6. 060  0 

1  •  5  G  0  0 

8.7100 

3.9100 

1.5000 

8 .7100 

3.9100 

1.6000 

6. 7300 

5.1900 

1  .  6  C  0  0 

6.7300 

5.19GC 

0.52000 

19.700 

3.6600 

C.52CCC 

19.700 

3.66CC 

1 .5800 

3. 3500 

6 .9700 

1 .5800 

3  .350  0 

6.97CC 

0.46400 

5. 5800 

6.52C0 

C.4640C 

5.580C 

6.5200 

0.36600 

4 .4000 

5.4900 

0.36600 

4  .4000 

5.49C0 

0.36600 

4.4000 

5.4900 

2.5600 

4.4000 

6.21  CD 

2  .5600 

4.4CCG 

6.2100 

2 .5600 

4.4CGC 

6.2100 

1 .2060 

7.8200 

7. 4300 

1 .2060 

7.82CC 

7.4300 

1 .2060 

7. 8200 

7. 4300 

1 .2060 

7.82CC 

7.4300 

1 .2060 

7.820D 

7. 4300 

2.4200 

6.160G 

6.5800 

2.4200 

6. 1600 

6.5800 

1 .0200 

1  C.OOC 

7.4300 

1 .0200 

1C.  000 

7. 4300 

1 .0200 

10.000 

7.43CC 

1.0200 

10. GOO 

7. 4300 

0.83600 

12.230 

6.C40C 

0.836C0 

12.23C 

6.0400 

0.63600 

12.230 

6.04CC 

0.6520C 

14.500 

6.52CC 

0.65200 

14.500 

6.5200 

0.652C0 

14.500 

6. 52 0  0 

0.27800 

18.900 

4. 8800 

0.27800 

18.900 

4. 8800 

7.2CCC 

3.9500 

5. 8000 

3.2500 

4 .6400 

3.4600 

1.3900 

3.95C  C 

6.9600 

2.C900 

3.7100 

3.4800 

2.0900 

3.710C 

3.4800 

2.0900 

3.7100 

3.4800 

77 

1* 

■  j 

i 

> 

si 

n 

Xl 

*2 

Y 

Xl 

x2 

Y 

i 

1  .  9  5  0  C 

3.S55C 

4  •  4  6  0  C 

3.5100 

4 . 9  7  C  L 

7 . 3 1  C  C 

-* 

1  .8400 

4.690C 

3.9100 

0.00030 

6.6830 

6;  •  £>  h  c  C 

L 

. 

1  •  8  4  C  G 

4  .  6  9  G  3 

3 . 9 1  C  Q 

0.97303 

6.54  :  c 

6.56CG 

j 

{ 

2.2000 

4.3100 

3.84CC 

1.9800 

7 .  o  8  : 

7 . 2  C  C  C 

3 

2.200C 

4.31CG 

3 . 6  4  C  0 

2  .66  30 

6.7801 

6. 78  GO 

t--: 

►  * 

1  •  8  fc  0  0 

3. 6  30  0 

5 .  l  a  o  o 

3 .9500 

6 . 2  5  0  C 

3.7600 

h~ 

1  .6600 

3  .630  3 

5.16,00 

4 .7901 

5.9233 

4 . 6  4  G  C 

s 

f-:- 

1  .4500 

6.4000 

6.92P0 

11.000 

29.000 

0.00000 

: 

6. 9600 

6  •  3  4  0  0 

4.540C 

0.03000 

16. CC  C 

5.8500 

1  .9700 

6.23CC 

5.96CC 

1.5000 

19.000 

7.250C 

•a 

-i 

1 

1  .9300 

b  •  3  2  G  C 

6.75 OC 

3.0000 

17.000 

8.5000 

2.9300 

6.0500 

5.14CC 

4.5000 

19.000 

2.05CC 

» 

3.9400 

6.3500 

3 . 0  4  0  C 

1.5900 

25.600 

3.97C  0 

► 

1  .9500 

6. 1300 

6.9100 

1.5800 

20.500 

2.8400 

3. 8300 

6.000C 

3.6300 

0.  CCOOC 

21.400 

3.29  0C 

■m 

5.7700 

6.9400 

1 .6300 

2.4800 

17.350 

6.71CC 

-4 

I 

1  .9400 

3 .5 2C  C 

8.54 OC 

1.9700 

13. 78  C 

6  •  1  5  0  C 

• 

6.4600 

5.15CC 

4.3700 

1  .4500 

10.160 

7.8600 

»  - 

3 .7700 

5. 4600 

5.0900 

1 .1500 

8.0600 

7. 2000 

* 

2.9900 

5.3000 

5.7600 

8  .COCO 

5.5000 

4.000C 

2.8500 

5.3000 

5. 4500 

4. COCO 

5. GOOD 

6. 65  CO 

Ni 

3.1000 

5.5900 

4.97C0 

2.0000 

5.1C00 

5.55CC 

. 

« 

5.8500 

2.92CC 

2.9200 

0.00000 

5.5000 

5.7500 

r 

2.9200 

2.92C0 

4.68CC 

6.0800 

2.4400 

4.070C 

1.17QG 

2 .9  20  C 

6.4400 

0.0  0  0  00 

2.21C0 

4.91 OC 

0.00000 

2.9200 

8.7800 

1  .8600 

2.0900 

5.3700 

\ 

0.71900 

7.0500 

7.C5C0 

3.7100 

2. 2100 

5.4100 

J 

i 

1  .G700 

8.0600 

7.4000 

5.7800 

2.1200 

4.1700 

1.4370 

6.91CO 

7.0500 

2.9500 

2.210C 

5. 5000 

1.4370 

8.4900 

8.7600 

1.4700 

2.2100 

5.3200 

- 

1.725C 

8. 6300 

7.610C 

0.00000 

2.2100 

5.3BCC 

2.3000 

8.49C0 

7.69C0 

9.4900 

6.520C 

0.00000 

2.675-0 

7.20C0 

7. 4700 

4.7400 

6.5200 

1 .3300 

i 

4.9600 

7.2000 

6.0400 

2.3700 

6.5200 

6.1200 

* 

■ 

0.GC000 

12.500 

6.2500 

0.0  00  00 

7.1100 

4.770C 

■ 

0.13500 

13.00C 

5.8800 

6.0000 

2.5000 

2.56C0 

i  * 

0.7250C 

12.500 

8.130U 

7.5600 

3.7800 

0.00000 

i* 

0.72500 

13.000 

8.1300 

9.5200 

4.7600 

0. 00000 

1.4500 

11.250 

8.5000 

12.000 

5.75C0 

C .COCOC 

• 

ii 

1.4500 

14.500 

10.250 

0.00000 

2.5000 

5. 2500 

"! 

1 

V 

kr 

2.0000 

12.25C 

7.8700 

0.21000 

3 .1500 

4.6500 

2.5000 

11.500 

7.0000 

0.00000 

3.9700 

5.1600 

2.500C 

12.750 

6.8600 

1.7500 

1  .6800 

7.46CC 

3.7500 

14.500 

3.8100 

0.39000 

6.2200 

5.0400 

k 

5.0000 

10.000 

3.6100 

0.48000 

13.270 

7.660C 

5.0000 

10.000 

4.5000 

0.30000 

10.540 

8.1200 

•• 

6.5000 

14.500 

4.5000 

0.28000 

9.5500 

6.5300 

8.0000 

14.000 

4.4400 

0.30000 

9. 9700 

6.9000 

10.000 

13.500 

2.6200 

O.OCOCC 

12. 12C 

6. 99  CO 

k  1 

10.000 

14.500 

5.9400 

0.00000 

9.2200 

5.8200 

r. 

0.00000 

22.500 

6.5000 

0.00000 

7. 9200 

6.2700 

1.4500 

25.000 

8.0000 

0.00000 

6.52C0 

6.4000 

t-* 

k4  * 

1.4500 

30.000 

3.50C0 

0. coooo 

12.580 

7.5000 

- 

2.5000 

22.000 

1.7500 

0.00000 

11.520 

6.5800 

2.5000 

29.000 

3.0000 

0.00000 

7.040C 

5.7500 

'' 

k, 

k: 
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r 

k 

- 

► 

L 

Xi 

X? 

Y 

4 .0000 

30.003 

1.8750 

5. 0030 

23.  OCO 

3 . 5  0  0  C 

5.0000 

23.000 

4.0000 

5.0000 

26. COO 

1 .5000 

5.0000 

30.000 

C. 87500 

6.0000 

3  0. CC  C 

2.2500 

7.5000 

29.000 

0.75000 

10.000 

27.000 

C .OUGOO 

0.0C00G 

6.98CC 

6.9200 

0.44000 

14.520 

5.6500 

0.28000 

9.62C0 

5.8200 

0.23000 

7 .920  0 

6.6200 

0.18000 

6.0500 

5.8600 

0.35C00 

12.  120 

6.1000 

0.28000 

9. 2200 

6.5000 

0.2900C 

9.3700 

6.4800 

0.2100C 

6.98C0 

6.5200 

0.41000 

14.030 

6.0700 

0.35000 

11.520 

6.0500 

0.27000 

8.7900 

6.2700 

0.21000 

7. 0900 

6.0500 

0.44000 

14.52C 

5. 7200 

0.35000 

11.520 

7.5100 

0.27000 

8.7900 

6.7500 

0.20000 

6.5200 

7.1000 

2.2200 

12.840 

5.8100 

1 .7400 

9.9800 

6.6300 

1 .5800 

9.0900 

6.3900 

1.1800 

6.740C 

6.4700 

2. 4800 

15.020 

5.5600 

1 .9700 

11.910 

6.1000 

1 .4800 

8  .9300 

6.6400 

1.1900 

7.2200 

6.3500 

2.8800 

5. 3700 

6.7900 

2.9600 

9.0900 

6.4000 

1.7000 

3.7200 

5.8100 

1.4700 

3.8000 

6.3800 

2.4400 

17.67C 

7.6000 

3.6400 

17.670 

4.4400 

4.8800 

17.090 

2.2400 

7.3200 

17.090 

3.1600 

2.2900 

10.760 

6.4800 

3.0700 

10.760 

5.7000 

3  .8300 

10.760 

3.2300 

Xi 

X2 

Y 

o.oococ 

5.59CC 

6. 4400 

0.  ooooc 

12.350 

5.02CC 

o.oococ 

11.520 

5.5400 

0.00000 

8 .3 5  0  C 

6.1400 

0. ocooo 

6 . 8  7  C  0 

6.900C 

0.00000 

1 4.52 C 

6.67C0 

0.00000 

11.520 

6.5500 

0.00000 

8.22CC 

7.2500 

4.6100 

10.760 

2.5500 

0.31000 

9 .4  ICC 

7. HOC 

0.31000 

9.4100 

5.9500 

0.31000 

9.4100 

7.1300 

0.31000 

9.4100 

7.2000 

0.31000 

9.4100 

6.36CC 

0.31000 

9  .4100 

6.5400 

1.4000 

12.410 

5.440C 

1.7300 

9.940C 

6.41  DC 

1 .4200 

8.2200 

6. 7300 

1.0700 

6.0400 

6.1300 

2.0600 

12.760 

6.3100 

1.7300 

9.9400 

5.7400 

1.3500 

8.2200 

5.5400 

1.1300 

6.8700 

6.7500 

2.4000 

13.90G 

3.1300 

1.9800 

11.240 

6.5800 

1.5700 

9.0700 

6.900C 

1 .2500 

7.2000 

5.7100 

1.0200 

8.9CCC 

6.4200 

2.2500 

12.990 

3.3100 

1.770C 

11.240 

4.1400 

2.0600 

11.040 

5.4200 

2.4100 

11.040 

5. 5000 

2.7500 

11.040 

5.8C0C 

3.0100 

11.040 

5.67CC 

3.4400 

11.04C 

5.5500 

3.8700 

11.040 

4.5200 

2. 2700 

10.710 

5.3100 

2.600G 

10.710 

5.6400 

2.8400 

10.710 

5.4000 

3  .2400 

10.710 

5.0300 

3.6500 

10.710 

4.5800 

2.4600 

11.380 

6.1200 

2.4300 

10.710 

6.1200 
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APPENDIX  C:  INITIAL  REGRESSION  ANALYSIS  USING  COMPLETE  POLYNOMIAL 


n  OF  OBSERVATIONS  =  291 
tt  OF  VARIABLES  =  1C 
CATA  FORMAT  =  <10013.5) 


THE  MATRIX  OF  CORRECTED  SUMS  OF  SQUARES  AND  CROSS  PRODUCTS 


VAR  1 

VAR  2 

V  Ak  3 

VAR 

4  VAR 

5 

VAR  6 

VAR  7 

VAR  6 

VAR 

9  VAR 

10 

VAR 

1 

1356.984 

VAR 

2 

316.4450 

13506.07 

VAR 

3 

11122.82 

3647. 5C9 

1 C  7  5  0  6 . 6 

VAR 

4 

17423.44 

37926.35 

151 712.9 

409D46.1 

VAR 

5 

12729.39 

4  4  C  2  4  5 . 8 

150691.2 

1193606. 

C.1626475E 

08 

VAR 

6 

96397.86 

0.1031309E 

06 

37844.24 

1022243. 

1362633  . 

1454390. 

VAR 

7 

146475.1 

233512.9 

1462337. 

3336613. 

7753955. 

C.1416536E 

06 

0 .3122C43E 

08 

VAR 

8 

336311. C 

11C9542. 

3C80C69. 

0. 1 01 8639E 

08 

C .3667886E 

08 

0.2841263E 

06 

0.840C481E 

08 

0.2795254 E 

0^ 

VAR 

9 

313493.7 

0 .1375857E 

08 

4340306. 

0 .3238C93E 

08 

C.5622751E 

09 

0.4307228E 

08 

0.22C2168E 

09 

0.1042 44 0E 

1  0 

0.21 0  982  2E 

11 

VAR 

1C 

-b77 .0793 

-621 .6602 

-5630.070 

-10561.18 

-35235.82 

-50474.70 

-85460.37 

-243978.2 

-1176479. 

902. 6241 

THE 

INVERSE 

OF  THE  ABOVE 

MATRIX 

THE 

GAUSSIAN 

MULTIPLIERS 

-  CIJ 

VAR 

1 

VAR 

2 

VAR 

3 

VAR  4 

VAR 

5 

VAR 

6 

VAR 

7 

VAR 

8 

VAR  9 

VAR 

10 

VAR  1  0.3164392E-01 

VAR  2  0.4200333E-Q2  0 .6523490E-02 

VAR  3-0.519C748E-02-0.1691900E-03  0 . 1 3 9 6 3 1 1 E -0 2 

VAR  4-C.1791126E-02-0.6130154F-03  C  .75 29 1 0 7E -04  0 . 2436C4 5E- C 2 

VAR  5-0.1291725E-03-0 .3769453E-03  C .3231 437E-05  0 . 2 262274E- 04  C .2406608E-04 

VAR  6  C .2597093E-03  0 . 1 1 C 84 69E - 04 - C .8 5 86747E -04  0 . 1 1 4 5565E - 05- 0 .8871 812E-06 
C. 62209751-05 

VAR  7  C.6171819E-04-C.1018447E-04-0.9329943E-05-0.7C26746E-C5  0.1565851E-05 
-0.2629901E-06  G  .  1 174 5 14E - 0 5 

VAR  8  0 .288771 3E-Q4  0 . 23 1  1 4 0 6E - 0 4  0  .2 7 84 674 E -06- 0 . 52 0 2772E - C 5- 0 . 1 30 4 508E -05 
0. 372 360 8 E- 07-0.1 253 193E-C 6  C .20 03564E-06 
VAR  9  0. 1222075 E- 0b  0 .56 C62 1 2E-C 5  -  0 . 2 3 69 765E -0 7- 0 . 237 64 1 4 E- C6- C . 3 76 72 0 1 E- 06 
0.1568884E-07-0.2940027E-07  0 . 1 8 58 839E -0 7  0 . 6 1 C 649  IE  -  08 
VAR  1C-0.3223029E-03-0  .7223418E-03  0 . 33 78 536 E -0 3  0 .3357886E-C5  0 .4 268880E-04 
-C.1755235E-04-C.4849C75E-05  0 . 1 1 585 73E-0 5- 0 . 5838990 E - 0 6  C.2212690E-02 


THE  SAMPLE  PARTIAL  CORRELATION  COCFF ICE  NT S J R I J 
RIJ  =  CIJ  /  (  CII  *  CJJ  ) *  * • 5 

THESE  REPRESENT  CORRELATIONS  BETWEEN  THE  VARIABLES 
THE  CLOSER  TO  ♦-  1.  THE  MORE  DEPENDENCE 
THE  CLOSER  TO  0  THE  MORE  INDEPENDENCE 

VAR  1  VAR  2  VAR  3  VAR  4  VAR  5 

VAR  6  VAR  7  VAR  8  VAR  9  VAR  10 

VAR  1  1.000000 

VAR  2  0.2923469  1.GC0C00 

VAR  3-0.7808969  -0 .56C5878E-01  l.OOGCOO 

VAR  4-C. 6451168  -0.4862826  0.129C951  1.G0  000  C 

VAR  5-0.1480206  -0.9513397  0 .1 762 796E -0 1  C. 2954609  1.000000 

VAR  6  0.5853466  0 . 55 0 242 7E - 0 1  -C .92 1 3 1 55  0 .294271 3E-C1-C .7250709E-01 

1.CCG000 

VAR  7  0.4238814  -0.1163507  -0.2303874  -0.4154153  0.2945227 

-0.9729284 E- 01  1.000000 

VAR  8  0.3626664  0.6393440  0 . 1 6 64 8 75E -0 1-0 . 74 4 7158  -0.5940766 

G.  33 35286E- 01-0.258 3374  1.0C0  0CC 

VAR  9  0  .8791379E-01  0.8882467  -0  .8 1 1555 7 E -02-0 . 1948425  -0.9826984 

0 .8049443E-01-0 .3471573  0.5314280  1.000000 

VAR  10-0. 3851756E-01-0. 1901264  0.1922106  0 . 4573651 E- 02  C. 1849912 

-0.1496048  -0.9511947E-01  0 . 55 0 25 1 9 E  -  0  1  -  0 . 1 58 8479  1.000000 


THE  STANDARD  ERRORS  OF  THE  ABOVE  RIJ  .*  SIJ 

SIJ  =  <  (  1  -  R I J  *R I J  )  /  DOF  RESIDUAL  I  *  * • 5 

VAR  1  VAR  2  VAR  3  VAR  4  VAR  5 

VAR  6  VAR  7  VAR  8  VAR  9  VAR  10 

VAR  1  0.0000000 

VAR  2  0 .570  488 1 E- 0 1  O.OOOOCOO 

VAR  3  0.3726409E-01  0 . 5956 1 1 9E-Q  1  0.00  00000 

VAR  4  0 .4558143E-01  0 .5212663E-01  0  .59 15582E -0 1  O.OOOOCOO 

VAR  5  0.5899786E-C1  0 . 1838234E-0 1  0  .59 6  4  573 E -0 1  0 . 56991 69E- 0 1  C.OCOOOOO 

VAR  6  0.4836727E-01  0  .5956462E-0 1  0 .23 19481E-01  0 . 596291 6E- 0 1  0.5949798E-01 
0.0000000 

VAR  7  0.5403059E-01  0 .5924983E-0 1  0 .58 05022E-01  0 . 54264 09E - 0 1  C .570 0897E-01 
0 .5937 198E- 01  0.0000000 

VAR  8  0.5559363E-01  0 .4586997E-01  0 .59 6 4673 E-0 1  0 . 398 1 266E- 0 1  0 . 4 798699E-0 1 
0.5962181E-01  0 .57630 0 OE- 0 1  0.0000000 

VAR’  9  0  «  5942402E-01  0 . 2 74  0352E-0 1  C . 596 5303E -0 1  0 . 585 1 1 68E -0 1  0 . 1 1 04887E-01 
0.5946142E-01  0 .5594487E-01  0 .5053396E-01  0.0000000 
VAR  10  0.5961073E-01  0 . 5856687E-0 1  0 .58 54265E-0 1  0 . 596543 7E -0 1  0 .5862536E-0 1 
0.5898363E-01  0 .5938451E-01  C .5956462E-01  0 .5889756E-0 1  0.0000000 


[* 

L 

1  ACTUAL 

PRE  CICTE  R 

RESIDUAL  ST 

.  RESIDUAL 

2.300C00 

4.322512 

-2.022512 

-1.594794 

K 

^  • 

5.000000 

6 .15677C 

-1.156770 

-0.9121377 

5.4800CC 

6.399461 

-0.9194613 

-0.7250149 

5  .799999 

6 .471604 

-0.6716051 

-0.5295751 

6.160000 

6.410108 

-0.2501 C78 

-0.1972153 

Jpi 

4.809999 

6.427693 

-1.617694 

-1.275586 

6.200000 

6.478167 

-0.2781668 

-0.21934C5 

« 

6.2000CC 

6.478167 

-0.2781 668 

-C. 2  193405 

;■  7 

6.690000 

6 .403580 

0.2864  199 

0.2258482 

.  ’ 

5.879999 

5.322294 

0.5577C49 

0.4397622 

m 

1.270000 

3.518832 

-2.248832 

-1.773252 

m 

1.2b0CGC 

4.524147 

-3.264147 

-2.573850 

■  ; 

6.000000 

5.194230 

0.8057699 

0.6353668 

■ 

6.129999 

5 .44  C  561 

0.6894178 

0.5436207 

* 

6.480000 

5.699246 

0.7807531 

0.6156405 

6.679999 

5.349C03 

1.330997 

1  .  C  4  9  5 1 9 

m 

7.919999 

5.917610 

2.002389 

1.578926 

W 

7.009999 

6 .1 6414  1 

0.8458586 

0.6669775 

8.109999 

6.208746 

1.901253 

1.499178 

6.2C00G0 

6.297116 

-0.97118380 

01-0.7657993E 

6.169999 

6.328126 

-0.1581268 

-0.1246864 

L* 

6.000000 

6.242813 

-0.2428131 

-0.1914634 

4.620000 

3.997135 

0.6228652 

C. 4911425 

4.620000 

3.997135 

0.6228652 

0.4911425 

4.6.09999 

4  .676984 

0.1330156 

0.1048857 

5.480000 

5.015120 

0.4648800 

0.3665678 

■  v 

5.480000 

5.015120 

0 .4648RC0 

0.3665678 

4.250000 

4.348463 

-0.98463C6E- 

C 1 -0 . 77640  24  F 

i 

4.250000 

4  .348463 

-0.9846306E- 

01-0. 77640 24E 

6.129999 

5.098335 

1.031664 

0.8134890 

6.129999 

5.098335 

1.031b64 

0.8134890 

•’  ’ 

5.820000 

5.627734 

C. 1922655 

0.1516055 

5.49000C 

5.834403 

-0.3444033 

-0.2715694 

5.490000 

5.834403 

-0.3444033 

-0.2715694 

p 

5.330000 

5.499974 

-0.1699743 

-0.1340284 

6.660000 

5.889662 

0.7703381 

0.6074280 

6.660000 

6.025513 

0.6344872 

0.5003066 

■  ‘ 

5.000000 

5.797777 

-0.7977772 

-0.6290644 

5.580000 

5.807294 

-0 .22  72940 

-0.1792262 

5. 450  00  0 

5.712016 

-0.2220163 

-0.1750647 

L  *• 

4 

5.469999 

5.895696 

-0.4256964 

-0.3356707 

3.260000 

3.799948 

-0.5399480 

-0.4257606 

3.300000 

5.738327 

-2.438327 

-1.922673 

3.300000 

5.738327 

-2.438  327 

-1.922673 

4.459999 

5.645411 

-1.185412 

-0.9347231 

4.459999 

5.645411 

-1.185412 

-0.9347231 

4 

3.910000 

5.807120 

-1.897120 

-1.495920 

-  * 

3.910000 

5.807120 

-1.897120 

-1.495920 

3.840000 

5.617790 

-1.777791 

-1.401826 

3.840000 

5.617790 

-1.777791 

-1.401826 

5.179999 

5.613799 

-0.4337997 

-0.3420604 

5.179999 

5.613799 

-0.4337997 

-0.3420604 

4 

3.350000 

6.022786 

-2.672786 

-2.107549 

r 

3.350000 

6.022786 

-2.672786 

-2.107549 

6.059999 

5.824012 

0.2359877 

0.1860813 

6.059999 

5.824012 

0.2359877 

0.1860613 

p 
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3*910000 

6.183818 

-2.273816 

-1.792954 

3.9100CC 

6.183818 

-2.273818 

-1.792954 

6.059999 

5.999811 

0. 60188290- 

Cl  0.474 59 75E 

6.  C59S99 

5.999611 

0.6018829E- 

Cl  C.4745975E 

3.910000 

6.291606 

-2.361606 

-1.877948 

3.910000 

6.291606 

-2.381606 

-1.877948 

5.1900C0 

6.133983 

-0.9439831 

-0.7443509 

5.19  0  OCC 

6.133582 

-0.9439831 

-0. 7443509 

3.660000 

5.759581 

-2.099581 

-1.655565 

3.660000 

5.759581 

-2.C99581 

-1.655565 

6.969999 

5.640507 

1 .329493 

1.048333 

6.969999 

5.640507 

1 .329493 

1.C48333 

6 • 52  0  0  C  0 

6.165604 

0.3543959 

0.2794487 

6.520C0C 

6.165604 

0.3543959 

0.2794487 

5.49C0CG 

5.985950 

-0.4959497 

-0.391C670 

5.490000 

5 .98555  C 

-0.4959497 

-C.391C670 

5.490000 

5.985950 

-0.4959497 

-0.3910670 

6.2C9999 

5  .483028 

0.7269707 

0.5732318 

6.209999 

5 .483028 

0.7269707 

0.5732318 

6.209999 

5.483C28 

0.7269707 

C. 5732318 

7.429999 

6.310980 

1.119C20 

0.8823707 

7.429999 

6.310980 

1.119020 

0.8823707 

7.429999 

6.310Ci80 

1.119020 

0.8823707 

7.429999 

6.310580 

1.119020 

C. 8823707 

7.429999 

6.310980 

1.119020 

0.8823707 

6 • 58  0  C  0  0 

5.791819 

0.7881813 

0.6214978 

6.580000 

5.791815 

0.7881813 

0.6214978 

7.429999 

6.423399 

1.006600 

0.7937258 

7.429999 

6.423399 

1.006600 

0.7937258 

7.429999 

6.423399 

1.006600 

0.7937258 

7.429999 

6.423399 

1.006600 

0.7937258 

6.040000 

6.419282 

-0.3792620 

-0.299C72C 

6.040000 

6.419282 

-0.3792820 

-0.2990720 

6.040000 

6.419282 

-0.3792820 

-0.2990720 

6.520000 

6.313617 

0.2063828 

0.1627372 

6.520000 

6.313617 

0.2063828 

0.1627372 

6.520000 

6.313617 

0.2063828 

0.1627372 

4.879999 

5 .904298 

-1.024299 

-0.8076814 

4.879999 

5.904298 

-1.024299 

-0 . 80  768 14 

5.799999 

2.317914 

3.482C85 

2.745698 

3.480000 

5.178669 

-1.698669 

-1.339437 

6.959999 

5.802749 

1.157250 

0.9125166 

3.480000 

5.549845 

-2.069845 

-1.632117 

3. 480  00  0 

5.549845 

-2.069845 

-1.632117 

3.480000 

5.549845 

-2.069845 

-1 .632117 

7.309999 

5.086012 

2.223988 

1.753662 

6.580000 

6.270692 

0.3093081 

0.2438960 

6.559999 

6.242435 

0.3175640 

0.2504060 

7.200000 

6.049043 

1.150957 

0.9075543 

6.780000 

5.647773 

1.132227 

0.8927851 

3.760000 

5.000769 

-1.240769 

-0.9783729 

4.639999 

4.416367 

0.2236328 

0.1763392 

6.919999 

6.139118 

0.7808809 

0.6157413 

4.540000 

2.984984 

1.555016 

1.226163 

5.959999 

5.970359 

-0.1035976E- 

01-0.8 168897E 

6.750000 

5.993528 

0.7564716 

0.5964940 

5.139999 

5.545798 

-0.4057989 

-0.3199812 
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3.C40003 

5.018438 

-1.978438 

-1.560041 

6.910000 

5.966016 

0.9439840 

0.7443515 

3.600000 

5.042548 

-1.442949 

-1.137795 

1.630000 

3.877679 

-2.247679 

-1.772343 

8.539995 

5.565712 

2.974287 

2.345269 

<.370000 

3.116086 

1.253513 

0.9887375 

5.089999 

5.0C6434 

C .83564760- 

Cl  0.6589259F 

5.759999 

5.418088 

0.3419113 

C. 2696044 

5.450000 

5  .487806 

-0.3780F 51E- 

Cl-C. 29811 24 F 

4.969999 

5  .4  01546 

-0.4315462 

-0.3402835 

2.920000 

3.045366 

-0.1293664 

-C.1C20082 

4.679959 

5.017040 

-0.3370409 

-0.2657640 

b. 440000 

5.639977 

0.8C0C  231 

0 .6306353 

8.779999 

5.667517 

3.112482 

2.454258 

7.049599 

6.317602 

C . 7323961 

C.57751C0 

7.400000 

6.349758 

1.050241 

0.8281378 

7.049999 

6.191698 

0 . 8  5  8  3  0 1  2 

0.6767888 

8.759998 

6.297883 

2.462115 

1.941431 

7.61C0C0 

6.226301 

1 .383658 

1 . 091C76 

7.6900C0 

6.020459 

1.669540 

1.316468 

7.469999 

5.687214 

1.782785 

1.405764 

6.040000 

4.454231 

1 .5R5669 

1.250334 

6.250000 

6.453784 

-0 .2037040 

-0.1606880 

5.879999 

6.438199 

-0.5581599 

-0.4401526 

8.129999 

b  .423483 

1.706516 

1.345624 

8.129999 

6 .4  00278 

1.729721 

1.363922 

8.5COOCO 

6 .324 19  C 

2.175010 

1.715672 

10.25000 

6.143608 

4.106392 

3.2375  78 

7.870000 

6.118374 

1.751626 

1.381194 

7,000000 

5 .950857 

1.C49143 

0.R272715 

6.679999 

5.892477 

0.9875221 

0.7786B22 

3.010000 

5.122766 

-1.312767 

-1.035145 

3.810000 

4.584194 

-0.7741942 

-0.6104687 

4.5000G0 

4.584194 

-C.8*194l8t- 

01-0 .6638892F 

4.50GC00 

3.512161 

0.9878387 

0.7789320 

4.44G0CQ 

2.812549 

1.627451 

1.283280 

2. 620 00  0 

2  .332145 

0.2678547 

0.2269796 

5.9400CC 

2 .33C984 

3.609015 

2.845786 

6.5GC0C0 

5.4564C9 

1.043591 

0.8228934 

8.000000 

4  .598853 

3.401147 

2.681877 

3.500000 

3.730584 

-0 .2,->05841 

-0.1818205 

1 • 75  0  0  C  0 

4.644468 

-2.89446 8 

-2.282351 

3. 0C0  00  0 

3.270216 

-0.2702165 

-0.2130715 

1.875000 

2.059444 

-0.1844444 

-0.1454384 

3.500000 

2  .988802 

0.5111976 

0.4030902 

4.000000 

2.988802 

1.011198 

0.7973509 

1.500000 

2.316247 

-0.616247C 

-0.6436262 

0.8750000 

1.389722 

-0.5147221 

-0.4058694 

2.250000 

0.7840884 

1.465912 

1.155903 

0.75000C0 

0.3426785 

0.4073215 

0.3211817 

o.oocooco 

0.6266217 

-0.6266217 

-0.4541C47 

0.0000000 

0.5972332 

-0.5972332 

-0.4709312 

5.849999 

6.237408 

-0.3874083 

-0.3054797 

7.250000 

5.588769 

1.661231 

1.309916 

8.500000 

5.223693 

3.276307 

2.583438 

2.050000 

4.069355 

-2.019355 

-1.592305 

3.970000 

4.388676 

-0.4186759 

-0.3301349 
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2.84C0CG 

3.290000 

6.709999 

6.15000C 

7.860000 

7.200000 

4.000000 

6.650000 
5.549999 
5.7500OC 
4.07000C 
4.9100Q0 
b. 370000 
5.410000 
4.169999 

5.500000 
5.320000 
5.379999 

0. 0000000 
1.330000 
6.120000 
4.770000 
2.560000 
0.0000000 
0.0000000 
o.ocacooo 

5.250000 

4.650000 
5.160000 
7.459999 
5.040000 

7.860000 
8.119999 
6.530000 
6.900000 
6.990000 
5.820000 
6.2700C0 
6.400000 

7.500000 
6.580000 
5.7500C0 
6.440000 
5.020000 
5.540000 
6.139999 
6.900000 
6.669999 
6.549999 

7.250000 
6.919999 

5.650000 
5.820000 
6.620000 

5.860000 
6.099999 

6.500000 


5.325740 

-2.485740 

-1.960059 

5.6C9191 

-2.319191 

-1 .828732 

5.427681 

1.282318 

1.011135 

6.032910 

0.1170893 

0.9232739E-C1 

6.334590 

1.525410 

1.202818 

6.336077 

0.8639231 

C.  6812217 

2.160582 

1.839418 

1.450420 

4  .795046 

1  .854954 

1.462671 

5.816327 

-0.2663279 

-0.2100052 

6.127997 

-0.3779974 

-0.2980590 

2.738162 

1.331838 

1.050183 

5  .502283 

-0.5922832 

-0.4670280 

5.272C45 

0.9795475E-G] 

C  •  7723941 E-0 1 

4.390880 

1.019120 

0.8035980 

2.875768 

1.294231 

1.020529 

4  .831560 

0.6684399 

0.5270791 

5.413430 

-0.9343052E-C1 

-0.7367197E-01 

5.502283 

-0.1222839 

-0 .9642351 E-Cl 

1  .646255 

-1.646255 

-1.298107 

4  .526943 

-3.196643 

-2.520858 

5.850915 

0.2690849 

0.2121792 

6.312851 

-1.542851 

-1.216572 

2.817309 

-0.2573090 

-0.2028936 

2.C2C229 

-2.020229 

-1.592994 

1.173855 

-1.173R55 

-0.9256097 

1.131901 

-1.131901 

-0.8925278 

5.571877 

-0.3218775 

-0.2538073 

5.739481 

-1.089481 

-0 .8590794 

5.880660 

-0 .7206602 

-0.5682560 

5.201477 

2.258522 

1.780893 

6.244391 

-1.204391 

-0.9496886 

6.412027 

1.447972 

1.141757 

6.488779 

1.631220 

1.286252 

6.473336 

0.5666351E-01 

0.4468039E-0  1 

6.482731 

0.4172688 

0.3290253 

6.464739 

0.5252609 

0.4141794 

6.448898 

-0.6288986 

-0.4959000 

6.378656 

-0.1086569 

-0 .8567828E-0 1 

6.253687 

0.1463127 

0.1153707 

6.451098 

1.048902 

0.8270812 

6.476042 

0.1039581 

0.8197320E-01 

6.306323 

-0.5563231 

-0.4386727 

6.140313 

0.2996864 

0.2363091 

6.458466 

-1.438466 

-1.134261 

6.476042 

-0.9360418 

-0.7380891 

6.406593 

-0.2665939 

-0.2102150 

6.289943 

0.6100569 

0.4810429 

6.352664 

0.3173351 

0.2502255 

6.476042 

0.7395744E-01 

0.5831702E-01 

6.398646 

0.8513536 

0.6713105 

6.300648 

0.6193514 

0.4883718 

6.335758 

-0.6857586 

-0.5407354 

6.475176 

-0.6551762 

-0.5166205 

6.394775 

0.2252245 

0.1775943 

6.216175 

-0.3561754  -0.2808520 

6.467246 

-0.3672466 

-0.2895818 

6.463087 

0 • 3691 292E- 0 1 

0.291 0662  E-0 1 
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6.4SC0C0 

6.466183 

C. 11816985-01 

0.93179405-02 

6.520000 

6.317667 

0. 2023325 

0.1595435 

6.C700C0 

6.372327 

-0. 3023272 

-0.2383915 

6.C49999 

6.481307 

-Q.4313C7B 

-0.3400955 

6.270000 

6.445546 

-0.1755466 

-0.1364223 

6.049999 

6.327826 

-0.2778263 

-0.2190720 

5.719999 

6.335758 

-0.6157589 

-0.4855391 

7.50999° 

6.481307 

1.028652 

0.8111458 

6.750000 

6.445546 

0.3044538 

0.2400683 

7.099999 

6.270913 

0.8290863 

0.6537522 

5.8C9999 

6.003121 

-0.1931219 

-0.1522806 

6.629999 

6.253039 

0.3769598 

C .29724  C9 

6.369999 

6.283586 

0.1064138 

0.639C959F-01 

6.469999 

6 .22960  3 

0.2403965 

C. 1895578 

5.559999 

5.713383 

-0.1533833 

-0.1209460 

6.099999 

6.143804 

-0.4380417E-C1 

-0. 34540535-01 

6.639995 

6.3047C9 

0.3352900 

C. 2643833 

6.349999 

6.270191 

0.79808245-01 

0.6 2930 49E-01 

6.790000 

5.793820 

0.9961796 

0.7855089 

6.4000CC 

5.748468 

C. 6515312 

0.5137463 

5.609959 

5.680187 

0.1298122 

0.1023597 

6.379999 

5.756392 

0.6236076 

0.4917279 

7.599999 

5.400577 

2.199423 

1.73429 2 

4.440000 

4.782805 

-0.3428059 

-0.2703098 

2.24000C 

4.149226 

-1.909226 

-1.505466 

3.16000C 

2.827332 

0.3326674 

0.2623153 

6.480000 

6.055430 

0.4245691 

0.3347818 

5.7C0000 

5.703501 

-0.3500938E-02 

-0 .2760565  E-02 

3.2300CC 

5.294816 

-2.064S16 

-1.628152 

2.550000 

4.829988 

-2.279989 

-1.797820 

7.110000 

6.469578 

0.6404219 

0.5049863 

5 • 95  0  0  0  C 

6.469578 

-0.5195780 

-0.4096984 

7. 129999 

6.469578 

0.6604214 

0.5207564 

7.200000 

6.469578 

0.7304220 

0.5759534 

6.360000 

6.469576 

-0.1095781 

-0.8640471E-01 

6.5400C0 

6.469578 

0.7C42217E-C1 

0.55529395-01 

5.440000 

6.297028 

-0.8570280 

-0.6757849 

6.410000 

6.255799 

0.1542006 

0.1215904 

6.730000 

6.288687 

0.4413128 

0.3479846 

6.129999 

6.173958 

-0.4395866E-01 

-0.34662355-01 

6.309999 

6.069370 

0.2406292 

0.1897413 

5.740000 

6.255799 

-0.5157995 

-0.4067190 

5.540000 

6 .304756 

-0.7647562 

-0.6030266 

6.750000 

6.250532 

0.4994678 

0.3938410 

3.130000 

5.852849 

-2.722849 

-2.147025 

6.580000 

6.162016 

0.4179840 

0.3295893 

6.900000 

6.285694 

0.6143055 

0.4843930 

5.709999 

6.257086 

-0.5470867 

-0.4313896 

6.419999 

6 .397839 

0.2216053E-C1 

0.1 747405E-01 

3.310000 

5.981867 

-2.671867 

-2.106824 

4.139999 

6.232156 

-2.092156 

-1.649710 

5.419999 

6.138083 

-0.7180834 

-0.5662242 

5.5COOGO 

6.001534 

-0.5015345 

-0.3954707 

5.799999 

5.852165 

-0.52165995-01 

-0.41134005-01 

5.669999 

5.728C03 

-0. 58003435-01 

-0.45736955-01 

5.549999 

5.506139 

0.43860445-01 

0.34584895-01 

4.520000 

5.267131 

-0.7471313 

-0.5891291 
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5.309999  6.064188 
5.639999  5.926G43 
5.4000CC  5.816261 
5.030C0C  5.617781 
4.580000  5.397246 
6.120000  5.971783 
6.120000  5.999149 


0.7541885  -0.5946939 

0.2860432  -0.2255512 

0.4162617  -0.3282313 

0.5877810  -0.4634779 

0.6172464  -0.6444163 

0.1482172  0.1168724 

0.12085C6  0.9529325F-01 


COEF  #  BEST  VALUE  CORRELATION  STD.  ERROR  T 

G  4.872201  0.1532939  0. 4965340  9.812422 

1  0.1456611  0 .3159698E-C1  0.2254287  0.6461514 

2  C. 3264541  0 .6287679E-02  G.10C5615  3.246313 

3  -0.1526691  0.1344724E-02  0 .46 5 0 53 7E -0 1  -3.283257 

4  -C.1517558E-02  0  .24 3 599 4 E - C3  0 . 1 9 79  359 E -0 1 -0 . 766691 7E - C 1 

5  -C.1929272E-C1  0 .2324250E-04  0 .6 1  1 4 035 E -0 2  -3.155480 


6  C.7932583E-02  0 . 6 08 1 74 0 E- 0 5  0 .31 27522E-02  2.536379 

7  0  .2191484E-02  0 . 1 1 63888 E - 0 5  0 . 1 3 68 1 76E -0 2  1.601756 

8  -0.5236040E-C3  0 .1997498E-06  0 .5667998E-03-0. 9237900 

9  C.2638865E-03  0 .5952408E-C8  0 .9 7 84 37 0 E-C 4  2.697021 


ANALYSIS  OF  VARIANCE 

SOURCE  OF  VARIATION  SUM  OF  SQUARES  DOF  MEAN  SQUARE 
TOTAL  902.6241  290 

REGRESSION  450.6854  9  50.07616 

RESIDUAL  451.9365  281  1.608322 

R  SQUARED  =  0.4993058  R  =  0.7066157 

THE  STANDARD  ERROR  OF  Y  ABOUT  THE  REGRESSION  PLANE  IS  1.268196 

THE  F  RATIO  TESTING  THE  NULL  HYPOTHESIS  THAT  ALL  THE 
COEFFICENTS  ARE  ZERO  «  I.E.  B1 =62  =  •  .  .  =BK = 0 «  IS  31.13565 
BASED  ON  THE  DEGREES  OF  FREEDOM  9/  281 


NOU  A  STEP-WISE  DOWN  REGRESSION  WILL  TAKE  PLACE 

THE  VARIABLE  WITH  THE  LOWEST  T  WILL  BE  DELETE 

THIS  WILL  CONTINUE  UNTIL  THERE  ARE  JUST  TWO  VARIABLES  LEFT 


COEF  #  BEST  VALUE  CORRELATION  STD.  ERROR  T 

0  4.895369  0 .9651806E-01  0.3933001  12.44691 

1  0.1345059  0.1B43446E-01  0.1718838  0.7825398 

2  0.3226420  0 .4750544E-02  0 .8725522E-01  3.697682 

3  -0.1522233  0 .1321769E-02  0 .46C2540E-01  -3.307375 

5  -0.1915219E-01  0 . 2 1 1 5 357E - 04  0 .58 22528E-0 2  -3.289325 

6  0 • 7939885 E- 02  0 . 6 0 760 99E- 0 5  0 . 31 2 0 557E-02  2.544381 

7  0.2147755E-02  0 .9616221E-06  0 .1241429E-02  1.730067 

8  -0.5560268E-03  0 .8855436E-0 7  0  .3767252E-03  -1.475948 

9  0.2624116E-03  0 .5722304E-08  0 .9576465E-Q4  2.740172 
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ANALYSIS  OF  VARIANCE 


SOURCE  OF  VARIATION  SUM  OF  SQUARES  DOF  MF AN  SQUARE 
TOTAL  902.6241  2°0 

REGRESSION  450.6760  6  56.33450 

RESIDUAL  451.9461  282  1.602653 

R  SQUARED  =  0.4992953  R  =  C.7C66083 

THE  STANDARD  ERROR  OF  Y  ABOUT  THE  REGRESSION  PLANE  IS  1.265959 

THE  F  RATIO  TESTING  THE  NULL  HYPOTHESIS  THAT  ALL  THE 
COEFFICENTS  ARE  ZERO,  I.E.  B 1  =  B2= .  .  .=BK =0  ,  IS  35.15076 
BASED  ON  THE  DEGREES  OF  FREEDOM  8/  282 


COE F  »  BEST  VALUE  CORRELATION  STD.  ERROR  T 

0  4.987683  0 . 8  78 34 78E  -  0 1  0.3749345  13.30281 

2  0.3255821  0 . 4 74 1 73 6E- 0 2  C.8711457E-C1  3.737401 

3  -0.1187190  0.17797B5E-03  C .  1 6 6 7 74  1 E - C 1  -7.C34195 

5  -0.1946494E-01  0 .21 05390E-C4  0 .58 C4815E -02  -3.353240 

6  0 .6C0  071  7E-02  0  . 2 24 4 5 1 2 E - 0 5  C . 1 8 95 322 E -0 2  3.166067 

7  0. 1933137 E- 02  0 .9146887E-C6  C .12 C9926E-02  1.597732 

6  -C .48881 72E- 03  0 . 6 39 5 1 69E- C 7  0 . 36 6 553 0 E -0 3  -1.333551 
9  0 .2668135E-03  0 .5702560E-08  0 .9 5 53379E -0 4  2.792871 


ANALYSIS  OF  VARIANCE 


SOURCE  OF  VARIATION 

SUM  CF  SQUARES 

DOF 

MEAN  SQUARE 

TOTAL 

902.6241 

290 

REGRESSION 

449.6945 

7 

64.24207 

RESIDUAL 

452.9295 

283 

1.600458 

R  SQUARED  =  0.4962080  R  =  0.7056385 

THE  STANDARD  ERROR  OF  Y  ABOUT  THE  REGRESSION  PLANE  IS  1.265C92 

THE  F  RATIO  TESTING  THE  NULL  HYPOTHESIS  THAT  ALL  THE 
COEFFICENTS  ARE  ZERO,  I.E.  B 1 =B 2= .  .  .  =  BK  =  0 ,  IS  40.13981 
BASED  ON  THE  DEGREES  OF  FREEDOM  7/  283 


COEF  #  BEST  VALUE  CORRELATION  STD.  ERROR  T 

C  4.809273  0 . 7665 1 4 1 E- C 1  0.35C7329  13.71207 

2  C. 3850728  0 .3498266E-02  0 .74 92 778E -0 1  5.139253 

3  -C. 1221848  0 .1737582E-03  0 .1 669895E-01  -7.316917 

5  -0.2425978E-01  0 . 1 297624 E- 0 4  0 .45 63425E-02  -5.316135 

6  0.7190097E-02  0 . 1  74 74 8 7E- 0 5  C . 1 6 74648E -02  4.293498 

7  0.4305980E-03  0 . 1 2 1 4 795E -0 6  0 .44 15379E-03  0.9752233 

9  0 .3456187E-03  0 .3520603E-08  0 .75 16661 E-04  4.598035 


ANALYSIS  OF  VARIANCE 


SOURCE  OF  VARIATION 

SUM  OF  SQUARES 

DOF 

MEAN  SQUARE 

TOTAL 

902.6241 

290 

REGRESSION 

446.8484 

6 

74.47472 

RESIDUAL 

455.7758 

284 

1.604644 

R  SQUARED  =  0.4950546  R  =  0.7036013 


THE 

STAN 

DARD  ERROR  OF  Y 

ABOUT 

THE  REGRE 

SS  ION 

PLANE  IS 

1  .266824 

THE 

F  RA 

TIO  TESTING  THE 

NULL 

HYPOTHESIS 

THAT 

ALL  THF 

COEFFICENTS  ARE  ZERO* 

I.E.  B1=B2=. 

..=BK=C,  IS 

46.40620 

BASED 

ON  THE  DEGREES 

OF  FREEDOM 

6/  284 

COEF 

#  BEST  VALUE  COR 

RELATION  STD. 

ERROR  T 

C 

4.806970 

0 .7664793E-01 

0.3506948 

13.70699 

2 

0.3756840 

0 .3440512E-C2 

0.743C031E-01 

5.056291 

3 

-0.1188613 

0.1665213E-03 

C.1634610E-01 

-7.271542 

5 

-0 .232310 1 E- 01 

0.1228282E-04 

C .4439439E-02 

-5.232871 

6 

0 .7416429E-02 

0 .1 713925E-05 

0 .1658345E-02 

4.472186 

9 

0.3276729E-C3 

0 .3309602E-08 

0.7287305E-04 

4.496490 

ANALYSIS  OF  VARIANCE 


SOURCE  OF  VARIATION 

SUM  OF  SQUARES 

TOTAL 

902.6241 

REGRESSION 

445.3221 

RESIDUAL 

457.3019 

DOF  KEAN  SQUARE 
290 

5  89.06441 

285  1.604568 


R  SQUARED  =  0.4933639 
THE  STANDARD  ERROR  OF 


R  =  0.7023986 

ABOUT  THE  REGRESSION  PLANE  IS  1.266716 


THE  F  RATIO  TESTING  THE  NULL  HYPOTHESIS  THAT  ALL  THE 
COEFFICENTS  ARE  ZERO,  I.E.  B 1  =  62= . .  •  =  BK  =  0  •  IS  55.50678 


BASED 

ON  THE  DEGREES 

OF  FREEDOM 

5/  285 

COEF 

U  BEST  VALUE  CORRELATION  STD. 

ERROR  T 

a 

4.423452 

0.7206470E-01 

0.3511621 

12.59661 

2 

0.4161844 

0 .33894 C0E-02 

0.7615663E-01 

5.464848 

3 

-0 .4793454E-01 

0 .9765912E-05 

0.4087923E-02 

-11.72589 

5 

-0.2572381E-01 

0.1206918E-04 

0.4548255E-02 

-5.655754 

9 

0 .3682456E-03 

0 .3258307E-08 

0.7466934E-04 

4  .931684 

ANALYSIS  OF  VA 

RI ANCE 

SOURCE  OF  VARIATION 

SUM  OF  SQUARES 

DOF  MEAN  SQUARE 

TOTAL 

902.6241 

290 

REGRESSION 

413.2301 

4  103.3075 

RESIDUAL 

489.3940 

286  1.711168 

R  SQUARED  =  0.4578C97  R  =  0.6766163 

THE  STANDARD  ERROR  OF  Y  ABOUT  THE  REGRESSION  PLANE  IS  1.308116 

THE  F  RATIO  TESTING  THE  NULL  HYPOTHESIS  THAT  ALL  THE 
COEFFICENTS  ARE  ZERO,  I.E.  B 1  =  B2  =  .  .  .  =BK  =  0  ,  IS  60.37251 
BASED  ON  THE  DEGREES  OF  FREEDOM  4/  286 
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T 


COEF  ft  BEST  VALUE  CORRELATION  STD.  ERROR 

0  5.702243  0.3C27507E-01  0.2366758 

2  C.77255OOE-01  0 .6292002E-03  0 .3012083E-01 

3  -0 .5168957L-01  0 . 90 2 7 1 1 2E - 0 5  C .0 1 76386E -02 

5  -0 .3780753E-02  0 .5239919E-06  0 .9806312E-03 


ANALYSIS  OF  VARIANCE 

SOURCE  OF  VARIATION  SUM  OF  SQUARES  DOF  MEAN  SQUARE 
TOTAL  902.6201  290 

REGRESSION  371.6119  3  123.8706 

RESIDUAL  531.0122  287  1.850217 

R  SQUARED  =  0.0117017  R  =  0.6016398 

THE  STANDARD  ERROR  OF  Y  ABOUT  THE  REGRESSION  PLANE  IS  1.36C227 

THE  F  RATIO  TESTING  THE  NULL  HYPOTHESIS  THAT  ALL  THE 
COEFFICENTS  ARE  ZERO*  I.E.  B1  =  B2  = .  .  ,=BK =0  *  IS  66.90923 
BASED  ON  THE  DEGREES  OF  FREEDOM  3/  287 


COEF  #  BEST  VALUE  CORRELATION  STD.  ERROR  T 

0  6.227016  0.5498795E-02  0.1015861  61.29792 
3  -0.5185799E-01  0 . 90 20 1 21 E- 0 5  0 .02 05532E-02  -12.33090 
5  -0.16920H2E-02  0 . 62 29 1 62E- 0 7  0 .30 19 128E-03  -0.908870 


ANALYSIS  OF  VARIANCE 

SOURCE  OF  VARIATION  SUM  OF  SQUARES  DOF  MEAN  SQUARE 
TOTAL  902.6201  290 

REGRESSION  362.1267  2  181.0634 

RESIDUAL  500.0973  288  1.876727 

R  SQUARED  =  0.0011933  R  =  0.6333982 

THE  STANDARD  ERROR  OF  Y  ABOUT  THE  REGRESSION  PLANE  IS  1.369937 

THE  F  RATIO  TESTING  THE  NULL  HYPOTHESIS  THAT  ALL  THE 
COEFFICENTS  ARE  ZERO,  I.E.  B 1  =  B2  =  . . . =BK =0  ,  IS  96.47827 
BASED  ON  THE  DEGREES  OF  FREEDOM  2/  288 

COEF  ft  BEST  VALUE  CORRELATION  STD.  ERROR  T 

0  5.983677  0 . 02 1 0 50 7E- 02  0 .9203528E-01  60.73369 

3  -0 .5022976E-01  0 . 93 0 1  73 5E- 0 5  0 . 03 00628E-02  -12.08203 


ANALYSIS  OF  VARIANCE 

SOURCE  OF  VARIATION  SUM  OF  SQUARES  DOF  MEAN  SQUARE 
TOTAL  902.6201  290 

REGRESSION  316.1632  1  316.1632 

RESIDUAL  586.0608  289  2.029276 

R  SQUARED  =  0.3502712  R  =  0.5918372 

THE  STANOARD  ERROR  OF  Y  ABOUT  THE  REGRESSION  PLANE  IS  1.024527 

THE  F  RATIO  TESTING  THE  NULL  HYPOTHESIS  THAT  ALL  THE 
COEFFICENTS  ARE  ZERO,  I.E.  B1  =  B2  =  . . .=BK =0 ,  IS  155.8010 
BASED  ON  THE  OEGREES  OF  FREEDOM  1/  289 


20 .26206 
2.260172 
-12.37663 
-3.80382R 


55 


b 


i 


W 
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«  OF  VARIABLES 


APPENDIX  D:  REGRESSION  ANALYSIS  WITH  TWO  COEFFICIENTS 
OF  THE  ORIGINAL  POLYNOMIAL  DELETED 

?C1 

8 


C  ATA 

FORMAT  =  ( 

8  0  13*5) 

THE 

MATRIX  OF  CCRREC 

TEC  SUMS  CF 

SGLARFS  AND  CRCS 

S  PRODUCTS 

VAR  1 

VAR  2 

VAR  3  VAR 

4  VAR 

5 

VAR  6 

VAR  7 

V  Afi  8 

VAR 

1 

13506. G7 

VAR 

2 

3847.509 

107506 .8 

VAR 

3 

440245.8 

152691.2 

C.1626475F  08 

VAR 

4 

37844.24 

1022243. 

1454390. 

0.1C31309E 

CF 

VAR 

5 

233512.9 

1462337. 

7753955. 

C.1416536E 

C  £ 

C .3 122043E 

08 

VAR 

6 

1109542. 

0 . 279  5  25 4  E  09 

3080069. 

0.3667886 E  08 

0 • 28  4 126  3E 

08 

C.8400481E 

08 

VAR 

7 

0.1375857E  C8 

434C3C6. 

C.5622751E  09 

0.4207228E 

C  8 

C.22C2168E 

09 

0. 1042  44 0E  10 

0  •  2109822E 

11 

VAR 

8 

-821.6802 

-5830.070 

-35335.82 

-50474.70 

-85460.37 

-243978.2 

-1176479. 

902 .6241 

THE 

THE 


VAR 

VAR 

VAR 

VAR 

VAR 

VAR 

VAR 

VAR 


THE 

RIJ 


INVERSE  OF  THE  ABOVE  MATRIX 
GAUSSIAN  MULTIPLIERS  -  CIJ 

VAR  1  VAR  2  VAR  3  VAR 

VAR  6  VAR  7  VAR  8 

1  G.4975777E-02 

2- 0.5676357E-04  0 .2G90964E-03 

3- G.3193991E-03  0  .  5 56 756 OE - 05  0  .2 1 £9 0 4 2 E -C 4 
4  0.1842201E-04-0.1891941E-04-0.1532943E-C5 
5-0 .2736755E-C4  0 , 3853 1 6 0 E - C £  0.2157947E-05- 

,4fc71C66E-06-C.8024814r-06 


VAR 


0.2324013E-C5 

0.6661299E-C6 


6  0.9865832E-05-0 
0  .B447924E-07 

7  0.4999476E-05-0.82C783CE-07-C.3535946E-06 
0.1324641E-07  0 .5859736E-08 

8-0.7188362E-03  0  .2621137E-03  0  .42 9 7S65E -04  -  0 . 1 32 4868E - 0 4 
0. 1 079  235 E-05-0. 5S9 084  IE-06  0 .22 C 7849E-C2 


C  .9  229395E-06 
0.1977934E-06-0.2601391E-06 
0.2442903E-C7-C.354C064E-07 
C.42te074E-05 


SAMPLE  PARTIAL  CORRELATION  COE F F ICE N T S  IR  I J 
=  CIJ  /  (  CII  *  CJ'J  >**.5 

THE  VARIABLES 


VAR 


THESE  REPRESENT 

CORRELATIONS 

BETWEEN 

THE 

CLOSER  TC  4- 

1  « 

THE  MORE 

DEPENDENCE 

THE 

CLOSER  TO  0 

THE 

MORE  INDEPENDENCE 

VAR  1 

VAR  2 

VAR  3 

VAR  6 

VAR  7 

VAR  8 

VAR 


VAR  1  1.000000 

VAR  2-0.5565010E-01  1.000000 


VAR 

3 

-0.9677791 

0.8229331E-01 

1.0CC000 

VAR 

4 

0.1713118 

-0.8582521 

-C. 2149214 

1.000000 

VAR 

5 

-0.4038488 

0.2773685E-01 

C.48C0948 

-0.4548341 

l.COOOOO 

VAR 

6 

0.4812029 

-0.1111393 

-0.5901099 

0.4463931 

-0.9316313 

1.000000 

VAR 

7 

0.9258808 

-0.7415084E-C1 

-C. 9872782 

0.2093379 

-C. 4813772 

0 .5953661 

1.000000 

VAR( 

8 

-0.2168776 

0.3857729 

0.1954837 

-0.1849560 

-0.9454981E-01 

0.79 02349 E-01 

-0.1637773 

1.000000 

THE 

STANDARD  ERRCRS 

CF  THE  ABCVE  R 

Iw  I  S  1  V 

SIJ 

2 

<  (  1  -  RIJ*RIJ  >  /  OOF  RES 

IOUAL  l**.5 

VAR  1 

VAR  2  V 

AR  3  VAR 

4  VAR 

5 

VAR  6 

VAR  7  V 

AR  8 

VAR 

l 

o.ooocooc 

VAR 

2 

0 .5935171E-01 

0.0000000 

VAR 

3 

0.1496799E-01 

0 .5924221E-01 

0.0000000 

VAR 

4 

0 .5856506E-01 

0.3050825E-01 

0 .58C5471E-01 

0.0000000 

VAR 

5 

0.5438075E-01 

8.5942096E-01 

0 .52 14511E-01 

0.5293926E-C1  C.OCOOOOO 

VAR 

6 

0.5210901E-01 

0.5907557E-01 

0.4799037E-01 

0.5319254E-01  0 .2 160218E-01 

0.0000000 
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1*1 


VAR  7  0  •  2245882E- 01  0 . 592 80 1 8E - 0 1  0  .945 1695E -02  0 .5812675E - C 1  C . 5 2 1 G 332E-0 1 
0.4776041E-01  C.CCOOOnC 

VAR  8  0.5802900E-01  0 . 5  4 8 4  2 5 1  E- 0 1  G  .582S69B*'-0  1  C  .58  4  1823E- C  1  C . 5 5 1 7 753E- 0 1 
0 .5925793E-C1  0 . 5 864  1 1 8E  -  0 1  O.COOGCOO 


1  ACTUAL 

PREC1 CTLD 

RFSIDUAL  ST. 

RESIDUAL 

2.300000 

4.304292 

-2.004292 

-1.584306 

5.000000 

6.222622 

-1.222622 

-0.9664294 

5.48000C 

6.478C02 

-0.9980  C2C 

-0.7888778 

5 . 799993 

6 .557627 

-0.7576275 

-C. 5988716 

6.160000 

6.507816 

-0.3478  165 

-0.2749238 

4  .  P  09999 

6.5C8A66 

-1.698468 

-1.242565 

6.200000 

6  .5  69  00  0 

-0.3690004 

-0.2916788 

6.200GSC 

6.569CCC 

-0.3690004 

-C .2916788 

6.690000 

6.501635 

P.1883650 

0.1488943 

5.879999 

5.264520 

0.6154795 

0.4865096 

1.270000 

3.581354 

-2.311254 

-1.827024 

A  1.260000 

A .529723 

-2.269722 

-2.584574 

6.000000 

5.156713 

C. 8432865 

0.6.665812 

6.129999 

5.388A37 

0.7415619 

0.5861722 

6.48000C 

5.651977 

0.8280230 

C. 6545160 

6.679999 

5.299209 

1.280751 

1.091455 

7.919999 

5 .85952  A 

2.060 A75 

1 .628716 

7.009999 

6.112702 

0 .8972569 

0.7092741 

8.109999 

6.165216 

1.9AA782 

1.527265 

6.200000 

6.258952 

-0.5895233E- 

01-C.4655925F 

6.169999 

6.308202 

-C.l?e2027 

-0.1092432 

6.000000 

6.242105 

-0.2421055 

-0.1913728 

4.620000 

3.961247 

C.6586523 

0.5206266 

4.620000 

2.961247 

0.6586523 

C.52C6266 

4.809999 

4 .632504 

0.1774960 

0.1403028 

5.480000 

4.967206 

0.5127535 

0.A053409 

5.480000 

4.967206 

0.5127525 

0.4053409 

4.250000 

4.341299 

-0.9139919E- 

C1-0.7224709F 

4.250000 

4.341295 

-0.9139S19E- 

01-0,72247  09  F. 

_  6.129999 

5.072552 

1.057447 

C. 8358660 

6.129999 

5.072552 

1.057447 

0.8358660 

5.820000 

5.552572 

0.2270265 

0.1794549 

5.490000 

5.795798 

-0.3057985 

-0.2417204 

5.490000 

5.795798 

-0.3057985 

-0.2 417204 

5.330000 

5 .466882 

-0.1268818 

-0.1CR1991 

6.660000 

5.848972 

0.8110275 

C. 6410818 

6.660000 

5.981642 

0.6783581 

0.5362124 

5.000000 

5.908144 

-0.9081440 

-0.7178483 

5.5P0C0C 

5.768806 

-0.1888065 

-0.1492434 

C  5.490C0C 

5.823181 

-0.3331814 

-0.2633654 

5.469999 

5.875707 

-0.4057074 

-0.3206940 

3 .260  COO 

3.805850 

-0.5458503 

-0.4314709 

3.300  CO  0 

5.659334 

-2.399324 

-1.896569 

3.300  OG  0 

5.699234 

-2.399224 

-1.856569 

4.459999 

5 .59462  C 

-1.134621 

-0.8968682 

4.459999 

5.594620 

-1.134621 

-0.89686 82 

3.910000 

5.760390 

-1.850390 

-1.462653 

3.9100C0 

5.76C29C 

-1.850350 

-1.462653 

3.840000 

5.560365 

-1.720365 

-1.359874 

3.840000 

5 .560365 

-1.720265 

-1.359874 

5.179999 

5.566005 

-0.3860054 

-0.3051205 

5.179999 

5.566CCS 

-0.3860054 

-0.3C51205 

3.350000 

5.980424 

-2.630424 

-2.079236 

3.350030 

5.980424 

-2.630424 

-2.079236 

6.059999 

5 .779962 

0.2800360 

0.2213562 

6.C59999 

5.779962 

0. 280026  C 

0.221356 2 

3.910000 

6.145840 

-2.235840 

-1.767334 

3.9100C0 

6.145840 

-2.235840 

-1.767334 

6.059999 

5.959550 

0.1C04496 

0.7940099E 

6.059999 

5 .9  5955  C 

0.1004456 

0.7940099E 

3.910000 

6.258507 

-2.348507 

-1.856392 

3.910000 

6 .258507 

-2.348507 

-1.856392 

5.190000 

6 .097091 

-0.907051] 

-0.717C160 

5.19C0C0 

6.C97C91 

-0.9C70911 

-C. 7170160 

3.660000 

5.767786 

-2.107786 

-1.666113 

3.660000 

5.767786 

-2.107786 

-1.666113 

6.969999 

5.604772 

1.265228 

1.079153 

6.969999 

5.604772 

1.365228 

1.079153 

58 


6.520030 

6.2  15G23 

0.3C99765 

0.241C706 

6.520000 

6.215023 

0.3099765 

0.2410706 

5.49000C 

6.098199 

-0.5581592 

-0.4412282 

5.990000 

6.098199 

-0.5581592 

-0.441 2282 

5.990000 

6 .098199 

-0.5581592 

-0.4412282 

6.209999 

5.920282 

0.7897167 

0.6242366 

6.209999 

5  .920282 

0.7897167 

0.6242366 

£  6.209999 

5  .920282 

0.7897167 

0.6242366 

7.929999 

6.259801 

1.135155 

0.8973249 

7.929999 

6 .299801 

1.135199 

0.8973249 

7.929999 

6.299801 

1.135199 

0.8973249 

7.929999 

6  .259801 

1.135159 

0.8973249 

7.929959 

6.259801 

1.135155 

C .8973249 

6.580000 

5.731690 

0.8983601 

0.6705916 

6.580000 

5  .731690 

0.8983601 

0.6705916 

7.929999 

6.916529 

1.C13C75 

0.8007914 

7.929999 

6.916529 

1.013C75 

C.8CC7914 

7.929999 

6  .916929 

1.013075 

0.8007914 

7.929999 

6.916929 

1.013C75 

0.8007914 

6.090000 

6.922196 

-0.3821569 

-0.3021097 

6.090000 

6.922196 

-0.3821569 

-0.3021097 

6.09000C 

6 .922196 

-0.3821969 

-0.3021097 

6.520000 

6  .325626 

0.1993712 

0.1536420 

6.520000 

6.325628 

0.1993712 

0.1536420 

6.520000 

6.325626 

0.1993712 

C. 1536420 

9.879999 

5  .939689 

-1.059685 

-0.8336823 

9.879999 

5 .939689 

-1.C59685 

-0.8336823 

5.799999 

2.912790 

3.387209 

2.677441 

3.980CC0 

5.118957 

-1.638957 

-1.295129 

C  6.959999 

5.777370 

1.182629 

0.9348162 

r  3.980000 

5 .999771 

-2.019771 

-1.592589 

3.9800C0 

5.999771 

-2.019771 

-1.592589 

3.980000 

5.959771 

-2.019771 

-1.552589 

5  7.309999 

5.030637 

2.279363 

1.801737 

6.580000 

6.373537 

0.2C69629 

0.1631999 

6.559999 

6.293627 

0.3161 726 

0.2499206 

7.200000 

E.998C86 

1.201519 

0.9500604 

6.780000 

5.586158 

1.193892 

0.9436799 

3.76000C 

9.958005 

-1.198005 

-0.9469709 

9.639999 

9.909923 

0.2355766 

0.1862130 

6.519999 

6.109731 

0.8102689 

0.6404817 

9.590000 

3.066562 

1.973938 

1.16*689 

U  5.959999 

5.919582 

0.4C41672E- 

Cl  0.3194765E 

^  6.750000 

5.993999 

0.8060064 

0.6371129 

5.139999 

5.983812 

-0.3438125 

-0.2717692 

3.090000 

9.975507 

-1.935507 

-1.529934 

6.910000 

5.915832 

0.9941683 

0.7858466 

3.600Q00 

9 .996951 

-1.396451 

-1.103834 

1.630000 

3.9C9931 

-2.279931 

-1.802186 

8.539999 

5.515058 

3.029941 

2.391084 

9.370000 

3.179690 

1.190310 

0.9408877 

5 . 089999 

9.957758 

0.1322412 

0.1045309 

5.759999 

5.355609 

0.4041548 

0.3194984 

I  5.950000 

5.929920 

0.2507973E- 

01  0 • 1982443E 

9.969999 

5.390929 

-0.3704243 

-0.2928042 

2.92C000 

3.086395 

-0.1663557 

-0.1315285 

t  9.679999 

9.952165 

-0.2721701 

-0.2151386 

•J  6.9900C0 

5.629003 

0.8109970 

0.6410577 

8.779999 

5 .779058 

3.000541 

2.372113 

7.099999 

6.339512 

0.7109874 

0.5616091 

7.900000 

6.392397 

1.057652 

0.8360281 

7.099999 

6.162768 

0. 8872309 

0.7013172 

8.759998 

6.268029 

2.491969 

1.969793 

7.610000 

6.183690 

1.426359 

1.127475 

7.690000 

5.96289C 

1.727110 

1.365205 

7.969999 

5.625905 

1.844094 

1.457676 

6.090000 

9.950609 

1.589396 

1.256348 

6.250000 

6.537173 

-0.2871733 

-0.2269979 

5.879599 

6.509018 

-0.6240187 

-0.4932595 

8.129999 

6.939383 

1.695616 

1.340310 

8.129999 

6.910065 

1.719934 

1.359533 

8.500000 

6.291588 

2.208412 

1.745653 

10.25000 

6.108129 

4.141876 

3.273972 

7.870000 

6.066779 

1.803221 

1.425367 

7.000000 

5.892526 

1.107474 

0.8754101 

59 


► 

6.879999 

5.834808 

1.045151 

0.8261777 

- - 1 

r\ 

3.61000C 

5 .085202 

-1.275202 

-1.007992 

3.810000 

4.584808 

-0.7748C84 

-0.6124523 

4.500000 

4  .584808 

-0.8480835E- 

01-0.6703731 E-Gl 

■  .* 

4.500000 

3.577902 

0.9220977 

0.7288779 

4.440000 

2.902813 

1.537106 

1.215078 

2.620000 

2  .335700 

0.2802557 

0.2215647 

11 

5.940000 

2.333645 

3.606250 

2.850662 

6.500000 

5.498426 

1.Q01 574 

0.7917001 

8.000000 

4.557915 

3.442CP5 

2.720818 

3.5COOOC 

3  .653275 

-0.1933794 

-0.1528580 

L 

1.750000 

4.596592 

-2.846552 

-2.250107 

wv- 

3.000000 

3.241274 

-0.2412729 

-0.1907165 

w\s 

1.875000 

2.09337C 

-0.2183700 

-0.1726120 

V 

3.500000 

3.026183 

0.4738165 

0.3745315 

4.000000 

3.026183 

0.9738169 

0.7697597 

1.500000 

2.370481 

-0.8704815 

-0.6880777 

0.8750000 

1.472952 

-0.5979524 

-C. 4726553 

» '• 

2.250000 

0.5111086 

1.238051 

1 .058335 

0.7500000 

0.4908113 

0.2591887 

0.2048773 

1 

0.0000000 

0.6249524 

-0.6249524 

-0.4939976 

r_  • 

J).  0000000 

0.4501428 

-0.4501428 

-0.3558183 

i 

5.849999 

6.3  06835 

-0.4568396 

-0.3611118 

7.250000 

5.548052 

1.701548 

1.345316 

L 

8.5000C0 

5.172577 

3.326423 

2.629392 

2.050000 

4.069414 

-2.019414 

-1.596259 

3.970000 

4.346027 

-0.3760276 

-0.2972334 

2.840000 

5.282257 

-2.442257 

-1.930498 

290000 

5.655828 

-2.365828 

-1.870084 

6.709599 

5  .373402 

1.336597 

1.C56522 

r 

M 

6. 150C00 

5.S81842 

0.1681576 

0.1329212 

7.860000 

6.3C2876 

1.557124 

1.230839 

7.200000 

6.323264 

0.8767357 

0.6930213 

4.000000 

2.270475 

1.729521 

1.367111 

6.65000G 

4.752628 

1.897271 

1.459789 

5.549999 

5.764395 

-0.2143955 

-0.1694703 

-• 

5 •  75  0  0  0  C 

6.233960 

-0.4839602 

-0.3825494 

• 

4.070000 

2.786915 

1.283085 

1.014222 

r 

4.51C000 

5.615028 

-0.7050285 

-0.5572544 

5.370000 

5.223456 

0.1465435 

0.1158362 

• 

5.410000 

4.337351 

1.072649 

0.84788 22 

4.169999 

2.909764 

1.260225 

0.9961606 

5.5C0C00 

4.766C63 

0.7239273 

0.5801454 

* 

H 

5.320000 

5.383271 

-0.6327 152E- 

01-0.5001338E-01 

5.379999 

5.615020 

-0.2350252 

-0.1857803 

0.0000000 

1  .731690 

-1.731690 

-1 .368825 

1.330000 

4.513255 

-3.183259 

-2.516307 

6.120000 

5.791569 

0.3284311 

0.2596105 

. 

4.77000C 

6.414475 

-1.644480 

-1.299890 

■ 

2.560000 

2.862C7C 

-0.3020697 

-0.2387729 

C. 0000000 

2.127C31 

-2.127C31 

-1.681326 

0.0000000 

1  .271621 

-1.271621 

-1.005161 

' 

0.0000000 

0.9473299 

-0.9473299 

-0.7488230 

1 

5.25000C 

5.684145 

-0.4341488 

-0.3431757 

4.650C0C 

5  .822532 

-1.172522 

-0.9268355 

5.160000 

5.990149 

-0.8301497 

-0.6561972 

*  ‘  - 

7.459999 

5.157095 

2.302500 

1.820342 

5.040000 

6.300698 

-1.260698 

-0.9965271 

7.86CCOO 

6.441965 

1.418C25 

1.120894 

8.119999 

6.544386 

1.575613 

1.245453 

6.530000 

6  .533924 

-0.3924370E- 

C2-0.3102043F-02 

6.900000 

6.539552 

0.3600473 

0.2846017 

*4 

6.990000 

6.549567 

0.4404225 

0.3481427 

5.820000 

6.543996 

-0.723996 2 

-0.5722874 

r— 

6.270000 

6.477884 

-0.2076848 

-0.1643239 

6.400000 

6.356975 

0.4302502E- 

Cl  0.3400940E-C1 

7.500000 

6.534182 

0.9656184 

0.7634373 

6.580000 

6.563098 

0.169019 7E- 

01  0.1336027E-01 

5.750000 

6.408152 

-0.6581516 

-0.5202402 

6.440000 

6.246C51 

0.1939487 

0.1533080 

4 

5.020000 

6.542422 

-1.522424 

-1.203410 

5.540000 

6.563098 

-1.023098 

-0.8087144 

6.139999 

6.50449C 

-0.36449C5 

-0.2881139 

,  -‘ 

6.900000 

6.392256 

0.5077438 

0.40134c3 

‘  ‘ 

6.669999 

6.428136 

0.2418623 

0.1911823 

60 

— 

T - T 


r-'.  -  w 


R ,,,, 


^ . 

tr? 

6. £49999 

6.563:98 

-0.13  C9872E- C 1 

-C • 1 C35397E -Cl 

7.250000 

6.496949 

0.7530508 

0.5952537 

6.919999 

6.402647 

0.5173521 

C  ->4089442 

5.650000 

6.365749 

-0.7157457 

-0.5657690 

V1. 

5.820000 

6  .535571 

-0.7155714 

-0.5656281 

6.620003 

6.465433 

0.1545668 

0.1221783 

5.660000 

6.297429 

-0.4374254 

-0.3457689 

% 

6.099999 

6  .513016 

-0.4130182 

-0.3264729 

,* 

6.500000 

6.524581 

-0.2458C56E-C1 

-C.1943C18F-01 

P.'- 

6.480000 

6.528156 

-0.4815674E-01 

-0.3806581E-01 

6.520000 

6.393055 

0.1269407 

0.1003411 

6.070000 

6  .406540 

-0.3365402 

-0.2660204 

6. G49999 

6  .528633 

-0.4788342 

-0.3784975 

6.270C00 

6.509317 

-0.2393179 

-0.1891704 

6.049999 

6.402545 

-0.35295C1 

-0.2789917 

| 

5.719999 

6.365745 

-0.6457500 

-0.5104373 

■ 

7.5C9999 

6.528632 

0.9811655 

0.7755687 

6.750000 

6 .509317 

0.2406826 

0.1902491 

-  . 

7.099999 

6.348620 

0.7513790 

0.5939323 

5.P09995 

5  .947769 

-0.1377657 

-0.1089009 

*  •• 

6.629999 

6.209372 

0.4206276 

0.3324877 

o 

6.389999 

6.246623 

0.1433764 

0.1133327 

»  »' 

6.469999 

6.216C37 

0.2539625 

0.2007463 

i 

5.559999 

5  .657181 

-0.9718132E-01 

-0.768176CE-01 

6.099999 

6 .05251  5 

0.7084846E-C2 

0.560G262F-C2 

6.639999 

6.272431 

0. 3675680 

0.2905465 

6.349999 

6  .255570 

0.9442902E-01 

0 .7464200E -01 

6.790000 

5  .734200 

1.05.5759 

0.8345634 

6.4C00CC 

5.689133 

0.7108669 

0 .5619092 

. 

5.809999 

5.638890 

0.1711092 

0.1352544 

3 

6.379999 

5.726476 

0.6535234 

0.5165817 

7.599999 

5.346547 

2.253452 

1.781256 

4.440000 

A. 748641 

-0.3086414 

-0.2439676 

2.240000 

4.158751 

-1.918751 

-1.516689 

3.160000 

2.916095 

0.2439051 

0.1927963 

6.48C000 

5.998654 

0.4813452 

0.3804823 

5.700000 

5.646690 

0.5330544E-C1 

0.421 38 79F-01 

3.230000 

5.253700 

-2.023701 

-1.599647 

2.550000 

4.815727 

-2.265728 

-1.790959 

1 

7.110000 

6.527246 

0.5827541 

0.4606417 

5.950000 

6.527246 

-0.5772457 

-0.4562876 

7.129999 

6.527246 

0.6027536 

0.4764504 

7.200000 

6.527246 

0.6727543 

0.5317829 

6.360000 

6.527246 

-0.1672459 

-0 . 1322006 

6.540000 

6.527246 

0.1275444E-C1 

0.1008183E-C1 

v 

5.440000 

6.265694 

-0.8256941 

-0.6526752 

6.410000 

6.212502 

0.1974974 

0.1561130 

■ 

6.730000 

6.259886 

0.4701138 

0.3716044 

(P* 

6. 129999 

6.168585 

-0.3858566E-C1 

-0.30500 29 E-01 

*  _ 

6.309999 

6.016543 

0.2934561 

0.2319642 

v  ■ 

5.740000 

6.212502 

-0.4725C27 

-0.3734928 

5.540000 

6.279718 

-0.7397184 

-0.5847152 

6.750000 

6.240209 

0.5097914 

0.4029678 

3.130000 

5.796247 

-2.666247 

-2.10755 2 

•*. 

:  r 

6.580000 

6.110955 

0.4690447 

0.3707593 

6.9000C0 

6.249181 

0.6508188 

0.5144439 

5.709999 

6.238689 

-0.5286SC3 

-0.4179067 

6.419999 

6.392898 

0.2710152E-C1 

0.2142256E-01 

3.310000 

5.926182 

-2.616182 

-2.067978 

" 

_4. 139999 

6.187001 

-2.047002 

-1.618066 

m  . 

5.419999 

6.085236 

-0.6652365 

-0.5258405 

5.500000 

5.943583 

-0.4435835 

-0.3506334 

ft 

5.799999 

5.793288 

0.6711006E-C2 

0.5304758E-02 

■# 

5.669999 

5.670807 

-0.8077621E-03 

-0 .63850 08E-03 

5.549999 

5  .4  5569  8 

0.9430122E-C1 

0.7454100E-01 

4.520000 

5.227554 

-0.7075548 

-0.5592912 

p 

5.309999 

6.007665 

-0.6976652 

-0.5514740 

5.639999 

5.866920 

-0.2269211 

-0.1793712 

•\ 

5.400000 

5.757475 

-0.3574791 

-0.2825717 

5.030000 

5.563276 

-0.5332766 

-0.4215319 

*  ■ 

4.580000 

5.351203 

-0.7712030 

-0.6096025 

6.120000 

5.913620 

0.2063799 

0.1631343 

6.120000 

5.940900 

0.1791000 

0.1415707 

*. 

COEF  ft  BEST 

VALUE  CORRELATION  STD.  ERROR 

T 

6! 


0 

1 

2 

3 

4 
c, 

f 

7 


4.987683 
0.3255821 
0.1187190 
0 .1946494E-01 
0 .6000717E-02 
0 .1933137E-02 
0.4888172E-03 
U .2668135E-03 


0.8783478E-01 
0 .4741736E-02 
0.1779785E-03 
O.210539OE-O4 
0.2244512E-05 
0 .S146887E-C6 
0  •  8395’169E-07 
0 .5702560E-08 


0.3749345 
0 .8711457E-01 
0.1687741E-01 
C .58  C  48 15E-C2 
0 .1895322E-02 
0.12C9926F-02 
C.3665530E-03 
C.9553379E-04 


13.30281 

3.737401 

-7.034195 

-3.35324C 

3.166067 

1.597732 

-1.333551 

2.792871 


ANALYSIS  OF  VARIANCE 

SOURCE  OF  VARIATION  SUM  CF  SQUARES  CCF  M  E  AN  SQUARE 

TCTAL  902.6241  290 

REGRESSION  449.6945  7  64.24207 

RESIDUAL  452.9295  283  1.600458 

R  SGUARED  =  0.4982080  R  =  0.7058385 

THE  STANDARD  ERRCR  OF  Y  ABOUT  THE  REGRESSION  PLANE  IS  1.265C92 

THE  F  RATIO  TESTING  THE  NULL  HYPOTHESIS  THAT  ALL  THE 
COEFFICENTS  ARE  ZERO.  I.E.  b  1  =  B 2 =  .  .  .  =  F*K  =  0  .  IS  40.13981 
EASED  CN  THE  DEGREES  OF  FREEDOM  7/  283 


NOW  A  STEP-WISE  DOWN  REGRESSION  WILL  TAKE  PLACE 

THE  VARIABLE  WITH  THE  LOWEST  T  WILL  EE  DELETE 

THIS  WILL  CONTINUE  UNTIL  THERE  ARE  JUST  TWO  VARIABLES  LEFT 


COEF 

C 

1 

A 

<- 

3 

4 

5 

7 


#  BEST  VALUE  CORRELATION  STD.  ERROR  T 

4.809273  0 .7665141E-01  0  . 3507329 

0.3850726  0 . 3 4982 66 E- 02  0  .74 9 2 778 E - 0 1 

-0.1221848  0.1 737582E-03  0 .1669895E-01 

-C. 24259781-01  0 .1297624E-04  0  .4563425E-02 
0.7190097 E- 02  0 . 1 747487E- 0 5  0 .1674648E-02 
0 • 43  0598  0  E- 03  0 . 1 2 1 4 795E- 06  0 .44 15379E-03 
0.3456 18 7E-03  0 .3520603E-08  0  .75 1666 1 E -0 4 


13.71207 

5.139253 

-7.316917 

-5.316135 

4.293498 

0.9752232 

4.598035 


ANALYSIS  OF  VARIANCE 

SCURCE  OF  VARIATION  SUM  CF  SQUARES  CCF  MEAN  SGUARE 
TOTAL  902.6241  290 

REGRESSION  446.8484  6  74.47472 

RESIDUAL  455.7758  284  1.604844 

R  SGUAREC  =  0.4950548  R  =  0.7026013 

THE  STANDARD  ERROR  OF  Y  ABOUT  THE  REGRESSION  PLANE  IS  1.266824 

THE  F  RATIO  TESTING  THE  NULL  HYPOTHESIS  THAT  ALL  THE 
COEFFICENTS  ARE  2ERC»  I.L.  B 1  =  B2  =  .  .  . =BK  =  C «  IS  46.40620 
EASED  ON  THE  DEGREES  OF  FREEDOM  6/  264 


COEF  #  BEST  VALUE  CORRELATION  STD. 
r  4.806970  0.7664793E-01 

1  0.3756840  0 . 2 4 4 0 5 1  2 E  -  0 2 

2  -0.1188613  0 .1665213E-03 

3  -0 .2323101E-01  0 . 1 2 2 8 28 2 E -04 

4  0.7416429F-02  0 . 1 7 1 3 925E - 0 5 

7  0 .3276729E- 03  0 .33 096 0 2E - 08 


ERROR  T 

0.3506948 

C.743C031E-01 

0.1624610E-01 

0.4435439E-02 

0.1658345E-02 

C.72873C5E-04 


13.70699 

5.056291 

-7.271542 

-5.232871 

4.472186 

4.496490 


ANALYSIS  OF  VARIANCE 

SOURCE  OF  VARIATION  SUM  CF  SQUARES 
TOTAL  9C2.6241 

REGRESSION  445.3221 

RESIOUAL  457.3019 


DCF  MEAN  SGUARE 

290 

5  89.06441 

285  1.604568 


R  SQUARED 


0.4933639 


R 


0.7022986 


THE  STANDARD  ERROR  OF  Y  ABOUT  THE  REGRESSION  PLANE  IS  1.266716 

THE  F  RATIO  TESTING  THE  NULL  HYPOTHESIS  THAT  ALL  THE 
COEFFICENTS  ARE  2ERO.  I.E.  B 1  =  B2  =  .  .  .=BK  =  0  •  IS  55.50678 
BASED  ON  THE  DEGREES  OF  FREEDOM  5/  285 


CCEF  #  EEST  VALUE  CORRELATION  STD.  ERROR  T 

0  A .425452  0 . 7 2 0 6 A  7 0E- 0  1  0.3511621 

1  0.4161844  0 .3389400E-02  0 .76 15663E -0 1 

2  -0 .4793454E-01  0 . 9 76 59 1 2E - 05  0 .40 8 7 92 3 E - C 2 

3  -0 .2572381F-C1  0 . 1 2 C 89 1 8 E - 0 4  C  .454P255E-Q2 

7  0 .3682456E-03  0 .3258307E-08  0 .7466934E-04 


ANALYSIS  OF  VARIANCE 

SOURCE  OF  VARIATION  SUM  OF  SGUARFS  CCF  FEAN  SGUARE 
TOTAL  902.6241  290 

REGRESSION  413.2301  4  103.3075 

RESIDUAL  489.3940  286  1.711168 

R  SGUAREC  =  C.4578C97  R  =  0.6766163 

THE  STANDARD  ERROR  OF  Y  ABOUT  THE  REGRESSION  PLANE  IS  1.30F116 

THE  F  RATIO  TESTING  THE  NULL  HYPOTHESIS  THAT  ALL  THE 
COEFFICENTS  ARE  2ERC,  I.E.  E 1 =B2  =  .  .  .  =BK = 0  ,  IS  60.37251 
BASEC  ON  THE  DEGREES  OF  FREECOM  4/  286 


COE F  #  BEST  VALUE  CORRELATION  STD.  ERROR  T 

:  5.742243  0 . 3 02 75 0 7E- 01  C. 2366758  24.26206 

1  0.7725544E-G1  0 . 6 29 24 0 2E- C 3  C  .34 1 2 0 8 3E- C 1  2.264172 

2  -0.5168957E-01  0 .942 711 2E- 05  0 .41 76386E-02  -12.37663 

3  -0 .3784753E-02  0 .5239919E-06  0 .9846312E-03  -3.843828 


ANALYSIS  OF  VARIANCE 

SOURCE  OF  VARIATION  SUM  OF  SQUARES  DOF  MEAN  SGUARE 
TOTAL  9  C  2 .624 1  290 

REGRESSION  371.6119  3  123.8706 

RESIDUAL  531.0122  287  1.850217 

R  SQUAREO  =  0.4117017  R  =  0.6416398 

THE  STANDARD  ERROR  OF  Y  ABOUT  THE  REGRESSION  PLANE  IS  1.36C227 

THE  F  RATIO  TESTING  THE  NULL  HYPOTHESIS  THAT  ALL  THE 
COEFFICENTS  ARE  ZERO.  I.E.  B  1  =  B2  = .  •  .=BK  =  0 «  IS  66.94923 
BASED  ON  THE  DEGREES  OF  FREEDOM  3/  287 


COEF  ft  BEST  VALUE  CORRELATION  STD.  ERROR  T 

0  6.227016  0 .5498795E- 02  C. 1015861  61.29792 

2  -0.5185799E-01  0 .9 424 1 2  IE- 05  0 .42 05532E -0 2  -12.33090 

3  -0 .1692082E-02  0 .6229162E-07  C .34 19 128E-03  -4.948870 


ANALYSIS  OF  VARIANCE 

SOURCE  OF  VARIATION  SUM  OF  SQUARES  DOF  MEAN  SQUARE 
TOTAL  902.6241  290 

REGRESSION  362.1267  2  181.0634 

RESIDUAL  540.4973  288  1.876727 

R  SQUARED  =  0.4011933  R  =  0.6333982 

THE  STANDARD  ERROR  OF  Y  ABOUT  THE  REGRESSION  PLANE  IS  1.365537 

THE  F  RATIO  TESTING  THE  NULL  HYPOTHESIS  THAT  ALL  THE 
COEFFICENTS  ARE  ZERO.  I.E.  B 1  =  B2  =  .  .  .  =  BK  =  0 «  IS  56.47827 
BASED  ON  THE  DEGREES  OF  FREEOOM  2/  288 


COEF  «  BEST  VALUE  CORRELATION  STD.  ERROR  T 

C  5.983677  0 .421 05 07E- 02  C.9243528E-01  64.73369 


12.59661 

5.464848 

11.72589 

5.655754 

4.931684 


63 


0.5422976E-01  0 . 9 3 0 1 7 3 5E - C 5  C  .4344626F.-G2  -12.48203 


e. 


ANALYSIS  OF  VARIANCE 

SOURCE  OF  VARIATION  SUF  CF  SQUARES  DCF  FEAN  SCUARE 
TOTAL  9C2.6241  290 

REGRESSION  318.1632  1  316.1632 

RESIDUAL  586  . 4608  289  2  .C29276 

R  SGUARED  =  0.35C2712  R  =  0.5918372 

THE  STANDARD  ERROR  OF  Y  AECUT  THE  REGRESSION  FLAKE  IS  1.424527 

THE  F  RATIO  TESTING  THE  NULL  HYPOTHESIS  THAT  ALL  THF 
COEFFICENTS  ARE  ZERO*  I.E.  B 1  =  B2  =  .  .  .  =  BK  -  C  ,  IS  155.8010 
EASED  ON  THE  DEGREES  OF  FREECOM  1/  289 


A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 


Mellor,  Malcolm 

Breaking  ice  with  explosives  /  by  Malcolm  Mellor. 
Hanover,  N.H. :  Cold  Regions  Research  and  Engineering 
Laboratory;  Springfield,  Va.:  available  from  National 
Technical  Information  Service,  1982. 

iv,  68  p.,  illus. ;  28  cm.  (  CRREL  Report  82-40.  ) 

Prepared  for  Office  of  the  Chief  of  Engineers  by 
Corps  of  Engineers,  U.S.  Army  Cold  Regions  Research 
and  Engineering  Laboratory  under  DA  Project  4A762730 
AT42. 

Bibliography:  p.  28. 

1.  Ice.  2.  Explosion  bubbles.  3.  Explosion  ef¬ 
fects.  4.  Explosions.  5.  Underwater  explosions. 
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